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Abstract 

For since the nature of our intellect is to abstract the essence of material things from 
matter, anything material residing in that abstracted essence can again be made 
subject to abstraction; and as the process of abstraction cannot go on forever, it must 
arrive at length at some immaterial essence, absolutely without matter; and this 
would be the understanding of immaterial substance. 

St. Thomas Aquinas, Summa Theologica, I, 88, 2. 

We address in this thesis two primary questions aimed at improving our ability to calculate 
reliably in the Standard Model of particle physics and probing possible new particles which 
may exist beyond it. 

First, we embark on an attempt to account for the abundance of matter in the present 
Universe if earlier in its history matter and antimatter were equally abundant. We explore 
whether baryogencsis at the elcctroweak phase transition could successfully account for the 
observed density of baryons in the Universe, using the closed-time-path (CTP) formalism 
of quantum field theory to calculate the buildup and relaxation of particle densities during 
the phase transition. For our model of the new particles and sources of CP violation 
necessary to account for the baryon asymmetry of the Universe, we adopt the Minimal 
Supersymmetric Extension of the Standard Model (MSSM). We look for regions of the 
parameter space in the MSSM that could give rise to sufficiently large baryon asymmetry 
without violating constraints on these parameters from existing experiments, in particular, 
constraints on masses of Higgs and supersymmetric particles from accelerator searches and 
precision electroweak tests, and on CP- violating parameters of the MSSM from searches 
for electric dipole moments of elementary particles. 

Next, we explore how to get around our ignorance of the dynamics of strongly interact- 
ing particles in the nonperturbative regime of Quantum Chromodynamics (QCD) by the 
clever use of effective field theories. Two applications are explored: the decay of Z bosons 
to hadronic jets using soft-collinear effective theory (SCET) and the radiative decays of 
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quarkonia to light hadrons using SCET and non-relativistic QCD (NRQCD). These tools 
facilitate the proof of factorization of decay rates into perturbatively-calculable and nonper- 
turbative parts. Universality of the latter among different observables provides predictive 
power even in our ignorance of the details of the nonperturbative physics. 
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Chapter 1 

Introduction 

Staunen kann nur, wer noch nicht das Ganze sieht; Gott staunt nicht. 

Josef Pieper, Gliick und Kontemplation 

Amazement is only possible for one who does not yet see the whole; God cannot be 
amazed. 

Josef Pieper, Happiness & Contemplation 

Today physicists eagerly anticipate the discovery of new phenomena beyond those already 
discovered in the Standard Model of particle physics. We await with high hopes the un- 
veiling of new particles and phenomena at the Large Hadron Collider, or perhaps even the 
Tevatron, within the decade. Meanwhile, we aim our current efforts at the improvement of 
the precision and reliability of both theoretical calculations and experimental tests of the 
properties of particles, not only those currently within the realm of speculation, but also 
of those within the Standard Model itself. On the one hand, to be able to extract evi- 
dence for new physics from future accelerator data, we must first calculate the predictions 
of the Standard Model as accurately and precisely as possible, in order to be able to find 
in the data deviations therefrom. On the other hand, we can already search for signatures 
of new physics in experimental data today, in high-precision low-energy experiments or in 
cosmological observations. 

In this thesis, we focus both on tests of models of new physics and improvement of our 
ability to calculate reliably in the Standard Model. First, we consider how extra sources 
of CP violation beyond those in the Standard Model (for example, in its supersymmetric 
extensions) might account for the baryon asymmetry of the Universe and how experimental 
tests, especially searches for electric dipole moments of elementary particles, constrain such 
possibilities. Then, backing up to the Standard Model, we attempt to improve calculations 
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of observables involving the strong interactions by use of effective field theories of Quantum 
Chromodynamics, studying the radiative decays of heavy quarkonia and the hadronic decays 
of Z bosons, attempting to understand the nonperturbative contributions to observables in 
these processes. Understanding hadronic jet production will be vital to the separation of 
QCD backgrounds at the LHC from signals of new physics, while heavy quarkonia provide 
a useful testing ground for the validity of the effective theories we use, establishing their 
reliability for other applications. 

1.1 Probing New Physics with Baryogenesis and Electric 
Dipole Moments 

Already today, suggestive evidence for physics beyond the Standard Model exists outside of 
accelerators, in the baryon asymmetry of the Universe (BAU) — the survival of more matter 
than antimatter in the Universe's early evolution. Big Bang Nucleosynthesis (BBN) and 
measurements of the cosmic microwave background tell us that: 



TlB 



3.4-6.9x10"^°, BBN m 

(1.1) 

5.9 - 7.3 X 10-1°, CMB 



where the number density of baryons is normalized to the number density of photons n^. 
(An updated version of these constraints from Ref. [3] is displayed graphically in Fig. 11.11 ) 
The Standard Model is unable to account for the (large!) size of these numbers, assuming 
the Universe begins in a symmetric state ((-B) = 0). Even if we assume a primordial excess 
of matter at the Big Bang, we could not easily account for its survival during the subsequent 
inflation, which would have diluted it nearly to zero, and which we must postulate to account 
for the flatness and uniformity of the Universe ^ and the density perturbations^ therein. 
To start in a matter-antimatter symmetric state and generate a baryon asymmetry requires, 
as pointed out by Sakharov j^, three conditions: 

• Baryon number violation — otherwise, the Universe would never depart from (B) = 0. 

• C and CP violation — otherwise, even with B violation, particles and their antiparti- 
cles are produced in equal amounts. 



^Which, of course, can be augmented after inflation too. 
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Figure 1.1: BBN and CMB constraints on the BAU, as analyzed in Ref. 0. The abundances 
of hehum-4 (Y), deuterium, hehum-3, and hthium-7 can be used to deduce the baryon 
density 0,b, or, equivalently, the baryon-to-photon ratio rj. The various curves are the BBN 
predictions for these abundances as functions of the baryon density, which are constrained 
by observations to he in the outhned boxes. The combination of the constraints from BBN 
and from WMAP measurements of the CMB is shown as the yellow, vertical band. (Figure 
courtesy of R. Cyburt.) 
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• Departure from thermal equilibrium — otherwise, the thermal average 



Tr(e-H/T^) 



never departs from zero, since 



TV [e-f^^B] = Tr [e-'^^{CPT){CPT)-^B] 

= Tr [{CPT)e-'^^{CPT)-^B] 

= -TV [e-'^^{CPT)-\CPT)B] = -TV [e-^^B] 
TV [e-f^^B] = 0, 



:i.3) 



assuming CPT invariance^ {\H,CPT\ = 0), and using cyclicity of the trace and the 
oddness of B under a CPT transformation. 

These three conditions may have been realized in electroweak baryogenesis^ , in which the 
BAU is generated during the electroweak phase transition — if it is strongly-enough first- 
order and there is enough CP violation in the electroweak sector of Nature. Neither of 
these conditions is satisfied in the Standard Model. The phase transition is not found to be 
strongly-enough first-order given present limits on the Higgs mass, and the complex phase 
in the CKM matrix is too small to account for the BAU. 

Before considering how to overcome the inability of the Standard Model to satisfy the 
last two Sakharov criteria, let us begin by examining how it satisfies the first. 



1.1.1 Baryon Number Violation 

Baryon number violation is already present in the Standard Model. Although B is conserved 
in all perturbative processes, there are nonperturbative topological transitions of gauge fields 
at high temperature which generate baryon and lepton number violation through triangle 
anomalies ^Hl- The baryon number current, 

fB = j-Y.'it^''i?' (1-4) 

i,a 

^We do not consider the possibility that CPT invariance may be broken by, for instance, Lorentz violation, 
perhaps spontaneously |51I^ I10| . Mechanisms violating CPT spontaneously have also been used to account 
for baryogenesis lllj . 

^This and other scenarios for baryogenesis are reviewed in Ref. |12|. 



Figure 1.2: Baryon number current anomaly in SU{2). 



summed over the np quark flavors i and Nq colors a, is not conserved at the quantum level 
in the Standard Model, due to the diagram shown in Fig. 11.21 

d,fB = (1-5) 

where W^^'^ is the SU{2) gauge field strength tensor, and W^^'^ = |e'^''"^Wa/3- The right 
hand side of this equation is minus the divergence of a topological current, —d^K>^, where 
[HI: 

= iLe'^-Z^TV (w,d^W/s - ^W.W^W^) , (1.6) 

where is the SU (2) gauge field. Integrating d^K^ over all space, we obtain the topo- 
logical charge, 

ncs = -jd''^K\ (1.7) 

called the Chern-Simons number. Different vacuum configurations of gauge fields having 
the same energy can have different values of ncs- These vacua correspond to gauge field 
configurations of the form: 

= -nd^n-\ (1.8) 

where i7 is a function from four-dimensional space into the gauge group SU{2). As SU{2) 
is topologically equivalent to the three-sphere 5^, the index ncs is essentially the winding 
number of the map Q. onto [in]. 

Transitions between these topologically distinct vacua induce a change in the baryon 
number, due to the anomaly equation H1.5() . as illustrated in Fig. 11.31 The lepton number 
current has the same anomaly. Thus, i? + L is violated while B — L \s conserved. These 
vacua are separated by barriers corresponding to solutions of the field equations called 
sphaleron configuratoins with energy Egph- The rate of tunneling between vacua is negli- 
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Figure 1.3: Topological transitions at high temperature. When SU{2) gauge fields tunnel 
from one vacuum configuration to another, characterized by different values of the Chern- 
Simons number, baryon and lepton number are violated through the quantum anomaly 
in the baryon and lepton number currents. The unstable solutions to the field equations 
between the vacua are called sphaleron configurations, whose energy is Sgph- The rate of 
the tunneling reactions is Fsph- 

gible below the temperature of the electroweak phase transition, the rate suprressed by a 
factor of exp(— Svr^/^l) ~ 10^^'^^ 16 . However, at higher temperatures, the tunneling is 
unsuppressed, and the sphaleron transitions induce a change in baryon number through the 
equation: 

dtPB{x) - DqV'^pBix) oc -npTsph[nL{x) + Rpb{x)], (1.9) 

where Dq is the diffusion coefficient for baryons, np is the number of families, Tgph is the 
rate of weak sphaleron transitions, R is a relaxation coefficient for baryon number, and np 
is the number density of left-handed weak-doublet fermions. The latter quantity appears 
because SU (2) gauge bosons couple directly only to left-handed fermions. 

We imagine a picture in which bubbles of the broken electroweak symmetry phase nu- 
cleate in the previously electroweak symmetric Universe. At or near the bubble boundaries, 
the generation of a nonzero due to CP violation in regions where weak sphaleron transi- 
tions are active leads to a nonzero baryon density pB, which, if it ends up inside the region 
of broken electroweak phase where sphalerons no longer change baryon number, is frozen 
in and survives until the present day in the Universe. The goal within this framework is to 
derive an equation for np dependent on the amount of CP violation present and solve for 
the final baryon density. The generation of n^, which tends to occur much faster than the 
sphaleron transitions themselves, is governed by another set of equations, involving both 
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CP-violating sources generating nonzero and CP-conserving reactions which tend to 
cause relaxation of n/, back towards zero. The presence of a nonzero ul then drives the 
weak sphaleron transitions to generate nonzero baryon number. 

We will consider the equations for tll in more detail after examining where we might 
find enough CP violation beyond the Standard Model to account for the observed size of 
the BAU. 

1.1.2 More CP Violation from Supersymmetry 

Supersymmetric extensions of the Standard Model^ are popular for their touted ability to 
resolve the hierachy problem (by canceling large corrections to the Higgs mass from heavy 
particle loops), achieve gauge coupling unification around 10^^ GeV, and provide candidate 
particles for the dark matter. What interests us is that they may provide a whole new 
set of CP violating parameters that may help account for the BAU. These may generate 
large enough and, thus, through Eq. H1.9|) . pB- Indeed, the Minimal Supersymmetric 
Extension of the Standard model (MSSM) by itself contains more than 100 free parameters 
including dozens of CP-violating phases. Adopting a specific model for supersymmetry 
breaking can reduce this menagerie to a more tractable set. In the minimal supergravity 
model, for instance, which we adopt in thesis, the only independent CP-violating phases 
remaining are the phase t/)^ of the /^-parameter in the superpotential: 

WussM D ^i{H+H;, - HlH^^) + • • • , (1.10) 

where p = l/ile*"^^, and the phases in the triscalar couplings in the soft SUSY-breaking 
terms in the Lagrangian: 

£soft ^ - {ua-uQIiu - da.dQHd - EaeLHd^ + c.c, (l-H) 

where the matrices au = AuoYu, = ^doYd, and ag = A^oye are proportional to the 
corresponding matrices of Yukawa couplings. Assuming that Auo = Ado = ^eo = at the 
scale of SUSY breaking, there is just one new CP-violating phase, (pA, where Aq = \ Aq\ e*"^^. 
These two CP-violating phases (j)^ and 0^ may be large enough to account for the baryon 
''We refer to Ref. 23 for details and notation. 
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asymmetry. They cannot, however, be so large as to conflict with constraints from precision 
tests searching for permanent electric dipole moments (EDMs) in electrons, neutrons, and 
atomic systems. These can possess an EDM only if time-reversal (T) symmetry, and, 
therefore, by the CPT theorem, CP symmetry, is violated in Nature. So far, no experiment 
has found a measurable EDM, thus placing stringent constraints on the size of any new 
sources of CP-violation beyond the Standard Model. 

Our task is to calculate the baryon density that could be generated from CP-violating 
phases small enough not to conflict with EDM constraints. 

1.1.3 Transport Equations for Electroweak Baryogenesis 

With this particular model for particles which can participate in the generation of the 
baryon asymmetry at the electroweak phase transition, we can write equations governing the 
generation of left-handed weak doublet fermions rii entering Eq. ()1.9() for pB, SiS introduced 
by Cohen, Kaplan, and Nelson ^HI and Huet and Nelson For example, the number 

density of left-handed third-generation quarks and squarks Q is governed by an equation of 
the form: 

a.Q - D,v^Q . r„ (X _ ^) + (£ - £ - ^) + 5f + . . . . ,1.12) 

where T is the density of right-handed third- generation quarks and squarks, H is the density 
of Higgs and Higgsinos, kx^Q^H are statistical factors, Sq is the CP-violating source for 
left-handed squarks, and T^^^y are the relaxation rates for the quantities in parentheses 
induced by interactions of quarks and squarks with the Higgs vacuum expectation value 
and with real Higgs particles, respectively. In this framework, supergauge interactions are 
assumed to be much faster, keeping superpartners in chemical equilibrium with each other, 
accounting for our use of combined particle and sparticle densities T, Q, and H. Equations 
similar to Eq. H1.12() for T and H round out a set of coupled equations whose solution gives 
the density ni which enters the equation for baryon density ()1.9() . 

In Refs. |181ll9j semi-classical techniques were adopted to estimate the size of the coeffi- 
cients like Tm,y and sources S^. Riotto Uni introduced into this framework the closed time 
path (CTP) formulation of quantum field theory, which incorporates finite-temperature 
and nonequilbrium effects, to derive from quantum field theory the sources but kept 



9 

the standard derivations of ^m,y He discovered large enhancements (of the order of 10^) 
over previous calculations of in certain regions of MSSM parameter space, namely, 
where certain particle masses (such as Higgsinos and Winos or left- and right-handed stops) 
become degenerate. 

In Ref. [2J and Chap. |21of this thesis, we continue Riotto's endeavor by deriving other 
terms appearing in the diffusion equations like (|1.12|) using the same CTP formalism as 
for the CP- violating sources. We discover similar enhancements in the coefficient Tm, and 
expect similar results for Ty, which has not yet been calculated with the CTP formal- 
ism. Enhancement of these rates tends to reduce the final baryon asymmetry compared to 
Riotto's result, but still leaving a significant enhancement oi ps- Beyond these enhance- 
ments, we also discover entirely new terms in Eq. ()1.12|) and its cousins with different linear 
combinations of the various particle densities. 

These findings, although still quite preliminary, make evident the importance of a fully 
consistent, quantum- field-theoretic calculation within the CTP formalism of all the terms 
entering the transport equations for particle densities. Each new application of this formal- 
ism to a different set of terms in the transport equations (e.g., by Riotto and by us) has 
significantly affected the size of the BAU predicted by these equations. The prospect of 
more powerful and more precise experiments that will constrain the parameters in models 
of new physics in the next several years make imperative the task of making as complete 
and reliable as possible these theoretical predictions. The work in this thesis, focusing at- 
tention just on a few of the terms in the transport equations within only the MSSM, lays 
the groundwork for future calculations, which ought also to extend to other extensions of 
the SM to prepare for the various possibilities for the new physics that may be discovered 
within the decade. 

1.1.4 Combining Constraints from BAU and EDMs 

In Fig. 11.41 we illustrate the present status of BAU and EDM constraints on (j)^ and (pA 
in the MSSM within the mSUGRA scenario for supersymmetry breaking, for a particular 
choice in the MSSM parameter space. Our choices are informed by precision electroweak 
constraints [121 combined with numerical simulations of the strength of the electroweak 
phase transition [2S1, which tell us that, in the MSSM, we must have a light right-handed 
top squark tR and heavy left-handed tL- We also choose degenerate Higgsino and SU{2)- 
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Figure 1.4: Combined constraints on CP-violating phases from baryon asymmetry and 
electric dipole moments. The colored bands pick out the region in the (j)fj,-(pA plane required 
for successful electroweak baryogenesis, based on WMAP (green) and BBN (green+blue) 
observations of the BAU. The wider white bands are limits placed on the size of these phases 
by electron EDM (de) and neutron EDM (d^) searches. There exists a region of overlap 
where BAU requirements are consistent with current EDM constraints. The narrow lines 
inside those bands denote the sensitivity of proposed searches for these EDMs at future 
experiments at Los Alamos National Laboratory [241 125j and Yale |26| . 

gaugino masses |^| = M2, which we find maximizes the BAU that can be generated from 
the phases (/>^, (pA of a given size (here we chose \fi\ = M2 = 200 GeV). 

We are on the verge of new EDM experiments in the next few years which will achieve 
sensitivites of two or more orders of magnitude greater than the present limits illustrated 
in Fig. 11.41 Further null results may succeed in ruling out the simplest scenarios of super- 
symmetric electroweak baryogenesis entirely. Indeed, achieving the situation illustrated in 
Fig. 11.41 already places the fairly stringent constraint of near-degeneracy between Higgsino 
and gaugino masses to generate a large enough BAU from the small phases implied by the 
EDM searches. Meanwhile, with the LEP-II constraint on the Higgs mass of mn ^114 GeV 
j27j . only a small window remains for a strongly first-order electroweak phase transition, 
requiring mn ^ 120 GeV (see Ref. |2H] and references therein). These constraints, however, 
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become relaxed with simple extensions of the MSSM^, to which we would turn if future ex- 
periments rule out the minimal scenario. More optimistically, discoveries of nonzero EDMs 
may pinpoint masses of superpartners that are indeed consistent with the BAU generated 
through electroweak baryogenesis (by picking out phases within the colored band in Fig. 11.41 
for instance), even before their possible discovery at LHC. 

All the data combined together would provide a powerful test of the scenario of elec- 
troweak baryogenesis in the MSSM, ruling it out for good (and leading us to consider 
extensions of the MSSM or other models of new physics) or perhaps even confirming it in 
a dramatic way. 

1.2 Effective Field Theories for Strong Interactions 

Before we can reliably analyze the data from experiments like those at the LHC, we have 
much work to do establish reliable calculations within the Standard Model itself, especially 
in the sector of strong interactions. The theory describing these interactions. Quantum 
Chromodynamics (QCD), possesses the important property of asymptotic freedom which 
makes the interactions between quarks at high energies amenable to relatively simple per- 
turbative theoretical calculations. However, at low energies, where quarks interact strongly 
and bind together to make hadrons, calculations are almost impossible, except numerically 
using lattice QCD. 

1.2.1 Nonperturbative Effects in Hadronic Jets 

One way to make progress analytically in the nonperturbative domain of QCD is to approx- 
imate the full theory with an effective field theory (EFT) which simplifies the separation 
of perturbative and nonperturbative effects. Given a particular class of physical processes 
to study, we can identify the relevant sector of the full theory which describes them. By 
integrating out degrees of freedom in the full theory which do not propagate for long times 
or distances and identifying small parameters characterizing the physics in which we can 
expand the full theory, we can form the appropriate EFT, which often possesses, at a 

■''For example, Ref. |29| found that in the NMSSM (Next-to-Minimal Supersymmetric Extension of the 
Standard Model) — which contains an extra gauge singlet chiral superfield — a strongly first-order electroweak 
phase transition can be achieved with heavy squarks and Higgs masses heavier than the MSSM limit of 120 
GeV. 
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given level of approximation, extra symmetries or other technical features facilitating the 
separation of perturbative and nonperturbative contributions to physical observables. 



One application of this strategy is to event shape variables in e~^e~ annihilation or Z 
decay to hadrons. Most events in these processes produce two or more jets of hadrons 
arising from the underlying partonic sub-processes. A number of variables conventionally 
used characterize the "jettiness" of an event. For example, in Z decay, the energy of a 
produced jet Ej is very close to Mz/2 for an event with two back-to-back jets. Another 
popular variable is the thrust: 



the sum of the projections of the three-momenta of all particles i in the final state onto 
an axis t, chosen to be the axis which maximizes this sum. (For exactly back-to-back jets, 
T = 1.) The distribution of events in thrust can be calculated in perturbation theory 
at the parton level, but hadronization effects will introduce additional nonperturbative 
contributions to this observable. In Fig. 11.51 taken from Ref. ^HU^, we find a comparison of 
a prediction of the thrust distribution da/dT in e^e^ annihilation in perturbation theory 
(the red dotted line) to experimental data from LEP-II. The prediction does not quite 
match the data, especially near the kinematic endpoint 1 — T ~ where two-jet-like events 
are found, until nonperturbative contributions are included, giving the black solid line. 
Of course, the nonperturbative contributions are not calculable; they are modeled using 
unknown parameters, which are merely fitted to the data. It would seem that no predictive 
power is truly gained. 

However, the nonperturbative contribution to these observables can often be expressed 
as the matrix element of some operator in QCD or an effective theory which approximates 
it, giving us some possibly useful information. For example, the total decay rate of Z to 
hadrons is given by: 




(1.13) 



r{Z ^ hadrons) = -^Yl J 



d^xe'P^''{0\J''{x)r*{0)\0) 



(1.14) 
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Figure 1.5: Thrust distribution in e+e~ annihilation. The dotted curve is a perturbative 
prediction for the thrust distribution, while the solid curve is a fit to the data in a model 
including the nonperturbative power corrections. (Figure from Ref. |3U).) 

summed over Z polarizations e, and where the current mediating the Z decay is: 

J^ = Y. ^^avr + 5A7^5)gf . (1.15) 
Using the operator product expansion, we find that the product of currents can be expressed: 

J^(x)j;(0) =C;i,(x)T¥l + C^;(x)TVG,^G"^ + -- - , (1.16) 

~ M^, C^^ ~ M'^, and the higher-order terms are suppressed by more powers of 
The first term gives rise to the purely perturbative prediction for the decay rate, while the 
second encodes the leading-order nonperturbative physics, in the matrix element: 

(0|TrG„^G"^|0)~AtjcD- (1-17) 

Thus the nonperturbative contribution to the total hadronic decay rate of the Z is sup- 
pressed relative to the perturbative contribution by a factor of order ~ (Aqcd/^z)^, or 
®This power counting would, of course, differ in six spacetime dimensions |31|. 



14 

~ 10~^. This means that the perturbative calculation should be highly accurate; on the 
other hand, it would be extremely difficult to extract the value of the matrix element p.l7jl 
from a measurement oiT(Z ^ hadrons). 

Two features of an observable would enhance our predictive power. First, its leading 
nonperturbative contribution should be large enough that its effect would be measurable in 
an experiment while being small enough not to overwhelm the perturbative contribution. 
Second, the nonperturbative contribution should be characterized as the matrix element of 
operator which also contributes to a different physical observable. This way, its extracted 
value in one experiment could be used to predict the results of another. 

The attempt to obtain these desirable features is the goal of our use of effective field 
theory to analyze event shape variables in hadronic Z decay. We use an EFT designed to de- 
scribe those processes involving strongly-interacting particles moving with large energy com- 
pared to their invariant mass, as in the hadronic jets produced in these events. In the soft- 
collinear effective theory (SCET), we expand QCD in a small parameter A ~ ^/^QCD/^z, 
the typical transverse momenta of particles inside a jet of hadrons, and keep only lightlike 
(collinear) and soft degrees of freedom. At leading order in A, the direct couplings in the 
Lagrangian between soft and collinear particles can be made to disappear in the effective 
theory by clever field redefinitions. This facilitates immensely the proof of factorization of 
event shape variables into hard and soft contributions. The latter are expressed as matrix 
elements of operators in the effective theory, and so we embark on the search for such matrix 
elements which are universal among different event shape variables. 

In Ref. IHSj, Bauer, Manohar, and Wise showed, at leading order in SCET, that the en- 
ergy Ej of a single observed jet depends on no more than two nonperturbative parameters, 
regardless of the number of jets in the event. In Ref. and in this thesis, the non- 
perturbative contributions to other event shape variables are similarly analyzed (and the 
calculation in Ref. [32] of the jet energy distribution in two-jet events is extended to 0{as) 
in perturbation theory.) Simple relations for different variables were proposed by Dokshitzer 
and Webber in Ref. [251 • Our analysis suggests, unfortunately, that only two of the several 
commonly-used event shape variables — thrust and jet mass sum (see Sec. I5.3.2|) — receive 
the same nonperturbative contributions, while those for the other variables are unrelated. 
The Dokshitzer- Webber model was tested by the DELPHI collaboration |3H]j whose results 
are summarized in Fig. 15.41 and appear to be consistent with our claims. 



15 

1.2.2 Reliable Predictions for Radiative Upsilon Decay 



Similar strategies can be pursued to analyze particular exclusive decays instead of inclusive 
or semi-inclusive processes like Z decays to hadrons. The appearance of universal nonper- 
turbative quantities in different exclusive decays would allow us to predict the branching 
ratios for particular decay channels once the nonperturbative contributions are measured 
in any one of them. 

In this thesis we apply SCET and non-relativistic QCD (NRQCD) to decays of the T 
meson to a photon and one or more light hadrons. While SCET is appropriate to describe 
the final hadronic state, NRQCD is required to describe the heavy quark-antiquark pair 
inside the T. 

In NRQCD, the T decay rate can be split into several different contributions arising from 
the possible spin and color configurations of the inital 66 pair. To leading approximation, the 
66 can be considered to be in a color-singlet, spin-triplet configuration which dominates 
the total T decay rate. For example, the photon energy spectrum in inclusive radiative T 
decay depends on the sum of matrix elements: 

(T ^ JX) r^^CnirPn |T) , (1.18) 



for operators On with various spin and color quantum numbers. In the approximation of 
the color singlet model, only one operator contributes: 

^(T ^ jX) ~ (TlxlpCrVp • xlp^V'p |T) . (1.19) 
This approximation is inadequate for some regions of the photon energy spectrum. Roth- 



66(8,1^0) I 66(1,3^1) 

— ^ — Uoooooa^ 

Figure 1.6: Color-octet and color-singlet channels in radiative T decay. The photon energy 
is more peaked near = My/ 2 in the color-octet channel and may need to be included 
for a reliable prediction of the decay rate in this kinematic region. 
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Figure 1.7: Photon spectrum in inclusive radiative T decay. This plot of the number of 
events with a photon energy E^^ = Myz/l from Ref. |SZ| compares data from CLEO 
j38j (the boxes) with predictions including only the color singlet contribution [39j (the blue 
dashed line) and including both color singlet and octet contributions j37j (the one solid 
and two dot-dashed lines). These latter three lines correspond to different choices for the 
renormalization scale /i. The predictions which include the color octet contributions have 
significantly improved agreement with data. 

stein and Wise [301 pointed out that color octet channels, for example, 

(TlxLpT^xV^IT), (1.20) 

could be equally important as the color singlet near the kinematic endpoint ~ Mx/2. 
As illustrated in Fig. 11.61 while the color-singlet configuration must decay to at least two 
gluons, the color-octet can decay to just one. The photon energy will be peaked closer to 
Afx/2 when recoiling against one gluon instead of two. 

In Fig. 11.71 we see that the color-singlet contribution is inadequate to account for the 
experimentally observed photon energy spectrum. Inclusion of the color-octet pieces are 
necessary to fit the data. Bauer, Fleming, et al. in Ref. jlUI and Fleming and Leibovich 
in Refs. ^Hl 1^ added SCET to the NRQCD analysis of radiative T decay, successfully 
introducing the collinear degrees of freedom necessary to analyze the dynamics near the 
kinematic endpoint E^^ = Afx/2 in the photon energy spectrum. 

In Ref. j42j and this thesis, we extend these analyses to the case of exclusive decays 
r(T ^ ^H) for a light hadron H. We show that the color-octet channel is in fact suppressed 
relative to the color-singlet even though we are in the region of large recoiling photon energy. 
The extra suppression comes from inclusion of interactions necessary to turn the single gluon 
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produced in the color-octet decay into a color-singlet hadron in the final state. Then, having 
thus dropped the color-octet contribution, we predict the sizes of the decay rates of the T 
to different hadrons in the final state. Working only to leading order in the effective theory 
and using its various symmetries, we deduce the operators which contribute dominantly to 
the radiative decays of the T and which hadrons those operators could produce in the final 
state. We find at leading order that the T should radiatively decay dominantly to the flavor- 
singlet, parity-even /2(1270). Furthermore, the appearance of universal nonperturbative 
parameters for various quarkonia decays to the /2 allows us to make predictions for the 
ratios of branching fractions such as 

5(T^7/2) _ 0.13 _ 0.18, (1.21) 



B{J/^ ^ 7/2) 

whereas experimentally this ratio is 0.06 it 0.03 1 . Although experimental and theoretical 
uncertainties are still fairly large to make robust comparisons of these ratios, no evidence 
has yet been found for any exclusive radiative decays of T to anything other than the /2. 
The observation of this channel before any others reassures us about the reliability of our 
prediction, whose further experimental scrutiny must await future data''. 



1.3 A Look Beyond 

New experiments promise an exciting era of particle physics in the coming decades, when 
we may verify current guesses about the nature of the physics that lies just beyond the 
Standard Model or flnd entirely other surprises. These prospects make imperative the im- 
provement of theoretical calculations in the Standard Model and more reliable calculations 
of the predictions of its extensions for baryogenesis, EDMs, and other such experimental 
observables. With luck, the LHC and EDM searches will begin producing relevant data 
within a few years. With hard work, the theoretical calculations will prepare scenarios 
such as supersymmetric electroweak baryogenesis for subjection to the rigorous and reliable 
scrutiny of those data. This thesis is an effort to make progress in these directions so as 
to be ready for any new, exciting physics that the next era of theory and experiments may 
reveal! 

^ Which appear to be forthcoming very shortly after the submission of this thesis |43|. Stay tuned! 
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1.4 Plan of the Thesis 



We begin our investigations beyond the Standard Model in Chap. |21 and calculate the 
baryon asymmetry of the Universe generated in electroweak baryogenesis in the MSSM and 
address the consistency of these results with constraints on CP violation in the MSSM from 
experimental searches for elementary particle EDMs. Much, though not all, of the material 
in this chapter was published in Ref. j21j . 

Then, returning to the Standard Model, in Chap. [31 we review the essential features of 
effective field theories in QCD, namely, soft-collinear effective theory and non-relativistic 
QCD, which will be applied to the radiative decays of T mesons in Chap. 0] (much of which 
was published in Ref. [l^]) and the hadronic decays of Z bosons in Chap. [5] (much of which 
was published in Ref. |33])) leading to the conclusion and future outlook in Chap.|BJ 
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Chapter 2 

Constraints on Supersymmetric 
Electroweak Baryogenesis 

In principio creavit Deus caelum et terrain. 

Terra autem erat inanis et vacua et tenebrae super faciem abyssi et spiritus Dei ferebatur 
super aquas. 

Genesis 1:1-2 

We begin our investigations beyond the Standard Model, considering the constraints that 
currently available data, especially the baryon density of the Universe and limits on electric 
dipole moments of elementary particles, might place on the parameters of the speculative 
but popular Minimal Supersymmetric Extension of the Standard Model, which we adopt 
as a theoretical testing ground for the application of the closed-time-path formalism of 
quantum field theory to the calculation of the baryon density generated in electroweak 
baryogenesis. (Much of this chapter appeared in Ref. ^21,.) 

2.1 Introduction 

The origin of the baryon asymmetry of the Universe (BAU) remains an important, un- 
solved problem for particle physics and cosmology. Assuming that the Universe was matter- 
antimatter symmetric at its birth, it is reasonable to suppose that interactions involving 
elementary particles generated the BAU during subsequent cosmological evolution. As 
noted by Sakharov [7], obtaining a nonzero BAU requires both a departure from thermal 
equilibrium as well as the breakdown of various discrete symmetries: baryon number {B) 
conservation, charge conjugation (C) invariance, and invariance under the combined C and 
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parity (P) transformations^. The Standard Model (SM) of strong and electroweak inter- 
actions satisfies tliese conditions and could, in principle, explain the observed size of the 
BAU: 



s 



(7.3±2.5) X 10-i\ BBN 



(2.1) 

(9.2± 1.1) X 10-i\ WMAP 



where pb the baryon number density, s is the entropy density of the universe, and 
where the values shown correspond to 95% confidence level results obtained from Big Bang 
Nucleosynthesis (BBN) P and the Wilkinson Microwave Anisotropy Probe (WMAP) [2j, 
respectively. In practice, however, neither the strength of the first-order electroweak phase 
transition in the SM nor the magnitude of SM C-P- violating interactions are sufficient to 
prevent washout of any net baryon number created by S-violating electroweak sphaleron 
transitions during the phase transition. 

The search for physics beyond the SM is motivated, in part, by the desire to find new 
particles whose interactions could overcome the failure of the SM to explain the BAU. 
From a phenomenological standpoint, a particularly attractive possibility is that masses 
of such particles are not too different from weak scale and that their interactions both 
strengthen the first-order electroweak phase transition and provide the requisite level of CP- 
violation needed for the BAU. Precision electroweak measurements as well as direct searches 
for new particles at the Tevatron and Large Hadron Collider may test this possibility, 
and experiment already provides rather stringent constraints on some of the most widely 
considered extensions of the SM. In the Minimal Supersymmetric Standard Model (MSSM), 
for example, present lower bounds on the mass of the lightest Higgs boson leave open 
only a small window for a sufficiently strong first-order phase transition, although this 
constraint may be relaxed by introducing new gauge degrees of freedom (see, e.g., [441 145j ) . 
Similarly, limits on the permanent electric dipole moments (EDMs) of elementary particles 
and atoms imply that the CP-violating phases in the MSSM must be unnaturally small (~ 
10^^). Whether such small phases (supersymmetric or otherwise) can provide for successful 
electroweak baryogenesis (EWB) has been an important consideration in past studies of 
this problem. 

In order to confront phenomenological constraints on the parameters of various elec- 



^ Allowing for a breakdown of CPT invariance relaxes the requirement of departure from thermal equi- 
librium. 
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troweak models with the requirements of EWB, one must describe the microscopic dynamics 
of the electroweak phase transition in a reahstic way. TheoreticaUy, the basic mechanism 
driving baryogenesis during the phase transition is well-estabhshed. Weak sphaleron tran- 
sitions that conserve B — L but change B and L individuahy are unsuppressed in regions 
of spacetime where electroweak symmetry is unbroken, while they become exponentially 
suppressed in regions of broken symmetry. Net baryon number is captured by expanding 
regions of broken symmetry ("bubbles"). Given sufficiently strong C and CP-violation as 
well as departure from thermal equilibrium, the non-zero B generated outside the bubble 
cannot be entirely washed out by elementary particle interactions that occur at the phase 
boundary. The baryon number density, ps, is governed by a diffusion equation of the form: 

dtPB{x) - DV^pb{x) = -T„sF„,{x)[nL{x) + Rpb{x)] , (2.2) 

where D is the diffusion coefficient for baryon number, F^s is the weak sphaleron transition 
rate, F^s{x) is a sphaleron transition profile function that goes to zero inside the regions 
of broken electroweak symmetry and asymptotically to unity outside, i? is a relaxation 
coefficient for the decay of baryon number through weak sphaleron transitions, and ni(x) 
is the number density of left-handed doublet fields created by "fast" chirality changing 
processes (see, e.g., US])- Thus, in order to obtain nonzero ps inside the bubble of broken 
electroweak symmetry, the left-handed density ni must be non- vanishing in the plasma at 
the phase boundary and possibly beyond into the region of unbroken symmetry. 

In effect, ni(x) acts as a seed for the i?-changing weak sphaleron transitions, and its 
spacetime profile is determined by the CP- violating sources and the quantum transport of 
various charges in the non-equilibrium environment of the plasma. Typical treatments of 
these dynamics involve writing down a set of coupled quantum transport equations (QTEs) 
for the relevant charges, estimating (or parameterizing) the relevant transport coefficients, 
and solving the system of equations under the appropriate boundary conditions. 

Among the developments in the past decade or so which have made significant impacts 
on this program, we identify two that form the basis of our investigation in this work. 
First, the authors of Ref. TH noted that diffusion of chiral charge ahead of the advancing 
phase transition boundary into the region of unbroken symmetry could enhance the impact 
of baryon number-changing sphaleron processes, thereby leading to more effective EWB. 
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The second, perhaps less widely-appreciated, development has been the observation by 
the author of Ref. [20] that the application of equilibrium quantum field theory (QFT) to 
transport properties in the plasma is not necessarily appropriate. In contrast to equilibrium 
quantum dynamics, the time evolution of quantum states during the phase transition is non- 
adiabatic. Consequently, scattering processes that drive quantum transport are no longer 
Markovian, but rather retain some memory of the system's quantum evolution. Using the 
closed time path (CTP) formulation of non-equilibrium QFT [17j to compute the CP- 
violating source terms in the plasma for the MSSM, the author of Ref. [20] found that 
these "memory effects" may lead to significant resonant enhancements (of order 10^) of 
the sources over their strength estimated in previous treatments (see, e.g., Ref. 19 and 
references therein) . The authors of Ref. |48l I49j subsequently found that performing an 
all-orders summation of scattering from Higgs backgrounds reduces the size of the CP- 
violating sources to some extent, but that the resonant enhancements nonetheless persist. 
Taken at face value, these enhancements would imply that successful EWB could occur 
with significantly smaller CP-violating phases than previously believed, thereby evading 
the present and prospective limits obtained from EDMs. 

To determine whether or not such conclusions are warranted, however, requires that one 
treat the other terms in the transport equations in the same manner as the CP- violating 
sources. Here, we attempt to do so, focusing on the terms that, in previous studies, have 
governed the relaxation of nL{x). In particular, chirality-changing Yukawa interactions 
with the Higgs fields and their spacetime varying vacuum expectation values (vevs) tend to 
wash out excess n/,(x). In earlier studies — including those in which non-equilibrium QFT 
has been applied to the CP-violating sources — these relaxation terms were estimated using 
conventional quantum transport theory jl91l2L)lll5ll49j . However, if the memory effects that 
enhance the CP- violating sources have a similar effect on these Yukawa terms, then the net 
effect on ps may not be as substantial as suggested in Refs. PUI 1151115] . 

The goal of the present study is to address this question by developing a more compre- 
hensive treatment of EWB using the CTP formulation of non-equilibrium QFT. In doing 
so, we follow the direction suggested in Ref. ^| and compute the transport coefficients of 
the chiral charges using the CTP formalism. To make the calculation more systematic, we 
identify the relevant energy and time scales that govern finite temperature, non-equilibrium 
dynamics and develop a power counting in the ratios of small to large scales (generically 
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denoted here as e). As we show below, both the CP-violating sources and the driving relax- 
ation terms first arise at 0{e^), and we truncate our analysis at this order. In contrast to 
the computation of the CP- violating sources, the derivation of the relaxation terms requires 
the use of finite density Green's functions. Given the resulting complexity, we consider here 
only the terms in the transport equations that previous authors have considered the domi- 
nant ones, and use our analysis of these terms to illustrate a method for obtaining a more 
comprehensive treatment of the QTEs. To make the phenomenological implications con- 
crete, we focus on the MSSM, realizing, however, that one may need to include extensions 
of the MSSM in order to satisfy the requirements of a strong first-order phase transition. 
Finally, we also attempt to identify the different approximations that have entered previ- 
ous treatments of EWB, such as the implicit truncation at a given order in e and outline 
additional calculations needed to obtain a comprehensive treatment. 

Based on our analysis, we find that under that same conditions that lead to resonant 
enhancements of the CP-violating sources, one also obtains a similar, resonant en- 

hancement of the driving chirality-changing transport coefficient, T. Since Yb ^ 
resonant relaxation counteracts the enhanced sources, though some overall enhancement 
of EWB still persists. Consequently, it will be important in future work to study the 
other transport coefficients whose impact has been considered sub-leading, since they may 
be enhanced under conditions other than those relevant for the leading terms. From the 
standpoint of phenomenology, we also illustrate how the implications of EDM searches for 
EWB depends in a detailed way on the electroweak model of interest as well as results from 
collider experiments and precision electroweak data. 

In presenting our study, we attempt to be somewhat pedagogical, since the methods are, 
perhaps, not generally familiar to either the practitioners of field theory or experimentalists. 
Most of the formal development appears in Sections I2.2H2.41 In Section 12.21 we review the 
CTP formalism and its application to the QTEs and discuss in detail the formulation of 
density-dependent Green's functions. In Section 12.31 we compute the CP-violating source 
terms, providing a check of Ref. [20], as well as the transport coefficients of the chiral charge 
densities. Here, we also enumerate the approximations used to obtain a set of coupled, linear 
differential diffusion equations, discuss their limits of validity, and identify additional terms 
(usually assumed to be sub-leading) that we defer to a future study. In Section ITU we solve 
these equations for the baryon density. A reader primarily interested in the phenomenologi- 
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cal implications may want to turn directly to Section [2.51 which gives illustrative numerical 
studies using the parameters of the MSSM. A discussion of the implications for EDMs also 
appears here. Section 12.61 contains a summary and outlook, while several technical details 
appear in the Appendices. 



2.2 Non-equilibrium Transport: CTP Formulation 

In what follows, we treat all CP-violating and non-topological chirality-changing interac- 
tions perturbatively^. In contrast to zero-temperature, equilibrium perturbation theory, 
however, the perturbative expansion under non-equilibrium, T > conditions requires the 
use of a more general set of Green's functions that take into account the non-adiabatic 
evolution of states as well as the presence of degeneracies in the thermal bath. Specifically, 
the matrix element of any operator 0{x) in the interaction representation is given by: 

(n|4,T{0(x)5int}|n), (2.3) 

where 

5int = Texpfi / d^xAnt) (2.4) 



for an interaction Lagrangian Ci^t, T is the time-ordering operator, and |n) is an in-state. 
Inserting a complete set of states inside Eq. (|2.3|) . we obtain: 

Y,{n\sL \m) HT{0{x)Si^t} \n) . (2.5) 

m 

In ordinary, zero-temperature equilibrium field theory, the assumptions of adiabaticity and 
of non-degeneracy of the states |n) imply that only the single state m = n contributes to 
this sum, so the only impact of S^^^ is the introduction of an overall phase, allowing one to 
rewrite Eq. 1)2. 5|) as: 

(n|T{0(x)5int}|n) 

/ ic I \ • 2.6) 

This simplification is no longer valid for non-equilibrium T > evolution, and one must 
take into account the action of 5^^^. appearing to the left of 0{x) in (|2.Hj) . Doing so is 



■^Sphaleron transitions, however, are manifestly non-perturbative, and we parameterize their effects in 
the standard way. 
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Figure 2.1: Closed time path integration contour. Fields (j) ^-^'s distinguished according to 
their placement on the forward (</>+) or backward ((/>-) portions of the contour. 



facilitated by giving every field in S'int and Sl^^ a "+" and "— " subscript respectively. The 
matrix element in 1)2. 3(1 then becomes: 



n) 



(2.7) 



where the path ordering operator V indicates that all "+" fields appear to the right of all 
"— " fields, with the former being ordered according to the usual time-ordering prescription 
and the latter being anti-time-ordered [here, 0{x) has been taken to be a field]. Note 
that the two integrals in the exponential in (|2.7j) can be written as a single integral along a 
closed time path running from —00 to +cxd and then back to —00: 



f dt [ (fx{C+ -£-)= [ dt [ d^xC, 
J -00 J Jc J 



(2i 



where the time t on the right-hand side is integrated over the contour C shown in Fig. 12.11 



Perturbation theory now proceeds from the matrix element ()2.7|) along the same lines 
as in ordinary field theory via the application of Wick's theorem, but with the more general 
V operator replacing the T operator. As a result, one now has a set of four two-point 
functions, corresponding to the different combinations of "+" and "— " fields that arise 
from contractions. It is convenient to write them as a matrix G{x,y): 

( GHx,y) -G<(x,y) \ 
G{x,y)= ^ : '^^ (2.9) 

\ G>{x,y) -G\x,y) ) 
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where 

G>(x,y) = (0_(x),/.V(y)) (2.10a) 

G<{x,y) = {(t>l{y)^+{x)) (2.10b) 

G\x,y) = {T{(l)+{x)<l,\{y)]) = e{xo - y^)G> {x,y) + e{y^ - xo)G<[x,y) (2.10c) 

G\x,y) = {f{(l)^{x)4>l{y)])=e{xo-yo)G<{x,y)+e{y^-XQ)G>{x,y), (2.10d) 

and where the ( ) denote ensemble averages, 

{0{x))^^Tt[pO{x)]. (2.11) 

Here p is the density matrix containing information about the state of the system, and 
Z = Tr (/5). In thermal equilibrium p is time-independent and is given byp = e-^(^-^»^«)for 
a grand-canonical ensemble. Note that the matrix G{x,y) may be written more compactly 
as: 

G{x,y)ab = {r[Mx)4iy)})i^3)bb- (2.12) 

The presence of the rs factor is a bookkeeping device to keep track of the relative minus 
sign between the and £_ terms in Eq. (|2.7|) . 

The path-ordered two-point functions satisfy the Schwinger-Dyson equations: 



G{x, y) = G°(x, y) + j d^w J d^z G°(x, w)^{w, z)G{z, y) (2.13a) 
G(x, y) = G°(x, y)+ f d^w [ d^z G{x, w)i:{w, z)G°{z, y) , (2.13b) 



where the "0" superscript indicates a non-interacting Green's function and where S(x,y) is 
the matrix of interacting self energies defined analogously to the G{x,y). An analogous set 
of expressions apply for fermion Green's functions, with an appropriate insertion of —1 to 
account for anticommutation relations. 



2.2.1 Quantum Transport Equations from CTP Formalism 

The Schwinger-Dyson Eqs. (|2.13|) are the starting point for obtaining the transport equa- 
tions governing ni(x). To do so, we follow Ref. [201 and apply the Klein-Gordon operator 



to G{x, y). Using 
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+ m2) y) = {Uy + n?) G°(x, y) = -i6^^^ (x - y) 



(2.14) 



gives 



(□^ + m^) y) = {x-y)-i J (fz S(x, z)G(z, y) (2.15a) 

(□y + m^) G(x, y) = -if^^'^) (x - y) - i f d*z G{x, z)S(z, y) . (2.15b) 



It is useful now to consider the (a, 6) = (1,2) components of these equations: 

{a, + rr?)G<{x,y) = -i J d^z \^\x,z)G<{z,y) -T.<{x,z)G\z,y) 
{ny + m^)G<{x,y) = -i f d^z \G\x,z)T,<{z,y) - G<{x,z)T^{z,y) 



Subtracting Eq. ()2.16b|) from Eq. H2.16a() and multiplying through by i gives 



i{n,-Uy)G<{x,y)l^y^^ 



id^{d^,-d^) G<{x,y)l__ 



--y=X ■ 



However, 



{d^,-d^,)G<{x,y) 



x=y=X 



(2.16a) 
(2.16b) 



(2.17) 



(2.18) 



where j^(x) = i{: 4^ {x)d^(t){x) :) = (^^(/,(x), j^(x)), since the "+" and "— " labels simply 
indicate the order in which the fields (/>^(y) and 0(x) occur and may be dropped at this 
point. Finally, expressing G*'*(x,y) and S*'*(x,y) in terms of 0-functions as in Eqs. 1)2.10(1 . 
we obtain from Eq. (|2.17p : 



Xo 



dzf) 



^>{X,z)G<{z,X) - G>{X,z)^<{z,X) 
+G<{X,z)J:>{z,X) -J:<{X,z)G>{z,X) 



(2.19) 



Following similar steps, but taking the sum rather than the difference of the components of 
the Schwinger-Dyson equations involving the 5^(x,y) component on the LHS, one obtains 
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the analogous continuity equation for Dirac fermions: 
dXn 



d^z / dzoTi ^>{X,z)S<{z,X) - S>{X,z)^<{z,X) 

J —oo 

+S<(X,z)S>(z,X) - S<(X,z)S>(z,X) 



(2.20) 



where 



and 



4{x) = {-4^^i:{x):), (2.21) 



S>^{x,y) = {i;.a{x)i^My)) (2-22a) 
Sa^i^, y) = -(^-/3(y)^+a(x)) , (2.22b) 

displaying explicitly the spinor indices a, (5. Note that the overall sign of the RHS of 
Eqs. 1)2.191 12.2U() differs from that in Ref. [20] since the definition of our Green's functions 
G{x, y) and S{x^ y) differ by an overall factor of —i. 

In many extensions of the SM, one encounters both chiral and Majorana fermions, 
which carry no conserved charge. It is useful, therefore, to derive the analogous continuity 
equation for the axial current j^^{x) = {ip{x)'yfj_j5ip{x)) . Doing so involves multiplying the 
Schwinger-Dyson equations by 75, performing the trace, and taking the difference rather 
than the sum of the components involving S''{x,y) on the LHS. The result is: 

1^ + V-j5(^) =2imP{X) (2.23) 



J d^z j "dzoTrj ^>{X,z)S<{z,X) + S>{X,z)T.<{z,X) 

- S<{X, z)S>(z, X) - S<(X, z)S>iz, X)] 75} , 



where P{x) = {'iIj{x)^^'iIj{x)) and m is the fermion mass. In principle, one could evaluate 
P{x) using path-ordered perturbation theory as outlined above. 



2.2.2 Power Counting of Physical Scales 

Evaluating the various terms in Eqs. (|2.191 I2.2U|) leads to a system of coupled quantum 
transport equations for the charges that ultimately determine nL{x). On the LHS of these 
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equations, it is conventional to parameterize j = —D[^n), in terms of the diffusion coef- 
ficient D (whose expressions we take from Ref. |5U|'). The RHS involves integrating the 
products of various Green's functions and self-energies over the entire history of the system. 
In practice, this integral depends on the various time and energy scales associated with 
non-equilibrium dynamics at finite temperature and density. Here, we observe that there 
exists a hierarchy among these scales that leads to a natural power counting in their ratios 
(generically denoted here as e) and that provides for a systematic expansion of the RHS of 
the transport equations ^T^^TM^T^. 

The changing geometry associated with the expanding region of broken symmetry and 
the spacetime variation of the Higgs vevs leads to a decoherence of states that have, initially, 
precise energy and momentum. The effect is analogous to the quantum mechanical evolution 
of a particle in a box of side L. If the value of L is changed to L + AL in some time 
interval At, a state that is initially a stationary state for the original box will become an 
admixture of the stationary states of new box. The shorter the interval At or the greater the 
wavenumber k of the initial state, the smaller the probability will be of finding the particle 
in the state with the same wavenumber in the new system. The time scale that characterizes 
this decoherence, r^, is naturally given by ~ where v = AL/At is the velocity of 

expansion of the box and k = p/h. In the present case, the relevant velocity is just Vw, the 
expanding bubble wall velocity, the relevant effective wave number k depends on |k| and 
the wall thickness, Ly^. The smaller the velocity or the longer the wavelength, the more 
adiabatic the dynamics of the expanding bubble become and the longer the decoherence 
time. Equilibrium dynamics are approached in the adiabatic limit: — > oo. The need to 
employ the CTP formalism follows from being in a situation with i;^^ > 0, or < oo. 

A second time scale that one encounters in quantum transport at the phase boundary 
arises from the presence of degeneracies among states in the thermal bath that vanish in 
the T ^ limit. At finite T, for example, a single, on-shell fermion may be degenerate 
with another state involving an on-shell fermion-gluon pair — a situation that is forbidden 
at T = 0. Interactions of strength g that cause mixing between such degenerate states 
give rise to thermal — or plasma — widths Tp of order aT with a = /A.-k, and transitions 
between the degenerate states take place on a plasma time scale Tp of order ~ l/Tp. Again, 
the use of the CTP formalism is necessitated when Tp < oo or T > 0. 

A third time scale, which we denote Tint, is associated with the intrinsic frequency ujk of 
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the quasiparticle states that characterize the plasma dynamics. This time scale is naturally 
given by Tint ~ l/t^fc- In the present case, we note that although the decoherence and 
plasma times are finite, they are typically much larger than Tint- For example, T\^t/T(i = 
Vwk/ujk < Vw/c. Numerical studies indicate that v^/c <^ 1. Similarly, Ti^t/Tp = oiT/uJk. 
Since quasiparticle thermal masses are of order gT or larger, one also has that the latter 
ratio is smaller than unity. Thus, one is naturally led to expand the RHS of the transport 
equations in these ratios: 



Finally, we observe that the generation of baryon number takes place in an environment 
of finite, but small particle number (or chiral charge) densities rij that are associated with 
chemical potentials /Xj. For the temperatures and densities of interest here, one has \l^i\/T <ti 
1, so that the latter ratio also provides for a natural expansion parameter. Denoting each of 
the ratios'^ in Eq. ()2.24|) and fJ-i/T by e, we show below that both the CP-violating sources 
and the relaxation term first arise at ©(e^), and we truncate our analysis at this order. We 
note that doing so introduces some simplifications into the evaluation of the RHS of the 
transport equations. For example, both the self energies and the Green's functions G^, 
etc. depend on thermal distribution functions f{T, fn) that differ, in general, from their 
equilibrium values, fo(T,fii). The difference 6f = f(T,fii) — fQ(T,fii) that characterizes 
the departure from equilibrium will be at least of 0(e), since it must vanish in the ^ 
limit. We find below that the effect of having 6f =/ contributes at higher order in e than 
we consider here, so that we may use the equilibrium distribution functions in the Green's 
functions and self-energies. 

2.2.3 Green's Functions at Nonzero Temperature and Density 

The computation of the various components of G{x, y) and y) appearing in Eqs. (|2.19| 
I2.2U|) at nonzero temperature and density requires knowledge of (T, /ij)-dependent fermion 
and boson propagators. The T-dependence of propagators has been studied extensively 

^For our purposes, it is not necessary to distinguish a hierarchy among the different scale ratios, as we 
work to leading nontrivial order in e. 



< Tint/Td < 1 



(2.24a) 



< Tint/Tp < 1 . 



(2.24b) 
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(see, for instance, Ref. |SJ and references therein), while the /ij-dependence of fermion 
propagators has been studied in Refs. Here we summarize the features of (T, ^Uj)- 

dependent propagators which are important for our subsequent application of the real-time, 
CTP formahsm of Sec. 12.21 and give some more technical details in Appendix [2^ 



For pedagogical purposes, we provide here a brief derivation of the non-interacting 
fermion propagator but only give final results for the case of interacting fermions and 
bosons. To do so, we start from the mode expansions for the field operators appearing 
in the free Dirac Lagrangian, ^(x) and tp^x): 



^(x) = / E [a^u'^i^)e-''-'' + bZ^v"ik)e^'-^] (2.25a) 

^(x) = / (03^ E [<^^"(k)e^'=-^ + 6£t;"(k)e-^-] , (2.25b) 



where = (cJkjk), U]^ = y |k| +m'^, the mode operators satisfy: 

{a^,a^J} = {b^yj} = (27r)35(3)(k-k02^k<5"^ (2.26) 

and 

= /(^fe,;u.)(2vr)3<5(3)(k - k')2uj^5^f' (2.27a) 

K^€') = fi^k, -^i^){27Tf5^^Hk - k')2u;k<5"^ (2.27b) 

with f(ij,fii) being the non-equilibrium Fermi distribution function. For our purposes, 
the relative change 5 f {lo , jii) / fo{uj , fii) enters the transport equations multiplying explicit 
factors of Tp and either v^j or fi, so that in working to second order in e we may replace 
/ by the equilibrium distributions /o(w,/ii) = npiuj — fii) = [e^^~^^^^^'^ + 1]~^. Using 
the mode expansion (|2.25|) it is straightforward to show that S^{x,y) = {'ip(x)'ip{y)) and 
5^(x,y) = —{'4}{y)il){x)) can be expressed as: 

5^(x, y) = j ^e-"'<^-y^g${ko,fii)p{ko, k) (^ + m) (2.28) 
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in terms of the free particle spectral density: 



p(A:o,k) 



2wfc 



1 



1 



/co — Wfc + ie kQ + ijjk + ie J \kQ — uj^ — ie ko + ujk — ie 



1 



. (2.29) 



and the functions: 



gpiko,Hi) = 1 - npiko - 



(2.30a) 
(2.30b) 



The propagators 5*'*(x,?/) can now be constructed from the S^{x,y) as in Eqs. ()2.10|) . 



In the presence of interactions (characterized by a generic coupling g), the fermion 
propagator becomes considerably more complicated than given by Eq. (|2.28j) . In particular, 
single fermion states can mix with other multiparticle states in the thermal bath, leading 
to the presence of additional poles (the "hole" modes) in the fermion propagator [531 154j . 
The general structure of the fermion propagator arising from these effects has been studied 
extensively at zero density jHS]- In Appendix A we generalize to the case of non-zero /ij. 
For massless fermions, the resulting propagators are given by: 



S<{x,y]Hi) 



d^k 



-ik-(x—y) ^ 



p+(/co,k,//i) H P-(A;o,k, m 



(2.31) 



where k is the unit vector in the k direction, and 



P+(^o,k,/Xi) 



Zp{k, Hi 



ko-£p{k,Hi) kf) - £p{k,Hi)* 



+ 



Zh{k, -HiY 



Zh{k, -Hi 



ko + £h{k,-Hi)* ko + £h{k,-Hi) 



+ F{kQ,k,Hi)* - F{ko,k,iJ,i 



(2.32) 



and 



P-{ko,k,Hi) = [p+{-kQ,'k,-Hi 



(2.33) 



Here, £p{k,fj,i) and —£h{k, —Hi)* are the two (complex) roots (in /cq) of the equation: 



= feo - k + D+{ko,k, Hi) + ie 



(2.34) 
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where iD±{kQ,k, are contributions to the inverse, retarded propagator proportional to 
(70 -F7"k)/2 arising from interactions. The function F{kQ,k, gives the non-pole part 
of the propagator, and k = |k|. We find that the resonant contributions to the particle 
number-changing sources arise from the pole parts of the propagators, so from here on we 
neglect the terms containing F{kQ, k, fii). 

In the limit 5 — > 0, one has Zh ^ and Zp ^ 1, recovering the form of the propagator 
given in Eq. H2.28() . For nonzero g, however, Z^ is not of order since the particle and 
hole modes arise from mixtures of degenerate states. In particular, at k = one has 
Zp = Zfi = 1/2. As k becomes large (of order the thermal mass or larger), Z^/Zp <^ 1, and 
the particle dispersion relation is well-approximated by £p = |k|'^-|-m^(r, fii), where m(T, fii) 
is the thermal mass. In our particular application to the MSSM, the gaugino Mi masses 
will typically be taken to be of order several hundred GeV, and for the SU(2)/,xU(l)y 
sector, thermal effects do not induce substantial mass corrections. We find that the gaugino 
contributions to the RHS of Eqs. (|2.191 [^^0]) are dominated by momenta of order Mj, so 
that the hole contributions to the gaugino S^{x, y) can be neglected. In contrast, for quarks 
we find non-negligible contributions from the low-momentum region, so we retain the full 
structure given by Eqs. (|2. 31112. 37?]) in computing their contributions. 

It has been noted in previous studies of quark damping rates that the one-loop thermal 
widths Tp^h = Iiii<fp,/i(A;, ^) are gauge-dependent (see Ref. (SHI and Ref. [3] therein), whereas 
the thermal masses m{T, fi) entering Sp^h are gauge-independent to this order. Gauge- 
independent widths can be obtained by performing an appropriate resummation of hard 
thermal loops (HTLs) [511 1561 IF7]. The latter are associated with momenta /co,A; ~ gT, 
for which the one-loop functions D±{kQ, k, fi) are of the same order in g as the tree- level 
inverse propagators. In what follows, we will estimate the widths Fp based on existing 
computations of damping |SSJ \^ 1^0], deferring a complete computation of the gauge- 
invariant, /ij-dependent contributions in the MSSM to a future study. In general, the 
residues Zp ^ also carry a gauge-dependence, and at this time we are not aware of any 
HTL resummation that could eliminate this dependence. In principle, elimination of this 
gauge-dependence requires inclusion of one-loop vertex corrections in the computation of 
the S^(x,y) and S^{x,y) appearing on the RHS of Eqs. ()2.19l ir^U]) . and we again defer a 
complete one-loop computation to a future study. 

The derivation of the finite-density scalar propagators proceeds along similar lines. 
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Starting from the mode expansion of the free scalar field (p{x) in terms of plane-wave solu- 
tions to the Klein-Gordon equation and following analogous arguments as for fermions, one 
arrives at the following scalar Green's functions: 

G^i^^y) = I ■0^e-''^<--y^g^{ko,fi^)p{ko,-k), (2.35) 

where the equilibrium distribution functions are: 

g'^{uj,fi) = 1 + nB{uj - fii) (2.36a) 
g^{uj,fj.) = nsiuJ - tJ-i) , (2.36b) 

with ub^x) = l/(e^/"^ — 1) and p(A;o,k) given by Eq. 1)2. 29() . As with fermions, one may 
include the effect of thermal masses and widths by replacing m? — > m?(T,fii) and ie 
ie + iT{T, m). 

2.3 Source Terms for Quantum Transport 

The expressions for G^{x, y) and S^{x, y) now allow us to compute the perturbative contri- 
butions to the source terms on the RHS of Eqs. (|2.19I2.2U() starting from a given electroweak 
model Lagrangian. Here, we work within the MSSM as an illustrative case, but emphasize 
that the methods are general. The Feynman rules giving the relevant interaction vertices in 
the MSSM are taken from Ref. ^7j) what follows, we only write down those relevant 

for the computations undertaken here. It is useful, however, to place our calculation in a 
broader context by considering the various classes of graphs that generate different terms 
in the QTEs. The simplest topologies are those involving scattering of particles and their 
superpartners from the spacetime varying Higgs vevs (generically denoted v) in the plasma 
[Fig. 12. 2j . These graphs give rise to both the CP- violating source terms discussed in Ref. 
j2Uj as well as terms proportional to chiral charge. The latter involve the number densities 
of at most two different species, such as the left- and right-handed top quarks [Fig. I2.2r a')] 
or their superpartners [Fig. I2.2f b')]. For purposes of illustration, we follow Ref. [201 and 
work in a basis of mass eigenstates in the unbroken phase, treating the interactions with 
the Higgs vevs perturbatively. This approximation should be reasonable near the phase 
transition boundary, where both the vevs and their rate of change are small, but it clearly 
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^'(^) V{X) y{y) 




Figure 2.2: Contributions to the relevant self-energies from scattering of particles from the 
spacetime varying Higgs vevs. 

breaks down farther inside the bubble wall, where the vevs become large (of order the phase 
transition temperature, Tc). In general, one would like to perform a resummation to all or- 
ders in the vevs, possibly employing the approximation scheme proposed in Refs. [481 161j . 
We postpone a treatment of this resummation to a future study^. 

Yukawa interactions involving quarks (squarks) and Higgs (Higgsinos) are illustrated in 
Fig. 12.31 (the self-energies are obtained by amputating the external legs). 

These interactions cause transitions such as / fH, f ^ fH, and / ^ fH. Con- 
tributions from gauge interactions appear in Fig. 12.41 The latter induce transitions of the 
type / <-> fV, f <-> fV, and / <-> fV. In general, one expects the Yukawa and gauge 
interactions involving three different species to depend on sums and differences of the cor- 
responding chemical potentials, as in fif — fi^ — Hy for the supergauge interactions. In 
previous studies, it has been assumed that the gauginos V are sufficiently light and the 
coefficients of the corresponding terms in the QTEs sufficiently large than one has ^Uy ~ 
and fif ^ nj. Although the quantitative validity of this assumption could be explored using 
our framework here, we defer that analysis to a future study and take fiy ^ 0, fif ^ fj,^. 
Consequently, one may, as in Ref. [T^, define a common chemical potential for SM particles 
(including the two Higgs doublets) and their superpartners. 

^The authors of Ref. |48| find that carrying out such a resummation reduces the resonant enhancements 
of the CP- violating sources, but they did not consider the CP-conserving, chirahty-changing terms that are 
our focus here. The consistency of the proposed approximate resummation with our power counting remains 
to be analyzed. 
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Figure 2.3: Contributions to the relevant self-energies from Yukawa interactions. 

In previous studies, it has also been assumed — based largely on simple estimates (see, 
for instance, Ref. JH]) — that the Yukawa interactions of Fig. 12.31 are sufficiently fast that 
they decouple from the set of QTEs, leading to relations between the chemical potentials for 
the Higgs (Higgsino) fields and those for matter fields. For example, Yukawa interactions 
that couple the Higgs doublet fields H with those of the third generation S\J{2)l doublet 
quarks, Q with the singlet top quark supermultiplet field, T, generate terms of the form: 

Ty (^Q - + m) ■ (2-37) 

To the extent that Ty is much larger than the other transport coefficients appearing in 
Eqs. ()2.19I2.2U() . one has hq = fix — /^H plus terms of ©(l/Py). The remaining terms in 
the QTEs will involve the CP-violating sources, sphaleron terms, and terms that couple 
left- and right-handed chiral charges, such as TMitJ-Q — Mt)- Again, this assumption could 
be tested using the current framework, but the computation of Ty is considerably more 
arduous than those discussed below, where we focus on the CP-violating sources and the 
FM-type terms that are generated by the diagrams in Fig. 12.21 

2.3.1 Bosons 

We consider first the scalar interactions in Fig. I2.2r a) . The largest contributions involve the 
L and R top squarks, tL,R owing to their large Yukawa coupling, yt. In the basis of weak 
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Figure 2.4: Representative contributions to self-energies from (super)gauge interactions. 

eigenstates, the relevant interaction Lagrangian is: 

£ = ythi*R{AtVu - fJ,*Vd) + h.c. , (2.38) 

where Vu,d are the vevs of H'^ ^, and we take v = ^Jv'^ + v"^ and tan /3 = Vu/vd- Note that in 
Eq. (|2.38j) we allow the Vu,d to be spacetime-dependent. In the region of broken electroweak 
symmetry and stable vevs, we have nit = UtVu- 

Using the Feynman rules for path-ordered perturbation theory, it is straightforward to 
show that the diagrams in Fig. I2.2r al generate contributions to T,ji(x,y) of the form: 

I^R{x,y) = -g{x,y)Gl{x,y), (2.39) 

where 

gix, y) = [Atvu{x) - iJ*Vd{x)\ [A*tVu{y) - /iWd(y)] . (2.40) 
Substituting Eq. (fT^ into Eq. leads to: 

d,i^^{x) = Si^{x) (2.41) 

for right-handed top squarks, where is the corresponding current density and the source 
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Si^{x) = - J d'z P dzo !^[g{x,z) + g{z,x)]Re[G>{x,z)G<iz,x) - G<{x,z)G>{z,x)] 

+i[g{x,z) - g{z,x)]lm[G^{x,z)G'^{z,x) - Gf{x, z)G%{z,x)]'^ , 

(2.42) 

where the L,R subscripts indicate the propagators for the L and R top squarks. 

The first term in the integrand of Sf^{x) is CP-conserving and leads to the rA/-type 
terms discussed above, while the second term in the integrand provides the CP-violating 
sources. We concentrate first on the former. Expanding g{x, z) about z = a; it is straightfor- 
ward to show that only terms involving even powers of derivatives survive in g{x, z)+g{z, x). 
Under the assumptions of gentle spacetime dependence of the Vi{x) near the phase boundary, 
we will neglect terms beyond leading order and take g{x,z) + g{z,x) 2g{x,x). Conse- 
quently, the CP-conserving source is: 



Sf^ix) 



Sz j dzQ [Gl{x,z)G^{z,x) -G-l{x,z)G'^{z,x)\ (2.43) 
X [5b(^o,^l)5b(9o,^k) - 5B(^o,AtL)c/B(9o,/^/?)] , 



with 



gix,x) =y^[\^l\^vj{x) + \At\\l{x)-2vdix)vuix)Re{^xAt)]. (2.44) 
Note the simplification 

5b(^o,/^l)5b(9o,^r) - 9B{ko, iiL)gB{(lo^ m) = nsiqa - m) -nsiko - pi)- (2.45) 

Performing the d^z integral leads to a (5 function in momentum space. After carrying out 
the d^q integral, we perform the k^,q^ integrals by contour integration^, expand to first 



^Here and in subsequent equations, we show only the terms arising from picking up the residues of the 
poles in the spectral functions such as pL(fc'', k), p_R(q'', q). A careful calculation would also include the 
residues from the poles in the thermal distribution functions such as nsik'^), hsik'^), etc. We relegate these 
to Appendix. 12. CI where we find their contribution to the final numerical results to be unimportant compared 
to the terms retained in the main text. 



order in fJ-L,R/T, and obtain: 
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\AtVu{x) - fl*Vd{x)\' 



X lni< — — h 



(2.46) 



where 



u;ln = \k\^ + Ml^^ (2.47a) 
£l,r = ^L,R - i'^L,R (2.47b) 

and A/j^ ^,^L,R are the thermal masses and widths for the Il^r, and the factor of Nc comes 
from summing over the colors. Note that, in arriving at Eq. (|2.46j) . we have neglected the /Uj- 
dependence of the pole residues Z{T, fiL^R), thermal frequencies, ujl,r{T, fiL^R), and widths, 
^l,r{T, [il^r). The effect on S'?^{x) of the /ij-dependence of the residues and thermal 
frequencies is sub-leading in the gauge and Yukawa couplings, whereas the effect from the 
thermal widths occurs at leading order. The /Xj-dependence of T l^r{T, ^l^r) is simply not 
known, however, so we do not include it here. A more explicit expression for the dependence 
of Sf^{x) on the thermal frequencies and widths is given in Eqs. (|2.128B2.13()j) of Appendix 



For purposes of future analysis, it is useful to rewrite Eq. ()2.46|) as: 



5?^ = r+(/ii + ^jlr) + rr(;xi - , (2.48) 



where 

* T 47r^ Jq ujRLOL [ ^L-c]^ CL+ £r 

(2.49) 

Before proceeding with the CP-violating source, we comment briefly on the structure 
of Eqs. (|2. 48112. 49|) . In particular, we note that 

(i) Terms of the type are absent from the conventional QTEs for EWB. It is straight- 
forward to see that in the absence of interactions that distinguish between tL and tR, 
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= 0, as the integrand of Eq. 1)2.49(1 is antisymmetric under L ^ R interchange. 
In contrast, the transport coefficient is nonzero in the hmit of exact ij^ <-> iji 
symmetry. This term corresponds to the usual damping term in the QTEs associated 
with scattering from the Higgs vevs. 

(ii) In the absence of thermal widths r^^/j, the quantity in brackets in Eq. H2.49() is purely 
real, and so the damping term would be zero. 

(iii) The structure of the energy denominators implies a resonant enhancement of the 
integrand for M? ~ M? . A similar effect was observed to occur for the CP- violating 
sources (see below) in Refs. [201 162| . The expression in Eq. (|2.49|) makes it clear that 
the relaxation terms display a resonant behavior as well. The resulting quantitative 
impact of this resonance on the baryon asymmetry is discussed in Sect. 12.51 

Properties (ii) and (iii) are shared by all source and damping terms, we discuss below. Note 
that the explicit factors of ^l,r/T and property (ii) imply that, away from the resonance 
region, Sf^ is 0{e^). 

The computation of the CP- violating source, given by the second term in Eq. (|2.42() . 
proceeds along similar lines. In this case, the coefficient [g{x, z)—g{z, x)] vanishes for x = z, 
so we must retain terms at least to first order in the expansion about x = z: 

g{x, z) - g{z, x) = 2iyj Im(/xAi) [vd{x)vu{z) - Vd{z)vu{x)] 

(2.50) 

= 2iy1 lm{^At){z - x)^ [vd{x)d\Vu{x) - Vu{x)d\Vd{x)] H , 

where the indicate higher order terms in the derivative expansion that we neglect for 

the same reasons as discussed previously. When the linear term in Eq. (|2.5UI) is substi- 
tuted in Eq. ()2.42|) . only the time component yields a nonzero contribution. The spatial 
components vanish due to the spatial isotropy of the spectral density: (^^(/co) Ai)p(^O) k) = 
g^{kQ, fj,)p{ko, |k|). We may then make the replacement: 

g{x, z) - g{z, x) 2iyf lm{fj,At) [vd{x)vu{x) - Vd{x)vu{x)] {z - xf 

(2.51) 

= 2iyj lm.{iiAt)v{xf ^{x) {z - xf . 

In general, we expect (5 to be of order fu,/c, so that the CP-violating source is first-order 
in one of the small expansion parameters discussed earlier. Consequently, when evaluating 
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this term, we may neglect the ;UL^i?-dependence of the g^{ko,fi). After carrying out the 
(,ko,qo) contour integrals and performing the time integration, we obtain: 

(2.52) 

Again, property (ii), in conjunction with the factor of /3 oc Vw, implies that 5"?^ is 0{e^). 
An expression giving a more explicit dependence on the widths and frequencies appears Eq. 
H2.13ip of Appendix 12. Bl which we note agrees with that of Ref. j2Uj except for a different 
relative sign in front of the cos2i;^> term of that equation and the overall factor of Nq- 



2.3.2 Massive Fermions 

The computations for fermions proceed along similar lines. We consider first the source 
terms for Higgsinos. We recall that it is useful to redefine the Higgsino fields to remove the 
complex phase from the Higgsino mass term: 

^T" = ^ (V'HoV'/fO - '^H-^Ht) + {i^H^^Hl - ^H-jHt) (2-53) 

via 

V,- ^ ^ e-'^^Hl^^ (2.54) 

leading to: 

= (h^.HI - H,Hi) + 1^1 (5°t^ot _ ^;t^+t) . (2.55) 
Defining the four component spinors, 




(2.56) 



for the Higgsinos, and 




(2.57) 
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for the gauginos, leads to the Higgsino-gaugino-vev interaction: 



1 



(2.58) 



where 



1 T 7^ 

Pl,r = (2.59) 



while the mass terms 1)2. 55p in the Higgsino Lagrangian become: 

CT' = -m^H,^H, + ^H,^H.)- (2.60) 

Note that the spinors ^'^o ^ h+ satisfy a Dirac equation with Dirac mass |//|, even 
though the H^,^ are Majorana particles. The ^^y± are Dirac particles of mass M2, whereas 
the ^^yo and ^ are Majorana particles with Majorana masses M2 and Mi, respectively. 
We also note that the construction of the Dirac spinor ^ ^0 allows one to define a vector 
charge and corresponding chemical potential, /x^O) for the neutral Higgsinos, even though 
they are Majorana particles. In contrast, there exists no such vector charge for the ^^^o and 
ipj^o- One may, however, study the quantum transport of the axial charge of the Majorana 
fermions using Eq. 1)2. 23(1 . An attempt to do so for the neutral Higgsinos was made in Ref. 
j48j . though only the CP- violating sources were evaluated using non-equilibrium methods. 
The impact of the corresponding axial charge density on the baryon asymmetry was found 
to be small. We will return to this issue in a future study, and consider only the vector 
densities below. 



The most straightforward computation is that of the source terms. For notational 
convenience, we rewrite the chargino interactions in Eq. ()2.58|1 as: 

-92^ H+ [9l{x)Pl + 9r{x)Pr] + h.c. (2.61) 
In this case, the self-energy generated by Fig. I2.2r al is: 

S^±(x,y) = -gl [gL{x)PL + 9R{x)PR] %±(x,y) [gdyT Pr + gRivT Pl] ■ (2.62) 
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Defining: 



9l 



gA{x,y) = — [gL{x)gL{yT + gR{x)gR{y)*] 



9l 



9Bix,y) = — [gL{x)gR{yy + gR{x)gL{y)*\ 



(2.63a) 
(2.63b) 



we obtain for the RHS of Eq. ^HTM : 



Sjr± {X) 



J J-oo ■_, n I 



j=A,B 

[gj{x,z) +gj{z,x)] ReTr S^^{x, z)S'^^{z,x) 



+i[gj{x,z) - gj{z,x)]lm.T:i {x,z)S'i^{z,x) 



S^^{x,z)S^^{z,x) 



S^^{x,z)S^^{z,x) 



3 



J. 

(2.64) 



where the subscripts "A" and "B" on the traces denote the contributions arising from the 
^ and m terms, respectively, in the spectral function in Eq. 1)2.28(1 (an overall factor of 1/2 
due to the presence of the chiral projectors Pl,r has been absorbed in the definition of the 
9a,b)- 



As in the case of the scalar fields, the leading density-dependent, CP-conserving contri- 
bution to {x) arises from the term in Eq. (|2.64|) containing the x ^ z symmetric factors 
[gj (x, z) + gj (z, x)] . To lowest order in Vw , we may set rc = z in these factors. Using the spec- 
tral representation of the S^{x,y) given in Eq. ()2.28|) . including gauge-invariant thermal 
masses and widths, and expanding to first order in /ij/T, we obtain the chirality-changing 
source term: 

(^) = {^^w^ + ^^H±) + r ~± (/^vk± - , (2-65) 

where 



r± = -^v(x? 



dk 



Im 



Sy^Sjj - k + M2 I /i| cos (/)^ sin 2/3 



w 



+ + A;2 - M2 cos <j)^ sin 2/?] 



f^FiS^) =F hpiSfj) 



S I- £~ 



(2.66) 
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where the definitions of ^ and ^ are analogous to those given in Eqs. ()2.47() and 



^Hx) = ^^^JT^^ . (2.67) 

Also, the factor of cos (j)^ is very nearly 1 for the region of small </>^ in which we find ourselves 
in subsequent sections. The explicit dependence of F ~ on thermal frequencies and widths 



is given in Eq. (|2.1.32|) of Appendix I2.B 



In the present case, we follow Ref. |19) and assume no net density of gauginos, thereby 
setting ^J'^r± = in Eq. 1)2. 65(1 and giving: 

S%l{x) = -r^±/i^± , (2.68) 

with = + ^^±- ^'^^^ case, it is straightforward to obtain the corresponding 

source term for the neutral Higgsinos, 

S'g^i^) = -r^^o/^^o , (2.69) 

where can be obtained from the formulae for r^_|_ by making the following replacements: 
92 — ^ 52 / for intermediate states and §2 ^ Qi/ \/2, — > ujj^ , and F^ for the 

B intermediate states. 

The Higgsino CP- violating source arises from the second term in Eq. 1)2. 64() . As before, 
we expand the gj{x, z) to first order about x = z and observe that only the xq — zq compo- 
nent survives when the d^z integration is performed. Also note that gA{x,z) — gA{z,x) = 
2i ImgAix, z) = so that only the terms proportional to the Higgsino and gaugino masses 
contribute. The result is: 

2 

S^i^) = ^vixfp{x)M2 |/i| sm<P,, 

dkk^ (nF{£^)-nF{£t) l-nF{£^)-nF{£^)\ (^.70) 







cyp 

The corresponding expression for (x) can be obtained by making the same replacements 
as indicated above for the CP-conserving terms. The correspondence with the results of 
Ref. ;2flj can be seen from Eq. (|2.i:-{4|) of Appendixlaill We again find essential agreement, 
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apart from a sign difference on the cos 2(j) term. 



2.3.3 Chiral Fermions 

The final source term associated with Fig. \2.2\ a) involves L and R top quarks. At this 
order, the latter only contribute a /ij-dependent CP-conserving term. In order to illustrate 
the structure of this term that arises when the terms of 0{g'^) are retained, we employ the 
interacting fermion propagators of Eqs. 1)2. 3112. 33(1 . The result is: 

sf^{x) = r+ (/xi^ + J + r^^ (/ii^ - ^lt^) , (2.71) 



with 



■-p ) 



(2.72) 

Here, the "p" and "/i" subscripts indicate contributions from the particle and hole modes, 
and "L " and "i?" refer to left- and right-handed quarks. We have not included in our 
calculation the effects of fit^ ^-dependence of the widths ^p'^iT, r)^ which in principle 
also enter at this order. For an expanded version of Eq. ()2.72|) . including these effects, see 
Eq. frn^ in Appendix 12751 

In the limit oi tR symmetry, F^ vanishes, and F^ simplifies to: 



1 Ncy^vuixf r,.,2, \ Zp{k)\ , Zt,{k)\ 

^tR = 7f 2 / hpiSp) + — hF{£h) 

+ ^-^^^^[hF{£p) + h,[£l)] 

We observe that all contributions to the CP-violating source terms and the F^ vanish in 
the limit of zero thermal widths. Since the widths are generically of order g^T (here, g 
denotes either a gauge or Yukawa coupling), the source terms for the QTEs are generally 
fourth order in the couplings. 
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2.4 Quantum Transport Equations and pB 



We now discuss diffusion equations for the particle species that significantly contribute to 
the density of left-handed doublet fermions nL{x) [cf. Eq. (|2.2|) ] that acts as the "seed" 
for baryogenesis. We subsequently relate pB to ni and solve explicitly the equations in the 
case of a simple geometry and profile for the bubble wall describing the phase boundary. 

2.4.1 Solving the Diffusion Equations 

Using the source terms computed in Section[T21 one can arrive at a coupled set of differential 
equations for the various particle number densities. These equations simplify considerably 
under the assumptions of approximate chemical equilibrium between SM particles and their 
superpartners (/i/ fij with fiy 0), as well as the between different members of left- 
handed fermion doublets {fJ-w^ ~ 0). In this case, one obtains transport equations for 
densities associated with different members of a supermultiplet. This approach is the one 
followed in Ref. jl9j . and for pedagogical purposes we summarize the development here. 
First, we define the appropriate supermultiplet densities: 

Q = nt^+ n^^ + + (2.74a) 

T^nt^ + nj^ (2.74b) 

B = nb^+ n-^^ (2.74c) 

H = njj+ + uho - njj- - njjo + ^/^+ ~ "/f- + ''^iio ~ ''^h^ > (2.74d) 

where the Higgsino densities arise from the vector charges = ^'^_|^7^\1'^_|^ and n^g = 
^^oT^^^o associated with the Dirac fields defined in Eq. 1)2. 56() . There are analogous 
definitions for the first- and second-generation (s)quarks. Although we do not consider 
them here, one may also define the corresponding axial charge densities. In the case of the 
Higgsinos, for example, it will involve the sum, rather than the difference, of the u- and 
d-type Higgsino densities^ 

The diffusion equation for a density rii has the structure: 

= 5^ + 5,^+5^^ , (2.75) 



®This density was considered in Ref. j48l, and its overall impact on the baryon asymmetry found to be 
small. 
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where J^" is the current associated with the density n^, Sf^ and sf^ axe the source terms 
computed above, and S^^^^ is the strong sphaleron transition term, arising from the QCD 
anomaly of the axial quark current: 



« = (2.76) 



where = Qi^^^5<lii summed over flavors i. Various derivations of the strong sphaleron 
term appear in the literature, so we do not reproduce them here. However, we note that 
the expressions in Refs. |19| IH.Sj have erroneously omitted a factor of I/Nq |H4j . 

The CP-conserving damping terms have been given in Eqs. 1)2. 48(1 . (|2.65|) . and 
1)2. 71() to linear order in the appropriate chemical potentials. Assuming local thermal equi- 
librium we relate the number densities to the chemical potentials via: 



(2^ 



where N{uj,fi) is the appropriate boson or fermion distribution function and gi counts the 
internal degrees of freedom (spin and color). Dropping terms of 0{fi^), one obtains: 

himjT)T^ 

rii = ^i^ , (2.78) 

D 

where the factors ki{mi/T) are exponentially small in the regime mi/T ^ 1, and reduce 
in the massless limit to /cj(0) = 1 for chiral fermions, /ci(0) = 2 for Dirac fermions, and 
ki{G) = 2 for complex scalars. In our analysis we keep the full dependence on mi/T: 



k,{m,/T) = h{0)^ dxx——^^/x^^^n^, (2.79) 

Jm/T (e^±l)^ 

where for fermions (bosons) c^(b) = 6 (3), and we choose the -|-(— ) sign in the denominator. 



Using Eq. (ITTHl) in Eqs. (|23H1 ESSl EUD) , and defining: 



r± = A (rf + Tf ) (2.80a) 
Th = ^{T^^+T^,) , (2.80b) 
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the resulting set of coupled transport equations is: 



\kT kff kg) \kQ kr k 



t 



(T H g \ _ /2g T ^ %Q + T)\ ^0> 



+ ry — - — - - 2r,, - _ + ' - s 



\kT kn kqj ykq kr ks 



t 



where Vss = e/c'fafT, with k' ~ 0(1). 

We comment briefly on the structure of these equations. In previous derivations of these 
transport equations, the terms on the right-hand sides containing the various reaction rates 
r^jFy, etc. were derived using semi-classical statisical mechanics. Consider a microscopic 
reaction which changes the number densities rii of particles i each by an amount Anj. Let 
fii be the chemical potentials for these particles. In one reaction, then, the free energy F 
changes by and amount AF = J2- HiArii. In a thermal ensemble, the probabilities to occupy 
a state with a number of particles ("old") or with rii + Aui ("new") are proportional to: 



(2.82) 

^ new ^ ) 

where /? = l/T. The net rates to transition from the old to the new state or vice versa are 
proportional to these occupation probabilities: 

^q.miPaA - Pnew) = rg.m.Poldil " 6^ '^''^"0, (2-83) 

where m is the quantum mechanical rate for the individual reaction rij ^ rij + Arij. 
Summing over all possible starting states "old", we have X^g^ates -^oid = 1, so the net rate, 
in the limit of small chemical potentials /ij/T ^ 1, is: 

rnet = -r,.„.i:^. (2.84) 
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Then, the rate of change of a particular density rij, since one reaction changes it by Anj, 
is given by: 

n, = -r,.„An, j;^. (2.85) 

i 

The minus sign indicates that the reactions wih tend to cause rij to relax back to zero. The 
chemical potentials /Uj are related back to the densities rii themselves by Eq. (|2.78jl . 



This sort of argument accounts for the structure of the T~j^j,TY,^hj and Tg s. terms 
in (|2.8H) . In the treatment of Ref. ^HI) these are the only terms that appear, with the 
thermodynamic rates being given simply by: 

r = (2.86) 

However, deriving the equations by starting from the Schwinger-Dyson Eqs. ()2.13|) in 
nonequilibrium quantum field theory, we obtain more general combinations of densities 
on the right-hand sides of Eqs. 1)2. 81() . such as the T~j^ terms, with more involved expres- 
sions for the rates T than Eq. (|2.86|) . In this work, we have only completed this computation 
for the T^j,Th terms, leaving the recalculation of Fy to a future study. For Ts,s. we adopt 
the semi-classical derivation. 



For now, though we have not computed Fy, we will follow the authors of Ref. who 
estimate Fy F^^. For k' ~ 0(1), one also has F^^ ^ F^j. These facts allow one to relate 
algebraically the densities Q and T to H, by setting the linear combinations multiplying 
Fy and Tgs equal to 5y = 0(1/Fy) and 6ss = 0{1/Tss), respectively. One then obtains: 

{kB - 9kT)kQ A , X \ 

Q = 7777 — -TT — , , H + agydY + aqsdss (2.87a 

(9fcT + 2kB)kT ^ . ^ r S7K^ 

{9kT + 9kQ + kB)kH 

with known coefficients aQY^Qs,TY,Ts- Taking 2 x [Eq. ()2.81a() ] + [Eq. I|2.81b|) ] + [Eq. 

(|2.81c|) ]. introducing the diffusion approximations T = —Dg'VT, Q = —Dg'VQ, H = 
-DhVH, and using Eq. ((TH7)) leads to: 



H - DV^H + rH-S = 0{6ss,Sy) 



(2.88) 
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where''' 



^ ^ {9kQkT + kskQ + AkTkB)Dq + %(9A;r + Q/cq + A:B)i^fe ^ 



- _ (9fcQ + 9fcT + kB){TM + Tfe) - (3fcB + 9fcQ - 9fcT)r+ 
'CT + kskq + 4A:TA;ij + /c//(9A:t + 9A:q + fc^) 



(2.89b) 

S = kH{^kQ+9kT + kB) / gQP^ gC/t\ .2 89c) 

()k^kT + kBkQ + 4:kTkB + kHi9kT + 9kQ + kB)\^ 7 



The subleading terms dY,ss can be determined by use of Eqs. (|2.87p in Eqs. H2.81al2.81b() . 
We include the effect of 5ss in our final expression for pB although its effect is negligible 
in the relevant MSSM parameter region. 

Equation 1)2. 88(1 can now be solved for a given set of assumptions about the geometry of 
the bubble wall. Again, for clarity of illustration, we will work in a framework that allows 
us to carry analytic calculations as far as possible, leaving to the future a numerical solution 
of the equations for a realistic wall geometry and profile. First, as commonly done in earlier 
studies, we ignore the wall curvature, thereby reducing the problem to a one-dimensional 
one in which all relevant functions depend on the variable z = \x + Vyjt\, where is the 
wall velocity. Thus, z < is associated with the unbroken phase, z > with the broken 
phase, and the boundary wall extends over < z < L^. Second, we take the relaxation 
term f to be nonzero and constant for z > 0. The resulting solution for H in the unbroken 
phase z < (related to as shown below) is: 

H{z)=Ae''^^/^ (2.90) 



with 

The above equation is valid for any shape of the source S{z). For simplicity, however, we 
assume a simple step-function type behavior for the source: S nonzero and constant for 
< z < L^. Specializing to this case of constant sources in < 2 < L^, using ADT ^ v'^, 



^Our expressions differ from those in Ref. )19| . which we beheve result from an algebraic error. The 
numerical impact of this difference, however, is not significant. 
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Lw \/r/D <^ 1, and taking F = rp (F/, + F^) from Eq. (|2.89ap . we arrive at: 



(2.92) 



M 



When evaluating the source terms S^,S?^ [see Eqs. (|2.52p . (l2.7UI) ] for this simple profile 
one has to use (5 = Vw/S.j3/Lw'- thus A is explicitly proportional to Vw and is only weakly 
dependent on Ly^. Solutions for Q and T are then obtained via Eqs. (|2.9U|) and (|2.87|) . 



2.4.2 The Baryon Density pB 

Neglecting the wall curvature and assuming a step-function profile for the weak sphaleron 
rate, the baryon density satisfies the equation [491 165j : 



-z) —T^sriLiz) , 



(2.93) 



where ni;' is the number of fermion families and the relaxation term is given by [HSl '■ 



n 



1 + — + - 
6 / 2 



(2.94) 



where rigq indicates the number of light squark flavors, and the weak sphaleron rate is given 
by Fws = GKalT, with k ~ 20 jHS]- 

The solution to Eq. (|2.93|) in the broken phase, eventually growing into the Universe, is 
constant and given by: 



PB 



2v. 



(2.95) 



w J — oo 



Neglecting leptonic contributions, ni \s given in the unbroken phase by the sum of left- 
handed quark densities over the three generations (Qil, Q2L1 Q)- Since appreciable densities 
of first and second generation quarks are only generated via strong sphaleron processes, 
it is possible to express Qil and Q2L in terms of Q and T, in such a way that ul = 
Q + QiL + Q2L = 5(5 + 4T praj. Using then Eq. one obtains : 



-H 



ri + r2 



FssD 



1 



D 



(2.96) 
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where 



n 



r2 



dkqkT - ^kqkB - SkxkB 
kHi9kQ + 9kT + kB) 

kl{5kQ + AkT){kQ + 2kT) 
kH{9kQ +9kT + kB)^ 



(2.97a) 
(2.97b) 



and finally, in the broken phase: 



Pb{z > 0) = 



A 



r 

WS 


-^^ 


(- 




r 

J- ss 




D 


r 

WS 


vl 




D, 


r 

J- ss 


D 


D 



2D 



D 



+ 2nD 



vl + 1Z{D + D, 



(2.98) 



where the second equality is true in the limit ^ ADqTZ, which holds for the parameters 
we have chosen in this calculation. The contribution from the first term in Eq. (|2.98|) 
is linear in v^, due to the linear dependence on contained in the (3 appearing in the 
CP-violating sources. The second term is suppressed by two additional powers of Vw and 
generally leads to a negligible contribution to pB in the MSSM case (see discussion below). 
It could, however, be dominant in the case of heavy degenerate ti and t/j, which leads to 

n ~ [in]. 



The central feature emerging from the above discussion is that the net baryon density 
is proportional to „4 ~ S^/\/T. A large relaxation rate T for the relevant charges will 
suppress the overall baryon asymmetry. While in Refs. |2U1 162j it was pointed out how 
a non-equilibrium quantum transport could result in a resonant enhancement of we 
observe here that similar resonance effects in the relaxation terms will mitigate the impact 
of the enhanced sources. In the next section we discuss the numerical impact within the 
MSSM, but caution that reaching definitive conclusions will require computing the other 
transport coefficients, such as Fy, within the same framework. 



2.5 
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Baryogenesis and Electroweak Phenomenology within the 
MSSM 



The results derived in the previous Sections allow us to perform an illustrative, prelimi- 
nary analysis of baryogenesis within the MSSM. This should be taken as an exploration, 
whose robustness will be tested once we implement the next steps in our treatment of 
the source terms in the transport equations. With this caveat in mind, we explore the 
connections between electroweak baryogenesis and phenomenology within the MSSM, fo- 
cusing in particular on the implications for EDM searches. Throughout, we assume — as in 
mSUGRA — that all the terms in the Higgs scalar potential and all gaugino masses are real, 
while all the ^-parameters (trilinear scalar couplings) are equal at the GUT scale, therefore 
sharing the same phase (pA- In this case, the baryon asymmetry and EDMs are sensitive to 
the two independent CP violating phases (j)^ and (pA-^ 

2.5.1 Dependence of the BAU on MSSM Parameters 

From the structure of Eqs. (|2. 9812. 92)1 and (|2. 5212 .701) we can write the baryon-to-entropy 
density ratio^ Yb = Pb/s as: 



where we have isolated the dependence on the phases 4' A- The first term that 

contains Fi stems from the Higgsino source, while the F2 term arises from the squark 
source. 

The functions Fi and F2 display a common overall dependence on bubble wall parame- 
ters {vw, Lm, A/3), while having distinct dependence on other MSSM mass parameters such 
as \fi\, the soft mass parameters for gauginos (Mi^2) and squarks {M^^, M^^), the triscalar 
coupling \At\, and tan/5. In order to assess the size of Yb and the impact on CP-violating 
phases, we must choose a reference region in the MSSM parameter space, and we follow two 

*One may, of course, work with a more general soft SUSY-breaking sector that contains additional CP- 
violating phases. 

®We evaluate the entropy density at the electroweak phase transitions via s — (27r'^)/45 x gcs{T)T^, with 
(/off = 130.75, resulting in s = 57.35 T''. Similarly, to convert the present ratio ps/rij to Yb, we use the 
relation s — 7.04 n^.. U None of these densities, of course, exhibits an equation of state with w = p/ p < — 1 




(2.99) 
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Figure 2.5: Left panel: CP- violating Higgsino source Sp^ = —S'^/{v'^$sin<f>n), as a func- 
tion of \fi\. Right panel: relaxation rate Rr = {Th + T^j)/{rh + 'i^m)H.N., normalized to the 
value used in as a function of We have taken M2 = 200 GeV, and the values of all 
other parameters as indicated in the text. 

obvious guidelines: (i) we require that v{Tc)/Tc > 1, so that the baryon asymmetry is not 
washed out in the broken symmetry phase; (ii) we require no conflict with precision elec- 
troweak physics and direct collider searches. Both criteria lead to non-trivial restrictions. 

The condition of a strongly first-order phase transition [v{Tc)/Tc > 1] requires light 
scalar degrees of freedom coupling to the Higgs sector. It has been shown j22(l23j that within 
the MSSM the only viable candidate is a light top squark, which should be mainly right- 
handed (tji) in order to avoid large contributions to the p parameter. Quantitatively, for 
lightest Higgs boson mass rrih < 120 GeV, one needs 100 GeV < < rrit, and sufficiently 
small stop mixing parameter \At — /i/tan/3| < 0.6M^^ '22 . Moreover, present experimental 
limits on ruh and the constraint v{Tc)/Tc > 1 jointly require either values of tan/3 > 5 or 
M^^ in the multi-TeV region ^J. Based on these considerations, for illustrative purposes 
we work with the following values of MSSM parameters at the electroweak scale: M^^ = 0, 
M^^ = 1 TeV, I Ail = 200 GeV, M2 = 200 GeV, tan/3 = 10. We also take for the CP-odd 
Higgs mass = 150 GeV, which translates into A/? ~ 0.015 |68j . We vary in the plots 
the scale \p\, in order to display the resonant behavior for |/i| ~ M2. Finally, for the bubble 
wall parameters we adopt the central values = 0.05 and = 25/T |68j . 

With the above choice of MSSM parameters, the stop-induced contribution to Yb is 
suppressed (P2 ~ 10~'^Pi), since one is far off the squark resonance [{M^^ — M^^) ^ ^ijil- 
On the other hand, the Higgsino-induced contribution Pi can account for the observed Yb 
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Figure 2.6: Left panel: Higgsino contribution to Yb {Cf Eq. ()2.99() '). normalized to the 
observed value. Fi displays residual resonant behavior for ~ M2. All other input 
parameters are given in the text. Right panel: Stop contribution to Yb { Cf Eq. 1)2. 99(1 ') 
normalized to the observed value. The upper curve is for = M^^, while the lower one 
is for M^^ > M^^. We have taken here M^^ = 100 GeV, \n\ = 200 GeV, and have allowed 
to reach unrealistically low values to explore the size of the squark resonance. For 
realistic input parameters F2 <^ Fi. 



even without maximal values of | sin (pf^\. We highlight below the salient results of our study: 

• The primary result of our analysis is that both the source and the relaxation term 
F/j display the resonant behavior |20[ I62j typical of quantum transport for \fi\ ~ M2. 
We illustrate this in Fig. 12.51 the left panel shows the behavior of the rescaled CP- 
violating higgsino source Sjj = — 5^/(f^/5sin(/)^) versus while the right panel 
displays the ratio i?r of the relaxation term (Tfi + F^^) as calculated in this work 
to the one used in previous studies, (F/j + T~^j)b.n. QH]- To our knowledge this 
is the first explicit calculation showing resonance behavior for the relaxation term 
f ~rr(F;, + F^). 



Since Fi is proportional to S'?f'/ + the baryon asymmetry retains a resonant 



behavior, albeit with an attenuation of the peak due to the enhanced relaxation term. 
This is shown explicitly in Fig. 12.61 In the left panel we plot Fi/Y^^^^ , normalizing 
to the baryon asymmetry extracted from CMB studies [2j: Y^^^^ = (9.2 it 1.1) x 
10"" (the quoted error corresponds to 95% CL). 

For completeness we also display in Fig. 12.61 (right panel) the behavior of the squark 
contribution F2/Y^^^^ as a function of M^^, with M^^ = 100 GeV. Within the 
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Figure 2.7: Allowed band in the | sin plane, obtained by requiring successful elec- 

troweak baryogenesis. All other MSSM parameters are given in the text. The light-shaded 
(green) narrow band corresponds to the experimental input from WMAP, while the two 
bands combined [dark (blue) + light (green)] correspond to input from Big Bang Nucle- 
osynthesis. 



MSSM, precision electroweak data and the requirement that v{Tc)/Tc > 1 force the 
masses to be far away from the peak region. However, in extensions of the MSSM 
where the phase transition is strenghtened by additional scalar degrees of freedom this 
contribution might be important (see, e.g., Refs. j44| I45j ) . 



• For given values of the MSSM parameter space explored here, successful EWB carves 
out a band in | sm(j)^\ centered at | sin(/>^| = Y^^''/|i<i| (whose width depends on the 
uncertainty in Y^"^). Due to the resonant behavior of Fi, the location of this band 
is highly sensitive to the relative size of M2 and As illustration, in Fig. 12.71 we 
plot the allowed band in the | sin i;^>^|-|^| plane determined by the baryon asymmetry, 
with all other MSSM parameters fixed as above. The bands in the plot combined 
together correspond to the baryon density determined from Big Bang Nucleosynthesis, 
yBBN ^ 3 ^ 2.5) X 10"^^ (the error corresponds to 95% CL Using WMAP 
input leads to the narrow, lighter-shaded band in our plot located at the upper edge 
of the BBN-induced band. 
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Figure 2.8: SUSY loop graphs inducing quark EDM. A quark may develop an electric dipole 
moment through one-loop effects involving superparticles with CP-violating couplings at 
one of the vertices. The external photon sees one of the charged particles in the loop. The 
vertex couplings may involve the phases (j)A in the MSSM. 



2.5.2 SUSY-induced EDMs 



We conclude this investigation with a brief account of the connections between the baryon 
asymmetry and EDM phenomenology. Since the Standard Model predictions are in general 
highly suppressed and well below present experimental sensitivity, limits on the electron, 
neutron, and atomic EDMs can be used to constrain the phases of a given new physics 
model. 

Let us first review how an EDM can be generated. The EDM of a spin-^ fermion / is 
the coefficient df in the low-energy effective Lagrangian^*' : 

Ce = -^d/V^a^.TsV'P^^ (2.100) 

Such a term could be induced by one-loop diagrams in SUSY containing a CP- violating 
coupling in one of the vertices, such as those in Fig. 12.81 For an elementary particle such 
as the electron, it is only the coefficient d/ in Eq. (|2.100|) which gives rise to the measured 
EDM. However, for composite particles such as the neutron or an atom, the EDM dn or d^ 
of the whole particle may arise from not only the EDMs of the constituent particles, but 
also other CP-violating operators involving these constituents. Quarks, for instance, may 
develop a chromoelectric dipole moment d^ : 

C^c = -Id^q'T'^V^qG^u, (2.101) 



The notational conventions for the terms in the effective Lagrangian are those of Ref. 
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Figure 2.9: SUSY loop graphs inducing quark chromo-EDM. A chromoelectric dipole mo- 
ment of the quark may be induced by SUSY loops containing squarks with CP-violating 
couplings to the other particles. 

or a gluonic dipole moment d'^: 



Co = —d^f^^^G^^G^PG^'^^'e^"'^'', (2.102) 



which is called the Weinberg operator j20i- The chromo-EDM can be induced by the MSSM 
through the graph in Fig. 12.91 



2.5.3 Combining Constraints from the BAU and Electric Dipole Mo- 
ments 

The present experimental limits on EDMs of interest to us are those for the electron |71j . 
neutron [72|, and ^^^Hg [7^1 EDMs: 

141 < 1.9 X 10~^'^e • cm 
\dn\ < 7.5 X 10~^^e-cm 
IdHgl < 2.1 X 10~2^e • cm, 

all given at 95% CL. New experiments under development promise to improve upon these 
bounds by up to two orders of magnitude or more. For example, an experiment using PbO 
molecules at Yale may achieve a sensitivity of ~ 10^^^ e • cm for the electron EDM, while 
an experiment at Los Alamos may reach ~ 10"'^^ e • cm. Meanwhile, another experiment 
at Los Alamos using ultra-cold neutrons may test the neutron EDM with a sensitivity of 
~ 10^2^ e • cm." 

Although a single EDM can be sufficiently small even for maximally large CP- violating 
^^These and other present and proposed EDM limits are reviewed in Ref. |74|. 
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phases (due to cancellations), constraints from more than one EDM can be very powerful. 
In Ref. [73], for example, it was pointed out how limits on electron and ^^^Hg EDMs single 
out a well defined region in the (p^-cpA plane, for given values of gauginos, squark and 
slepton masses. As shown above, for each point in the MSSM paramter space, electroweak 
baryogenesis also selects a band in the 4>^-(j)A plane. This implies in general non-trivial 
constraints on the MSSM parameter space, as the EDM-allowed region need not in general 
coincide with the one required by the baryon asymmetry. 

To illustrate this situation, we have evaluated the bands in (j)^-(j)A allowed by present 
limits on electron, neutron, and mercury EDMs, and EWB for several representative points 
in the MSSM parameter space (see Figs. I2.1UI and 12.111) . In our analysis we take the ex- 
pressions for the electron EDM and quark chromo-electric dipole moments from Ref. f5!?| . 
In relating the ^^^Hg EDM to the quark- level CP- violating couplings, we follow the treat- 
ment of Ref. l75j where it was shown that the dominant contribution arises from the 
chromo-electric dipole moments of quarks {dq) according to 



For the neutron EDM, QCD sum rule techniques were used in Refs. [77j to derive the 
expression in terms of quark EDMs and chromo-EDMs: 



There are also in general contributions from the Weinberg operator and also four-quark 
operators, but Ref. [7S] demonstrated that, in the MSSM with large tan/3, these contribu- 
tions are only about ~ 10% the size of those in Eq. 1)2.104(1 . In this work, since we will take 
tan (3 = 10, we ignore these extra contributions for our purposes. We also neglect the renor- 
malization group evolution of (j)^ and (j)A from the weak scale to the atomic scale, having 
assumed a common, flavor-independent phase for the tri-scalar coupling at the former. 

The plots in Fig. 12.1(71 correspond to taking the first and second generation sfermions, 
along with the gluinos, to be degenerate with masses equal to 750 GeV; the gaugino mass 
Ml = 100 GeV; and the triscalar coupling A = 200 GeV. We consider then values of M2 
and \^\ corresponding to the peak of resonant baryon generation, M2 = |//| = 200 GeV. In 
^^For a recent reanalysis of hadronic EDMs in SUSY see Ref. |76|. 




(2.103) 



dn = 0.7{dd - 0.25(i„) + 0.556(7, (drf + 0.5(i„). 



(2.104) 
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Figure 2.10: Allowed bands in the (p^-cpA plane implied by consistency with the 95% C.L. 
limits on current and proposed limits in EDMs. In the left panel, we use the current 
limits on the electron, mercury, and neutron EDMs. In the right panel, we illustrate future 
limits that could come from improved sensitivities to electron and neutron EDMs. These 
constraints are shown together with the phases required by baryogenesis. The shaded [dark 
(blue) and light (green) combined] EWB band corresponds to BBN input 0, while the 
narrow light-shaded (green) band on the left corresponds to WMAP input [7. In these 
plots we use \fi\ = M2 = 200 GeV (resonance peak). The the other supersymmetric masses 
are as specified in the text. 

the left panel we show the current EDM constraints on the phases (p^^A^ and in the right the 
limits that could come from proposed electron and neutron EDM experiments. For these 
choices of MSSM parameters, Eq. (|2.99j) predicts for Yg: 

M2 = |)u| = 200 GeV : = -1.3 x 10"^ sin 0^ + 1.7 x 10"" sin(</.A + 0^). (2.105) 

In the left panel of Fig. 12. Ill we consider lowering the masses of heavy sfermions and gluinos 
to 500 GeV, which does not change Eq. (|2.1()5|) for Yg, but tightens the EDM bands. In 
the right panel of Fig. 12.111 we move the heavy sparticles back to 750 GeV, but move off 
the peak of resonant baryon production, to M2 = 200 GeV, = 250 GeV. In this case, 

M2 = 200 GeV,|/^| = 250 GeV : Ys = -2.0x 10"^ sin 0^+4.6x 10"" sin((?:)A+</'/.). (2.106) 
In both cases in Fig. 12. IH EDM constraints already rule out successful baryogenesis for the 



61 




-D.04 -a. 02 0.02 Q.04 -0.04 -0.02 o!o2 0.04 



(pti (rad) (j)fj, (rad) 

Figure 2.11: Choices for MSSM parameters ruled out by EDM constraints and requiring 
successful electroweak baryogenesis. Left panel: EDM and EWB bands in (pfj.-4>A plane on 
the resonant peak M2 = = 200 GeV, but with first- and second-generation sfermion and 
gluino masses at 500 GeV. The EDM limits are tighter in this case, ruling out successful 
baryogenesis. Right panel: For heavy sparticles at 750 GeV, but with M2 = 200 GeV, = 
250 GeV, which lies off the resonant peak for the creation of the baryon asymmetry. The 
required phases are already ruled out by present EDM limits. 



chosen parameters. Thus, we find the general trend that for M2 ~ |/i| (and some relatively 
heavy sparticle masses), electroweak baryogenesis requires only small phases, consistent 
with the constraints from EDMs. As one moves off resonance (or lowers heavy sparticle 
masses), then larger phases are needed to generate the observed baryon asymmetry, and 
this requirement tends to conflict with the EDM constraints. 

Ultimately, if supersymmetry is discovered at collider experiments, spectroscopy will 
dictate the input for mass parameters. Then joint constraints from low-energy EDM mea- 
surements and collider searches could be used to tightly test the scenario of baryogenesis 
at the electroweak scale, especially as a new generation of lepton, neutron, and neutral 
atom EDM searches will likely tighten the constraints in Fig. 7 by two or more orders of 
magnitude (for a recent discussion, see Ref. jJl] and references therein). 



2.6 Conclusions 
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It is instructive to consider the essential physics leading to the enhanced sources and re- 
laxation terms discussed in this work. The propagation of quasiparticles in the plasma is 
modified by scattering from the spacetime varying Higgs vevs that causes transitions to 
intermediate states involving other quasiparticle species. The system retains some memory 
of each scattering due to the presence of thermal widths, Tj, that reflect the degeneracy 
of states in the thermal bath. For Fj = 0, the oscillating exponentials appearing in the 
Green's functions wash out any memory of the scattering. For Fj =^ 0, the Green's functions 
now contain decaying exponentials as well as oscillating terms, and the memory washout 
is incomplete. The impact of quantum memory effects are, thus, characterized by the ratio 
of time scales, Tint/rp ~ Fj/wj, where Tint is the characteristic propagation time associated 
with a quasiparticle of frequency coi and Tp ~ 1/Fj, the plasma time, is the time scale on 
which transitions between the quasiparticle and other, degenerate states may occur. To the 
extent that the quasiparticle thermal mass and/or three- momentum is large compared to 
Fj, this ratio T\at/Tp is 0{e). 

A special situation arises, however, when the spacetime variation of the Higgs vevs is 
gentle and the thermal mass of an intermediate state is close to that of the initial state. 
Under these conditions, the scattering event injects essentially zero four-momentum into the 
initial state leading to resonant production of the intermediate state j. The characteristic 
lifetime of the latter is no longer Tint ~ l/i^ij but rather the resonance time scale 

^ (2.107) 



Aa;2 + Tl 



where Aw = — ujj and Fjj = Fj + Tj [see, e.g., Eqs. 1)2. 128112. 13U|) and 
Appendix 12. Bj . In this case, the impact of quantum memory is characterized by the ratio 
Tres/Tp. For \Auj\ <^ Tij, this ratio becomes of 0(1), and the impact of quantum memory 
is resonantly enhanced^^. On the other hand, for |Aa;| S> Fjj, the ratio is 0{e) and one 
returns to the more generic conditions. 

In this study, we have shown how this effect can enhance both the particle number- 



^''An examination of Eqs. 12.12812.1301 and 12.13212.13411 of Appendix 12 . B I indicates the presence of an 
additional, dynamical enhancement factor ~ uj/./Auj^ + T^- in the relevant integrals. 
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changing relaxation terms as well as the CP-violating sources that enter the transport 
equations relevant to electroweak baryogenesis. Importantly, the effect of resonant relax- 
ation tends to mitigate the impact of resonantly-enhanced sources, as both enhancements 
occur under the same conditions for the electroweak model parameters (in this case, those 
of the MSSM). We suspect that analogous resonant effects occur in other transport coef- 
ficients, such as the Ty Yukawa terms discussed above, but that the conditions on model 
parameters leading to enhancements — owing to simple kinematic considerations — will be 
different. It may be, for example, that the Yukawa interactions are no longer fast compared 
to the Higgs vev induced transitions when the latter are resonantly enhanced, and in this 
case, the solution to the differential equations will differ from the general structure obtained 
here and by other authors. This possibility is one that should be explored in future work. 

Additional refinements of the present analysis are clearly in order, including some form 
of all-orders resummation of the Higgs vev insertions (possibly along the lines proposed in 
Refs. [481 161j ') and a treatment of the axial charge transport equations via Eq. (|2.23|) . In 
principle, one would also like to study the density dependence of the thermal frequencies and 
widths, the impact of nonzero gaugino densities, variations in bubble wall geometry, and 
possibly higher-order effects in e, such as the departure of 6f of the thermal distribution 
functions from their equilibrium values. In short, it is apparent that EWB is not yet a 
solved problem, but rather one that calls for additional study. 

Undertaking this effort will be important for electroweak phenomenology. As illustrated 
here as well as in other studies {e.g., |2H1)) determining the viability of EWB within a 
given electroweak model involves a detailed interplay of collider phenomenology, precision 
electroweak data, EDM searches, and a careful treatment of the dynamics of the electroweak 
phase transition. In particular, in light of the open questions pertaining to the latter, it 
is too soon to draw definitive conclusions about the implications of the next generation 
of EDM searches for the baryon asymmetry. One hopes, however, that by the time these 
searches obtain their first results, the context for their theoretical interpretation will have 
been further clarified. 
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2. A Appendix: Propagators at Finite Temperature and Den- 
sity 



In this section, we derive some useful properties of propagators at finite temperature and 
density, using derivations based on those for the case of finite temperature and zero density 
in Refs. [MIES]. 



2.A.1 General Structure of Fermion Propagators 

We begin with the spectral function for fermions at temperature T = 1//3 in the presence 
of a chemical potential fi: 

Papix) = ixr [e-^(H-MN)|^^(^)^^^(o)|J ^ (2.108) 
where Z = Tr [e~^('^~^'^)]. It is convenient to define the retarded and advanced propagators: 

S^{x) = 9{x^)p{x) (2.109a) 
S^{x) = -e{x°)p{x), (2.109b) 

supressing spinor indices. The Fourier transforms of S^'^{x) and p{x) are related by: 

S-^t»,k) = ,/_"^-jje(5;iM_. (2,110b) 

It is possible to express the momentum-space spectral function in terms of a single product 
of ipa{x) and 'ipp{x) instead of the anticommutator in Eq. H2.108|) . whose Fourier transform 
is: 



9«^(u;,k) = j d4a;e*M-k-x)p^^(i,x) 



(2.111) 



X n) . 
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Now insert a complete set of states between the fermion fields: 



Pal3 



(^,k) = j d^xe^^-'-^-^^YM^ 



(2.112) 



We can rewrite the second term by switching summation labels and translating from x 
to 0: 



^g^(E„-E,)^g-^(k„-k,).Xg-/3E,g/3/.(iV„+l)^^|^^(Q) |^.^ 01V'/3(0) |n) , 



(2.113) 



which after integrating in Eq. ()2.112() . becomes 



i 5](27r)^<5(a; + - E,)6^{k + k„ - k,)e-^(^"+-)e^^(^"+^) (711^^,(0) | j) H 



(2.114) 

where we used the first delta function to replace Ej with En + a; in the exponential e~^^J . 
This can now be written: 



f ^4^giM-k.x) ^(n|e-^(H-^^)Va(x) |j) (j|V^;3(0) |n) (2.115) 



which is e ^^'^ times the first term of Eq. 1)2.112(1 . so we conclude: 



/9a/3(w,k) 



1 + e 



j (fxe'^'''-^-''^^Ti[e~'^'^^~>'^^^{x)i}p{Q)]. (2.116) 



Similarly, we could have manipulated the first term of Eq. (|2. 112(1 in the same way, and 
derived the companion relation: 



1 + e''^(^-^) 



(2.117) 
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Appearing on the right-hand sides of Eqs. H2.116I2.117() are the Green's functions S^{k^, k) 
and — S'^(A;'^, k), giving the relations: 

S>(A;°,k) = [1 - npik^ - fi)]p{k^ ,'k) (2.118a) 

S<(fc°,k) = -npik^ - fi)p{k^,'k), (2.118b) 

where nF{x) = 1/(1 + e^'). 

The various Green's functions satisfy the identities: 

S\x, y) = 5^(x, y) + S<{x, y) = S^{x, y) + S>{x, y) (2.119a) 
S\x,y) = S>{x,y) - S''{x,y) = S<{x,y) - S^{x,y) , (2.119b) 

which follow directly from the definitions in Eqs. H2.1UI2.1U9() . Thus, using Eq. ()2.118() . 
the time- and anti-time-ordered propagators can be expressed in terms of the retarded and 
advanced propagators: 

5*(A:0, k) = [1 - npik"^ - ^)]5^(A;°, k) + n^(A;° - ^)S^(A;°, k) (2.120a) 
5*(/t°, k) = -npik^ - fi)S^{k^, k) - [1 - npik^ - p)]S'^{k^, k) . (2.120b) 

Also note that p = - S"^ = S> - S< . 
2. A. 2 Bosonic Propagators 

Similar results may be derived from scalar bosonic propagators, for which the analog to 
Eq. is: 

G>(A;°,k) = [l + nB(A;°-^)]/>(A;°,k) (2.121a) 
G<(fc°,k) =nB(fc° -^)p(fcO,k), (2.121b) 

where the momentum-space spectral function p{k^,]<.) for bosons is the Fourier transform 
of: 

pix) = ItV {e-^(H-MN)[^(^)^^*(o)]| (2.122) 
The bosonic propagators also satisfy the identity p = — = G^ — G^. 
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2. A. 3 Tree-Level Propagators 



At tree level, the propagators S^'^ for fermions are given by: 



i{lt + m) 



(2.123) 



and G^'^ for bosons are given by: 



G^'^(A:0,k) 



(2.124) 



(fcO ± ie)2 - El ' 



where -E^ 



2 9 

k| + m . Note that these propagators are independent of the temperature 



and chemical potential, which only enter in the thermal distribution functions appearing in 
the relations of the retarded and advanced propagators to the other Green's functions, for 
example, in Eq. (|2.12U|) . 



2. A. 4 One-Loop Corrections to Massless Fermion Propagators 

Resumming the one-loop self-energy into the fermion propagator at finite temperature 
changes the pole structure of the propagator dramatically, introducing a new collective 
"hole" excitation of the plasma [531 154j . In fact, this structure can be shown to hold even 
beyond perturbation theory [SS]. Extending the results of Ref. [HSI to include dependence 
on a chemical potential, the propagator takes the form given in Eqs. H2.31fE3n|) . Recall 
that in those equations Sp^h = ^p,h — i^p,h are the complex poles of the spectral function, 
and Zp^h are the corresponding residues. At leading order in the "hard thermal loop" ap- 
proximation (see Ref. 51 ), calculating the poles only to order £ ~ gT, one finds F = 0, 
and Zp h{k,fj,) and ujp^hik, fJ-), where k = |k|, depend only quadratically on fi/T, which we 
thus neglect in our analysis in the present work, where we keep only effects linear in n/T. 
In this limit, and including only a single gluon loop in the quark self-energy diagram, the 
poles of the spectral function are given by the solutions to the equation: 




(2.125) 



68 

where Cp = 4/3 is the Casimir of the fundamental representation of SU{3). The solutions 
to this equation give the poles = Ep{k), —Eh{k). The residues satisfy: 

Zp,h{k) = ^^'^ ~ ^ , (2.126) 



where 

~2 



m) = ^ . (2.127) 



Calculation of the imaginary parts ^-p^h of the poles, since they begin at order g^T, requires 
a resummation of hard thermal loops in self-energy diagrams j56| I57| I58j . We are also 
interested in their dependence on the chemical potential [i. We leave the calculation of 
these effects to a future study. 



2.B Appendix: Expanded Source Terms for Quantum Trans- 
port 



2.B.1 Bosons 



The CP-conserving source term for right-handed stops in Eq. 1)2. 46p can be expanded by 
explicitly taking the imaginary part of the integrand: 



1 iVcyi 



T 27r2 



\AtVu{x) - ii*Vd{x)y 



OO 7,2 



k^dk 

^Ri^L 



(2.128) 



' 1 1 

— (sin cl)lmh^ + cos (p Re /i^j) — - (cos Re /ij — sin 6 Im h~^) 

— (sin (plmh'l^ — cos Re /i^) -|- - (cos 9Keh^ — sin 9 Im /i^) | 
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where 



uoL,R = ^\\^V + Ml^ (2.129) 
A = ^{Tl + VrY + {UJL - ojrY 



tan 6 
tan 6 



{eM{^L,R±irL,R)/T]-iy 



and where Tl^r are the thermal widths for the tL,_R. The rates defined in Eq. 1)2.49(1 can 
then be expressed: 

rf = -^1^ l^tM^) - f^*Mx)\' (2.130) 



X / < — - [sin(/)Im(/i^ ± /i^) - cos<?!>Re(/i^ =F ^J)] 

+^ [cos e Re{hl =F h^) - sin 6 Im(/i^ =F h^) 

Meanwhile, the CP-violating source given in Eq. (|2.52|) can be expanded: 

sf{x) = NcynMf-^AtHxfm (2.131) 

^ [Re (l + n+ + n+) sin 29 + Im (n+ + n+) cos 261] 

+ ^2 (""R ~ ""l) ^^'^ '^^ ~ (^i? + ""l) ^'^l I ' 

where nf^ r = nsiuJij^ ^±.Tl^r). Our result agrees with that of Ref. except for a different 
relative sign in front of the cos2(/) term and the overall factor of Nc- 



2.B.2 Massive Fermions 
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The CP-conserving rates for Higgsino-gaugino interactions given in Eq. (|2.66|) can be ex- 
panded: 



k^dk 



27r2 \uj^^^ 



1 



(2.132) 



A 



[w^cj^ + r^r^^ - k"^ + M2\n\ cos<?!>^sin2/3(x)] 



where 



cosdiReih— ^ ht) — smchlmih— ± ht) 



sin 6 Ref/ii =F ^ ~ ) + cos d> Im( /ii it ] 



'^\\[^H^W ~ ^H^W + ~ l/^l cos(/)^ sin2/3(x) 



cos e Re(/ii =F ~ ) - sin 6* Im(/i± =F ^ ~ ] 



[^H^w + ^vK^//] [cos 9 Im(/i± ^ /i± ) + sin 6 Re{h± ^ h±] 



A 

5 

tan 
tan 6 



W,H 



w 



T <^T~ 



(2.133) 



The CP-violating Higgsino source in Eq. 1)2. 70p can be expressed: 



sf^{x) = 2glM2lm{ti)v{xffi 



k'^dk 



1 

A2 

1 

^ L' 



sin2(/> Re ( A^i - iV± ) + cos 26 Im ( iV± + 

W HI ^ \ W H 



sin 26* Re 1 - iV- 



w 



Nt 

H 



cos 29 Im iV+ + Nt 



(2.134) 
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where ~=nB(w^^^± iF^ Our result agrees with that of Ref. PU] except for the 
sign of the cos 2(/> term. 



2.B.3 Chiral Fermions 

For chiral fermions, the CP-conserving chirality-changing rates in Eq. 1)2.72(1 can be ex- 
panded: 



I ^ 1 Ncytvl 



,2 /-oo 



k'^dk 



(2.135) 



sin0p{Re(A^/i+^ T A^/i+j) - Im(/i+^ ^ /i+j)} + cos0pRe(/i+^ ^ 



+ f cos 20p(A^ T A^) Re(l - iV+^ - 7V+j) 

Or, 



^sm^hp{MX^ph;^ ± X^h+^) - lm{h;^ ± htj,)} - cos0,pRe(/i+^ T Kr) 



A 



■hp 



+ 



- COS 



^hp 



2</>,p(A^±Af)Re(iV+^-iV+j)' 



where 



^hp 

tan ^„ 



tan ( 



(o.,^ + u:^f + (F^ + F^)2 
{ul; - u^f + (Ff + F^)2 



P 

hpiuJp ± iFp ), etc. 



nF{u;p ± iFp ), etc. 



^P+^P 
F^ + F« 



-^P h 



(2.136) 



and where the 



A 



L,R 



L,R 

p,h 



(2.137) 



parameterize the linear shifts in the thermal widths due to non-vanishing chemical potential. 
As noted at the end of Appendix 12. Al in a fully resummed calculation of the fermion self- 
energy, such shifts which are linear in Hi/T may arise, and thus have to be included in 
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our calculations, which we defer to future work. Also note that we have approximated the 

L Ft 

residues Z \ to be purely real, which is true at the order we are working. 



2.C Appendix: Residues of Thermal Distribution Functions 

The expressions for the CP- violating and conserving sources presented in Sec. 12.31 omit the 
terms arising from the residues of the poles of the thermal distribution functions appearing 
in the thermal Green's functions. 

For example, the CP-violating source for squarks, Eq. (|2.52|) . arises from the second 
term in Eq. (|2.42|) . which at an intermediate stage in the derivation takes the form: 

X [nB(A;0) -nB(gO)]p^j(A:0,k)pi((?o,q). 

The hP,q^ integrals can be done by contour integration. The exponential factor 6*'-'^°"'^°^* 
determines that the contour should be closed in the lower half-place, while the contour 



should be closed above. The terms in Eq. (|2.52)) come from picking up the residues of the 
poles in the spectral functions />ij(A;'^,k) and PL{<f-,^)- However, nsik^) and nsiq^) also 
contain poles. The function nsix), 

= gx/T_i ' (2.139) 

has poles at the points = 27rmT, where n is any integer. These are illustrated in Fig. 12.121 
Near one of these points, the Bose-Einstein function behaves as: 

nB{x) = ^^ + --- , (2.140) 

X Xfi 

where the ellipses denote non-singular terms in the series expansion of nsix). The residues 
of these poles generate new contributions to sf^. First, note that the residues at A;'', g*^ = 
contribute nothing, as pj^(0, k) = pL{0,q) = 0, as seen by inspection of Eq. (|2.28j) . The 
remaining residues give, after completing the time and integrals in Eq. (|2.138|) . the 
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Figure 2.12: The integral in Eq. (|2.138|) is evaluated along the contour closed in the 
upper half-plane. The integrand contains poles at the locations from the spectral 

function and at x„ = 2mnT from the Bose-Einstein function nsiq^)- The other integrands 
appearing in the various source terms all have this basic analytic structure. 



extra contributions: 



-%lm{fiAt)v\X)${X)T 



dk r 



oo . 



n=l 



2ujR 



1 



{8l + 2i:inTY {£L-2mnTY 

1 1 
{£r - 2-KinTY ~ {£*j^ + 2mnT f 



PR{-2-ninT, k) 
PL{2TTinT, k) 



(2.141) 



In the region of resonant enhancement of the original terms in S-^ , where M^^ ~ ^ta^ 
these new terms are numerically about 10^ times smaller (primarily because of the large 
factors of T in the denominators), while they become comparable only far away from the 
resonant region. Since these regions are phenomenologically unimporant in our analysis, 
we can saiely ignore A^P^for our purposes, although they should be included in future 
calculations for complete consistency. 

Corrections to the CP- violating Higgsino source in Eq. come similarly from 

poles in the Fermi-Dirac function nF{x) at the points = (2n — l)z7rT for integers n: 



np{x) 



1 

e^/T + 1 



T 



+ 



(2.142) 
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The extra contribution to the source generated by the residues of these poles: 



A5g(X) = '^lmifi)viXfPiX)M,T J 



dk 



+ 



{£y/ l^nY {£— + l^r, 



('^^~ + /^ra)^ {£y^r 



PH±{f^n, k) 



(2.143) 



Again, the size of is numerically insignificant compared to given in Eq. (|2.7Uj) . 

except far away from the region |//| = M2, thus leaving our phenomenological studies 
substantially unaffected. 

The CP-conserving sources S*^^ are also modified. For instance, to Sf^ in Eq. 1)2. 46() 
must be added the terms: 



^Sg''(X) = \A,MX) - ^'vdXtT J dkk^ 



n=l 



p'^{—27rinT, k) 
PR{-2TrinT, k) 



1 



+ 



1 



£l - 2TTinT £1 + 2mnT 
1 1 



{£l - 2-KinTf {£1 + 2mnTf 



-ML 



n=l 



p'^{2TTinT, k) 
PL{2'irinT, k) 



1 



+ 



1 



£r — 2TrinT £^ + 2TrinT 
1 1 



(^K - 27rinr)2 (£:* + 2TTinTY 



(2.144) 



where 



d 



+ 



(2.145) 
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These corrections are also numerically insignificant compared to Eq. H2.46() . except far away 
from the region M^^ ~ -^ijj- Similar expressions should hold for the C-P-conserving sources 
for quarks and Higgsinos, which we also expect we can safely ignore, though we do not 
include their explicit expressions here. 

2.D Appendix: Towards the Yukawa Source 

In this Appendix we take initial steps towards a computation of the source terms arising 
from the Yukawa interactions illustrated in Fig. 12.31 These describe the interactions, for 
instance, of squarks with real Higgs bosons: 



^nt = vMAtHl - fi*H^/)i*ji + h.c 



(2.146) 



These interactions contribute to the squark source given by Eq. (|2.19|) : 



oo 



G<{X, z)T>{z, X) - G>(X, z)^<(z, X) 
+S>(X, z)G< (z, X) - S<(X, z)G>{z, X) 



by inducing the self-energies: 



(2.147) 



E>(x, y) = -y^G>{x, y) [\At\' G>o (x, y) + \^,\' G>o. (x, y)] 

" d 

S<(x, y) = -y^G<{x, y) [\At\' G<, (x, y) + j/ij^ G<o. (x, y)] . 



(2.148a) 
(2.148b) 



Note that these self-energies contain no C-P- violating phases, unlike Eqs. 1)2. 3912. 4U|) . which 
contained cross-terms between the Vu and Vd vevs. Thus, these Yukawa interactions involv- 
ing real Higgs particles contribute only to the CP-conserving part of the squark source. 
The Hu and Hd contributions to the self-energy have essentially the same structure, so for 
clarity, we include only the former in the following calculations. The contributions can 
be restored straightforwardly. 

Every Green's function appearing in Eqs. (|2.147l I2.148() contains dependence on the 
corresponding chemical potential //j. Expanding each one to first order in Hi/T, 



G^iX,z) = G''i{X,z)+^i^6G^{X,z). 



(2.149) 
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Let us Wigner transform each of the Green's functions to momentum space: 

G,iX,z) = I ^e-^^'-(^-^)G,(p;/.,(X)), (2.150) 

where we have approximated the dependence of the chemical potentials fj,i{X+z) on the col- 
lective coordinate X + zhy simply /u(X), which assumes the variation in is slow compared 
to the typical scale of individual interactions between particles. With this assumption, after 
plugging Eqs. H2.148|) into Eq. H2.147|) . using the Wigner transforms ()2.15U() . and expanding 
to first-order in the chemical potential /Uj/T as in Eq. H2.149() . the Yukawa- type source for 
squarks can be rearranged into the following useful form: 

^^^^^^ = jy Iwfl |^[e^^^°-°^* + e-(^»-°)*]/..(/) 

x|-^pMp°,p)[Sj^(^p)-S<(g°,p)] 

^ (2.151) 

+ f ^L(p°,p)[S^(?°,p)-S<(g°,p)] 
+ ^PH(Ap)[S^(g°,p)-S<(gO,p)]}, 

where we have used 

G>{p) = [1 + nB(/ - ^i^)]p^{p) (2.152a) 
G<{p) = nB(/ _ i^,)p,{p), (2.152b) 

so that 

5G>ip) = 5G<ip) = -^^hB{p')p^{p). (2.153) 

The momentum-space self-energies are evaluated at zero chemical potentials and are given 
by: 

^^G>{k)G>{q-k) (2.154a) 
^>{q) = -y^\Atf I ^^G>{k)G<{k-q) (2.154b) 

/H^k 
^^G>{k)G<{k-q), (2.154c) 
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where, again, all Green's functions here are evaluated with zero chemical potentials. The 
self-energies T,f are obtained by flipping all ^ signs. 

The evaluation of Eq. ()2. 151(1 is considerably complicated by the presence of finite widths 
Ti in the Green's functions appearing inside the integrand. Naive contour integration as 
for the FM-type sources derived earlier produces a result which is ultraviolet-divergent. 
Investigations into the proper regulation of these terms or a correct procedure for integra- 
tion (which is even further complicated by the poles in the thermal distribution functions 
as described in the previous Appendix) is still underway at the time of this writing. The 
quantitative analysis of Sj^ and the comparison of its size to the T jvf-type sources derived 
earlier is essential to check the consistency of the approximations used in solving the trans- 
port equations that give the left-handed weak doublet fermion density and, thereby, the 
baryon density, ps- If the coefficient Ty appearing in the transport equations is not con- 
siderably larger than T^^, then the assumption that Ty S> must be discarded, changing 
the solution of the transport equations entirely. This scenario is particularly likely to oc- 
cur in the regions of MSSM parameter space where Tjj is enhanced. The evaluation of 
the Yukawa-type sources and their impact on the phenomenological analysis presented in 
Sec. 12.51 is one of the most urgent tasks handed to us by the basic foundational analysis 
presented in this chapter. 
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Chapter 3 

Effective Theories of Strong 
Interactions 

Theoria bczeichnet die rein cmpfangendc, von aller >praktischcn« Bezweckung des 
tatigen Lebens durchaus unabhangige Zuwendung zur Wirklichkeit. Man mag diese 
Zuwendung »uninteressiert« nennen — wenn hiermit nichts anderes ausgeschlossen sein 
soil als jegliches auf Dienliclikeiten und Belange gcrichtcte Abschen. Im iibrigen ist 
hier auf hochst entschiedene Weise Interesse, Beteiligung, Aufmerksamheit, 
Zielsetzung. Theoria und contemplatio zielen mit ihrer voUen Energie dahin — freilich: 
ausschliefilich dahin — , dafi die ins Auge gefafite Wirklichkeit ofFenbar und deutlich 
werde, dafi sie sich zeige und enthiiUe; sie zielen auf Wahrheit und nichts sonst. 

Josef Pieper, Gliick und Kontemplation 

Theoria has to do with the purely receptive approach to reality, one altogether 
independent of all practical aims in active life. We may call this approach 
"disinterested," in that it is altogether divorced from utilitarian ends. In all other 
respects, however, theoria emphatically involves interest, participation, attention, 
purposiveness. Theoria and contemplatio devote their full energy to revealing, 
clarifying, and making manifest the reality which has been sighted; they aim at truth 
and nothing else. 

Josef Pieper, Happiness & Contemplation 

In this chapter we review the basic features of several effective theories for the strong 
interactions — heavy quark effective theory, soft-colUnear effective theory, and non-relativistic 
QCD — following developments of these theories in recent years. This forms the background 
for the discussion of the applications of these effective theories pursued in the subsequent 
chapters. 

3.1 Quantum Chromodynamics 

It is a lovely language, but it takes a very long time to say anything in it, because we 
do not say anything in it, unless it is worth taking a long time to say, and to listen to. 

Treebeard, in The Two Towers, by J.R.R. Tolkien 
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Quantum Chromodynamics (QCD) is the theory of the strong interactions between quarks 
and gluons, which bind together to make protons, neutrons, and other hadrons. The theory 
accounts successfuhy for many hadronic phenomena, especially those occurring at relatively 
large energies, for example, in the bombardment of protons by highly energetic electrons 
in deep inelastic scattering, revealing the pointlike substructure of the proton. QCD is 
governed by the Lagrangian density: 

Cqcd = - - Ig^uG^'^", (3.1) 

where we sum over quark flavors q and SU{3) generators A. The covariant derivative is 

QCD allows for reliable quantitative predictions at large energies due to the phenomenon 
of asymptotic freedom |79[ l8nj. The observed strength of the interaction between quarks and 
gluons is characterized by a coupling constant Qs, which is a function of energy, becoming 
small at large energies and large at small energies. Physical quantities can be calculated as 
a perturbation series in powers of Qs-, or, rather, = g1/A.-K, which is a reliable procedure 
as long as as <C 1. 

At small energies, however, such as the scale at which quarks and gluons bind together 
into light hadrons (protons, pions, etc.), the coupling constant is large, > 1, and per- 
turbation theory breaks down completely. Even in processes involving strongly-interacting 
particles at large energies, because the particles used or observed directly in experiments 
are hadrons, not free quarks and gluons, these low-energy binding effects contaminate the 
analysis, preventing completely precise calculation of the observable quantities. However, it 
is often possible to separate the perturbatively-calculable large energy phenomena from the 
low-energy hadronization effects in a way that preserves much predictive power. Imagine 
calculating the total rate for Z bosons to decay to hadrons, T{Z — > hadrons). In perturba- 
tion theory, we begin with the Z coupling to quarks: 

We can calculate the decay rate for Z to qq, T{Z qq). This process produces a qq pair 
moving back-to-back with energy Mz/2. This is much larger than the scale of hadronization. 
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^QCD- The quark and antiquark move far apart before they hadronize. We do not know 
how to calculate the dynamics of this hadronization, but since we know it must happen 
with probability 1, we can approximate the total hadronic decay rate of the Z by: 



to leading order in as(M^). We can compute to next order in as{Mz) by including the 
rate T{Z — > qqg). And so on. Because we have separated the parton-level physics from the 
longer-range, lower-energy physics of hadronization, we are able to make a reliable quanti- 
tative prediction that is a very good approximation to reality. We say we have "factorized" 
the total hadronic decay rate of the Z — into T(Z — > partons) and the probability for partons 
to hadronize, which is 1. 

For more complicated observables, we seek to perform similar factorizations, although 
the low-energy, nonperturbative quantities which appear will not, in general, be so simple 
as "1". These may be objects such as parton distribution functions, meson light-cone wave 
functions, etc. Although such quantities are not calculable in perturbation theory, they may 
appear in more than one physical observable, thus providing some remnant of predictive 
power. They may also be calculable numerically in lattice QCD. 

Proving factorization and identifying the relevant nonperturbative quantities is a major 
thrust in the direction of modern research in QCD. There are two main camps of research. 
There are those who attack the problem directly in QCD, the so-called perturbative QCD 
or QCD factorization approach Then there are those who quail at the enormity of 

full QCD and attempt to work instead with a simpler version of the theory. This is the 
approach of effective field theory [821 183j . which we adopt in this thesis. 

3.2 Example: Heavy Quark Effective Theory 

In an effective field theory one identifies a small parameter determined by the relevant 
physics in a given problem in which to expand the full theory to a given order of approx- 
imation. Equivalently one integrates out the degrees of freedom living at an energy scale 
much larger than those relevant in the physical problem at hand.^ 

^The presentation in this and the following section are heavily influenced by the lectures presented in 




(3.3) 



Ref. 
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As an illustrative example and a prelude to the effective field theories we consider later 
in this thesis, we overview the heavy quark effective theory (HQET) jSS], an effective field 
theory describing hadrons containing one heavy quark, namely, charm (c) or bottom (b)? 
In such a hadron, the heavy c or 6 quark moves with a relatively slow velocity compared to 
the other light quark(s) accompanying it. The typical momenta of the constituent partons is 
of order p ~ Aqcd- Thus, the velocity of the heavy quark Q is of the order v ~ Aqcd/'tt-q- 
We can expand the Lagrangian of full QCD in Eq. p.lj) for heavy quarks in powers of this 
small velocity v. 

Physically, we imagine that the interaction of a heavy quark with low-energy gluons 
causes only small fluctuations of its momentum, p = rngv + k, with k ~ Aqcd- We want 
the effective theory to describe these small fluctuations. First, we remove from the full 
QCD heavy quark field Q{x) the dependence on the large momentum p = rnqv, defining a 
new field Q^: 

Q(x) = j;e-*'"«'"^'Q,(x), (3.4) 

V 

summing over labels v. In terms of the new fields, the quark part of the QCD Lagrangian 
becomes: 

= Qv{x){ip + mQ]^ -mQ)Q^{x), (3.5) 
summing implicitly over u.^ Now write Qv{x) as the sum Qv{x) = h^{x) + Hy{x), where: 

hv{x) = ^-^Qv{x), H,{x) = i^Q,(x). (3.6) 

Then, the heavy quark Lagrangian (|3.5() becomes: 

-Cq = hv{x)iIpK{x) + H.f,{x){ilp - 2mQ)Hy{x) 

(3.7) 

+ H^{x)iIpK{x) + K{x)iIpH.,{x). 

Derivatives acting on the fields hy,H^ produce momenta of order Aqcd/w-Q; so the kinetic 
term for Hy{x) is suppressed relative the leading term quadratic in H^, —2mQHy{x)Hy{x) . 
therefore is not a dynamical field and can be integrated out by using the equation of 

^We limit our attention, not surprisingly, to heavy hadrons containing two (Qq) or three (Qqq), not five 
|86l 187) . quarks and antiquarks. 

^Strictly speaking, we should have summed over separate labels v,v' for the two heavy quark fields, but 
interactions with soft gluons cannot change these label velocities; hence we sum over only a single velocity. 
This simplification will not occur in SCET. 
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motion (obtained by varying with respect to Hv)'- 



{ip - 2mQ)H^ = -iJpK. (3.8) 

Substituting into Eq. p.Tj) . we obtain 

Co = hv{x)iphy{x) - K{x)ilp ^ ilj)hy{x). (3.9) 

lip - 2mQ 

Expanding in powers of l/mq: 

Cq = K{x) {ivD - I^JPP + • • • ) K{x), (3.10) 

having used P^^^Py = v^, where P„ = (1 + ^)/2. The dots denote higher-order terms in 
the l/mq expansion. 

Keeping terms only to a fixed order in l/mq in Eq. (|3.1U() defines the Lagrangian of 
heavy quark effective theory. The leading order term is very simple: 

CnQET = hvix)iiv-D)hy{x), (3-11) 

giving the heavy quark propagator 

2 vk + ie ^ ' 

and the heavy quark-gluon interaction vertex 

-igT^v^. (3.13) 

This leading-order Lagrangian exhibits more symmetries than evident in the full QCD La- 
grangian. First, it is independent of the heavy quark mass, giving rise to a flavor symmetry 
between h and c quarks. Second, containing no Dirac matrices, it also exhibits spin sym- 
metry. This heavy quark spin-flavor symmetry is thus an approximate symmetry of full 
QCD made evident only by the effective theory expansion in l/mg. It allows for powerful 
predictions for processes involving different hadrons containing heavy quarks which are re- 
lated by spin-flavor symmetry. Violations of this symmetry at higher order in l/mq can be 
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systematically calculated by keeping more terms in the HQET Lagrangian p.lUf) ^. 

This relatively simple example of HQET illustrates two key features we will seek in 
formulating any effective field theory: 

• Extra, approximate symmetries not obvious in the full theory but easily identified at 
leading order (s) in the the effective theory. 

• The simplification of interactions between heavy (hard) particles and soft ones, as in 

Eq. (nm|) . 

The first property grants us more predictive power in relating apparently different physical 
processes to one another, while more sophisticated exploitation of the latter will simplify 
proofs of factorization of perturbative and nonperturbative contributions to physical ob- 
servables. We now embark on this task in the effective theory truly of interest in this 
thesis. 

3.3 Soft-Collinear Effective Theory 

Consider a process in which there are hadrons moving with very large energy compared to 
their invariant mass, for instance, in the decay of Z bosons to light hadrons as considered 
in Sec. 13.11 or in B decays such as S — > Dtx or B ^ vrvr. We can formulate an effective 
field theory for quarks and gluons moving on such collinear trajectories and the soft partons 
with which they interact — the soft- collinear effective theory |89l l9Uj . 

3.3.1 Effective Theory for Inclusive Decays— SCETi 

We first identify the relevant energy scales and small parameters to use in formulating the 
effective theory. Recall that in HQET, we took this to be the heavy quark mass mq and 
velocity v ~ Aqcd/'"T'Q- Here, we begin by defining two light-cone vectors, n and n, which 
satisfy 

= n^ = 0, n-n = 2, (3.14) 
"'Or estimated by being more clever |88|. 



84 

for example, n = (1,0,0,1) and n = (1,0,0,-1). Vectors can be decomposed along these 
two light-cone directions and the orthogonal transverse directions: 

= ^ri'' + + v^. (3.15) 

We thus define the light-cone components V'^ = n-V, V~ = n-V. For a collinear particle, 
one light-cone component of its momentum will be large, the other small. Its momentum 
will scale in powers of some small parameter A as: 

(P+,P",P±)~Q(A^1,A), (3.16) 

where Q is the large scale governing the physical process being studied, on which the 
definition of A also depends. The collinear particles also interact with soft, or ultrasoft, 
particles, with momenta scaling as: 

soft: Ps ~ Q\ (3.17a) 
ultrasoft: p^s ~ QA^. (3.17b) 

The effective theory for collinear and (ultra)soft particles with these scalings is called SCETj. 
Another choice of scalings gives the theory SCETn, to be introduced in the next section. 



We form the SCETj Lagrangian essentially by expanding the QCD Lagrangian in powers 
of A. We start with the Lagrangian for collinear quarks, which we take to be massless. In 
QCD, 

Cq = qipq. (3.18) 

We seek to describe fluctuations in the momenta of the collinear particles about its collinear 
trajectory. So, for a collinear quark of momentum p, we write: 

p = p + k, (3.19) 

where p contains the large {label) components of the momentum, 

V^=r^+P±, (3.20) 
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and the residual momentum k ~ Aqcd represents the fluctuations about the label momen- 
tum. 



As in HQET, we extract the large momentum fluctuations from the full QCD field, 
defining a new field qn,p- 

q{x)=Y,^-'P"qn,p{x). (3.21) 



The QCD Lagrangian (|3.18|) then becomes: 



(3.22) 



Again as in HQET, we project out the large and small components of Dirac field, writing 
qn,p = in,p + Sn,p, where: 



in,p — . Q.n,pi ^n,p 



-Qn. 



p, 



(3.23) 



in terms of which the Lagrangian 1)3.22(1 becomes: 



p,p 



in,p'{x)-in-Din,p{x) + H„y (x)-(p" + in-D)'~.n,p{x) 



+ in,p>{x)(§^ + i^_L)S„,p(x) + H„y(x)(^_L + l_)in,p{ 



(3.24) 



Again, as in HQET, we have a derivative, n ■ d, acting on H„^p, which is suppressed relative 
to the term containing the label momentum p~ . So the field is not dynamical, and we 
eliminate it using its equation of motion: 



p +m-D 2 



Inserting this into (|3.24|) leaves us with the Lagrangian: 



(3.25) 



'i(p—p')-x 



P,P 



in-D + {p± 



1 



p +in-D 



in,p{^). (3.26) 



To sort out the interactions between gluons and collinear quarks, we must first distinguish 
between the collinear and (ultra)soft gluon fields. The gluon fields appearing in the covariant 
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p + k 



.Tfi n - p 

2 n ■ kfi ■ p + p'^ + ie 



P 



n - p n ■ p' n ■ p n ■ p' I 2 



V 



Figure 3.1: Feynman rules involving collinear quarks in SCETi. The collinear particles are 
shown with label momenta and residual momentum k. 



derivatives in the collinear quark Lagrangian ()3.26|) can be split up: 



A = ^ A' ^ A"" 



(3.27) 



where the collinear, soft, and ultrasoft fields are assigned the power countingas: 



^'=~Q(A^1,A), A^~Q(A,A,A), A'^^ ~ Q(A2, A^, A^^ 



(3.28) 



to match the scalings of the corresponding momenta. Note that interactions of soft gluons 
with collinear quarks leave the quark with the momentum scaling (5(A, 1,A), which does 
not exist in the effective theory. So soft gluons should not appear in the Lagrangian 1)3. 26() . 
Collinear quarks do interact with ultrasoft gluons. The effective theory collinear gluon fields 
q are defined by: 

A\x) = Y,e-^-^-^ K.,,i^x\ (3.29) 



factoring out the large label momentum q. 

In the Lagrangian ()3.26|) . the components Z)j_ and n ■ D of the covariant derivative 
contain both ultrasoft and collinear gluons. Due to the power counting in Eq. (|3.28|) . the 
ultrasoft gluon fields give a subdominant contribution compared to that of the collinear 
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gluons. Thus, the correct cohinear quark Lagrangian in SCETi at leading order in A is: 



|e„,p(x), (3.30) 



where = d — igA^ contains only the collinear gluon field, while n-D still contains both 
collinear and ultrasoft fields. Some of the Feynman rules arising from this Lagrangian are 
shown in Fig. 13.11 

The Lagrangian (|3.3()j) may be written in a still more compact form by the introduction 
of the label operators pjTj : 

n-Vin,p = n-pCn,p (3.31a) 
KUp = pIUp, (3-31b) 

and similarly for label operators acting on collinear gluon fields. In the second, messier, 
term of the Lagrangian (|3.3()|) . the ordinary derivatives which appear hit everything to their 
right. Thus, they bring down extra factors of the label momenta of the collinear gluons, 
because of the exponential factor in Eq. 1)3. 29() . For example, 

n-dAUx) = = ^ e-«'-n.g^^|g(x). (3.32) 

g g 

We can account for derivatives on the exponential factors and eliminate the explicit label 
momenta appearing in the Lagrangian 1)3. 30(1 by making use of the label operators in (|3.31|) : 



-CsCETi = 

where now 



in-D + Hpl-J—ilf)'^, 



|en(x), (3.33) 



Dl = V,-igA'^, (3.34) 

and 

en(x) = X;e-'^"en,p(x). (3.35) 
p 

Finally, note that this Lagrangian contains the collinear gluon field n-A'^ in the denominator 
of the second term, giving rise to couplings of collinear quarks to an arbitrary number of 
n- A'^ gluons. These infinitely many couplings can be resummed into the form of Wilson 



lines: 



Wn{x) = Pexp 



ig / dsn-A'^{x) 



where P denotes path ordering. Using the property 



(3.36) 



Wnix)n-VWl{x) = n-V - ign-A^ix) = n-D" 



(3.37) 



we can express the SCETj Lagrangian 1)3. 33() as: 



in-D + iIp\Wn{x)^—Wl{x)iIp' 



in-V 



(3.38) 



It is also possible to argue that invariance under gauge transformations of collinear field 
requires this precise form of the Lagrangian, which we do not go through here. 
Similar analysis leads also to the effective Lagrangian for collinear gluons j92j: 



where 



^(s) = ^ e-(^--^")-^Tr { \^V^ + gA^,,,, iV^ + 5^,,,] ]' 



72 T}/ 

= n-V— + + in-D—. 



(3.39) 



(3.40) 



Ghost fields and gauge-dependent terms have been ignored in Eq. 1)3. 39(1 . 



3.3.2 Decoupling Ultrasoft Fields 



We could take Eq. (|3.38|) as the final form of our Lagrangian for SCETj. However, one more 
manipulation simplifies greatly the separation of hard and soft (nonperturbative) physics 
in the effective theory |92j . First, introduce the ultrasoft Wilson lines: 



Yn{x) = Pexp 



ig I dsn-Aus{ns) 



Rewrite the fields in (|3.38|1 as: 



(3.41) 



e„(x) = y„(x)^f (x) 

An{x) = Yn{x)A^^\x)Y^{x) 

Wn{x)=Yn{x)Wi°\x)Y^{x) 



(3.42a) 
(3.42b) 
(3.42c) 
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(The last rule actually follows from the second.) From the property 

[{n-d-ign-A^,)Ynix)]=0, (3.43) 

we find that upon making the replacements (|3.42|) in the SCETi Lagrangian ()3.38|1 . the 
term containing the ultrasoft gluon disappears! That is, 

Ux)in.D^Ux) ^ ^^^\x)in . {d - igA(^^)^e^\x). (3.44) 

Thus, in terms of the fields , An^ , there are no couplings of collinear particles to ultrasoft 
gluons at leading order in SCETj. The same decoupling occurs in the gluon Lagrangian 
(|3.39|) . Since nonperturbative physics is governed by interactions with ultrasoft gluons, this 
decoupling greatly simplifies the separation of perturbative and nonperturbative contribu- 
tions to physical observables in the effective theory. 

Ultrasoft gluons have not entirely disappeared from the theory, however. When we 
match currents or operators mediating various decays from QCD onto SCETi, we must 
include ultrasoft Wilson lines Yn{x) whenever we have a collinear field in the operator, 
according to Eqs. (|3.42j) . (Equivalently, we must ensure that all operators are invariant un- 
der ultrasoft gauge transformations.) These rules will become apparent in the applications 
presented in subsequent chapters. 



3.3.3 Effective Theory for Exclusive Decays — SCETn 

We now have all the tools in the effective theory required to describe inclusive decays such 
as in Z decays to hadrons in Chap. El Before proceeding to this application, however, let 
us introduce the novel features of a second version of SCET — SCETn [HH] — required to 
analyze exclusive decays as in Chap. 21 on radiative exclusive T decays. 

To see the inadequacy of SCETi in analyzing an exclusive decay producing energetic 
light hadrons, such as T — > 'ymr, consider the typical invariant mass of these light hadrons. 
For a particular light hadron, such as the pion, the invariant mass squared is of the order 
~ Aqqj-). The typical momentum of a collinear particle in SCETj, however, scales as: 



(3.45) 
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where A = y^Aqcd/Q- The invariant mass squared of such a particle is: 



P^=P+p--plr.Q^X^r.QA^CB, 



(3.46) 



which is too large to represent a light particle with ~ Aq 



QCD- 



The correct effective theory is formed by using instead the parameter r] = Aqcd / Q and 
collinear fields which scale as: 



giving an invariant mass squared ofp'^ ^ Q'^'if = Aqqj-,, which is now correct for an exclusive 
light hadron. These collinear fields interact with soft fields, whose momenta scale as: 



with invariant mass squared pi ~ Aqqj-, . The effective theory for collinear and soft degrees 
of freedom with these scalings is called SCETn- 

The soft fields of SCETn are in fact the same as the ultrasoft fields in SCETj. The 
collinear fields, however, live at a lower energy scale in SCETn than they do in SCETj. 
For this reason, the collinear fields in SCETi are often called hard- collinear fields to dis- 
tinguish them from the collinear fields of SCETn- The introduction of another degree of 
freedom, the soft- collinear messenger mode, was introduced in Refs. |941 195j so that all 
infrared divergences in one-loop graphs in SCETn and QCD would match, which can also 
be accomplished by introducing a proper infrared regulator in the effective theory, as in 
Ref. 96 . We do not evaluate any loop graphs with these divergences in the subsequent 
chapters, and so blissfully ignore these modes or regulators. 

To work in SCETn, one can imagine matching directly from QCD onto SCETn- How- 
ever, it is more convenient to match QCD first onto SCETi, and then match SCETj onto 
SCETn, since we have already done the hard work of completing the first step- In the sec- 
ond step, we take the collinear fields of SCETj, and lower their off-shellness from the scale 
\/QAqcd to Aqcd; putting them in SCETn- Interactions with ultrasoft gluons are already 
removed from the SCETi collinear quark and gluon Lagrangians via the field redefinition 
(|3.42|) . We simply change ultrasoft Wilson lines Yn{x) in SCETi into soft Wilson lines Sn{x) 



{p'^,p ,p_l) ~ Q(7?2,l,r/), 



(3.47) 



(3.48) 
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in SCETn which look the same as in SCETj but contain the soft gluon fields of SCETn 
with the scaling (|3.48jl . 

Note that in SCETn, there could not be interactions of collinear quarks with soft gluons 
anyway, because the soft gluon knocks the collinear quark off shell, due to the scalings in 
and imHl) . 

We will make use of SCETn in Chap. \^ where the procedures described above will be 
illustrated by explicit example. 

3.3.4 Reparametrization Invariance 

The introduction of the light-cone vectors n, n along which to decompose felds and momenta 
in the effective theory breaks the Lorentz invariance of the original theory of QCD. The 
Lorentz symmetry manifests itself in the effective theory by means of invariance under 
redefinitions of the vectors n, n and of the label momenta assigned to collinear particles. 
These features are known as reparametrization invariance (RPI) |97j . 

We focus here on the theory SCETj for simplicity. The first invariance requirement 
arises from the fact that the decomposition of collinear momenta into label and residual 
components is not unique. We wrote in Eq. ()3.19|) for a collinear momentum p, 

p = p + k, (3.49) 

where p is contains the large 0{Q) and 0{QX) label momenta, and k is the 0{QX'^) residual 
momentum. Shifting an amount of momentum of order QX^ between the label and residual 
momentum should yield an equivalent description of the collinear physics. The main conse- 
quence of this invariance, together with gauge invariance, is that operators in the effective 
theory containing derivatives must appear in the combination: 

Vf^ = Di! + D^„ (3.50) 

where Dc = — igAn^q, and Dug = — igAus- This condition relates operators at 
different orders in A, since the terms contained in Eq. ()3.5U|) have different scalings in A. In 
this thesis we will focus on the constructions of operators only at leading order in A, so this 
constraint will not be crucial. 
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More relevant constraints arise from the second type of invariance, that of redefinitions 
of the hght-cone vectors n,n themselves. These vectors enter the very definition of the 
SCET Lagrangian, and the theory must remain invariant for equivalent choices of n, n, 
which are those that leave invariant the conditions: 



= = 0, n-n = 2, 



(3.51) 



must yield an equivalent theory. There are three classes of transformations we can make. 
In the first two, we fiddle just with the transverse components of one or the other of n, ft: 



Type I: { 



n 



Type II: < 



\^ri^ rill 



(3.52) 



and in the last class we mutually rescale the vectors: 



Type III: < 



^ (1 + a)n^ 
^ (1 - a)n^ 



(3.53) 



The parameters A"*-, a are infinitesimal, but can be assigned a particular scaling in 
A. Starting with type-II transformations, consider the change induced in the light-cone 

components of a vector F^. We start with the decomposition: 



2 2 



(3.54) 



Under the transformation, the components become: 



n-V ^ n-V 



n-V ^ {n + £±)-V = n-V + £±-Vi_ 



77^ 

yi" = V'' -n-V— 



n^ 
■ n-V — 
2 



V^ -{n-rei_)-V^-n-V^ ^ ^ 



= y'^-e^-V^\-n-v'-^ 



(3.55a) 
(3.55b) 

(3.55c) 
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If V is, say, a collinear momentum or field with scaling {V'^ ,V~ ,V±) ~ Q{X'^,1,X), then 
we see from the shifts in h-V and V± that we can assign a scaling of 1 to the parameter 
e± without messing up the power counting of p. Similar examination of type-I and type-III 
transformation reveals that cannot have a scaling larger then A, while a can also be 
order 1. Therefore, type-II and type-III RPI can constrain operators at the same order in 
A, while type-I will relate operators at different orders in A. We will make use of type-II 
and type-III RPI to restrict the number of operators which can appear in the description 
of T decays in Chap.^ 



3.4 Non-Relativistic QCD 

We have described effective theories in QCD for hadrons containing a single heavy quark 
(HQET) and for light, energetic hadrons (SCET). We make use in this thesis of one more 
effective theory, that for mesons containing a heavy quark-antiquark pair, QQ. In such a 
meson, both heavy partons move with relatively small velocity. Thus, we are led to expand 
QCD about its non-relativistic limit — hence, the effective theory of non-relativistic QCD 

(NRQCD) piFMirnnj. 

In NRQCD, we separate the fields in the QCD Lagrangian which create quarks and 
antiquarks, into the two-component spinor fields ip and respectively. The small expansion 
in NRQCD is the heavy quark velocity v. This parameter now determines two separate 
physical scales, not just one: the heavy quark three-momentum is of the order |p| ~ itlqv, 
while the kinetic energy is of order E ~ rnqv^. 

To form the NRQCD Lagrangian for heavy quarks, begin again in QCD: 

Cq = ^Q{i]p-mQ)^Q 
The Dirac field \I'q contains fields creating/annihilating both quarks and antiquarks: 




(3.57) 



(3.56) 
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in terms of which the Lagrangian H3.56() is: 



Lq =Tjj\iDo - mQ)^p + xHiDo + mQ)x 



(3.58) 



having used the conventions for Dirac matrices: 

7° = I, f = I " I- (3.59) 



1 o) 




( a'\ 




|, 7-1 




[o -1) 







We can decouple the ^ and x fields by the following transformation: 



which turns (|3.56j) into: 



^«^(# ..ll-^^'""^^ » * , (3.01) 

up to terms suppressed by higher powers of v. Removing the dependence on large label 
momenta p from the full-theory fields: 

^(x) = e-^(™«*-P-^)V'p(a;), x(.x) = e^('"«*+P-^Vp(2;), (3.62) 

and keeping the dominant terms in v, we obtain the leading-order Lagrangian for quarks 
and antiquarks in NRQCD: 

/:nrqcd = ^U^) - ^) M^) + xiix) (iD^ + ^) Xp{x). (3.63) 

The fields ipp, Xp describe fiuctuations about the label momentum p of order rnqv^: 

P = p + k, E = k^, (3.64) 

where P is the total heavy quark or antiquark three-momentum, and the residual momentum 
k = (A:0,k) ~ Mt>2. 
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The construction of operators to describe quarkonium decays in this effective theory is 
performed in Chap.^ 
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Chapter 4 

Exclusive Radiative Decays of T 

"For I am Saruman the Wise, Saruman Ring-maker, Saruman of Many Colours!" 

I looked then and saw that his robes, which had seemed white, were not so, but 
were woven of all colours, and if he moved they shimmered and changed hue so that 
the eye was bewildered. 

"I liked white better," I said. 

Gandalf, in The Fellowship of the Ring by J.R.R. Tolkien 

We now turn to the application of SCET to exclusive decays — the radiative decays of quarko- 
nia to an exclusive light hadron in the final state. This chapter is based on work published 
in Ref. 02]. 

4.1 Introduction 

In a recent series of papers the differential decay rate for the decay T — > has been studied 
in the "endpoint" region where the decay products have a large total energy of order the 
T mass (Mx), and a small total invariant mass squared of order AMx, where A ~ 1 GeV 
is the typical hadronic scale j4Ul 1411 1391 1371 llUlj . An important tool in this analysis is the 
soft-collinear effective theory (SCET) [891 19U1 IHTl I92j . which is a systematic treatment of 
the high energy limit of QCD in the framework of effective field theory. Specifically SCET 
is used to describe the highly energetic decay products in the endpoint region. The heavy h 
and b quarks which form the T are described by non-relativistic QCD (NRQCD) [981 llUUj . 

The soft-collinear effective theory is not limited to applications involving inclusive pro- 
cesses. In fact SCET has been extensively applied to exclusive decays of B mesons into 
light mesons ^ IMl EHl CUSl IMl HEZl HDHl Unni ma Here we use similar 

techniques to study exclusive radiative decays of the T. We make use of some of the results 
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derived in the analysis of inclusive radiative decays in the endpoint region, but the analysis 
of exclusive decays is complicated by the existence of two different coUinear scales. This 
necessitates a two-step matching procedure jHS]- In the first step one matches onto SCETi 
which describes collinear degrees of freedom with typical offshellness of order -y/AMx, as is 
appropriate for inclusive decays in the endpoint region as discussed above. In the second 
step SCETi is matched onto SCETn, which is appropriate for exclusive processes since it 
describes collinear degrees of freedom with typical offshellness of order A. 

The analysis of T decay is further complicated by the existence of two types of currents: 
those where the bb is in a color-singlet configuration and those where it is in a color-octet 
configuration. The octet operators are higher-order in the combined NRQCD and SCET 
power counting, so one might suppose that they can be dropped. However, the octet 
currents have a Wilson coefficient which is order as(Mx) while the singlet current has 
a Wilson coefficient of order as(Mx). The additional suppression of the singlet Wilson 
coefficient is enough so that both color-octet and color-singlet operators must be included 
as contributions to the inclusive radiative decay rate in the endpoint region |41l I39j . 

In this work we show that in exclusive decays the octet currents are truly suppressed 
relative to the singlet current and can be neglected. We then determine the minimal set of 
color-singlet currents which can arise and fix their matching coefficients in SCETj. We run 
this current to the intermediate collinear scale ~ \/AMj- and match onto SCETn. Our 
expression for the decay rate agrees with that derived in QCD using a twist expansion jll21 
11131 1114j . Finally we use our results to make a prediction for the ratio of branching fractions 
B{T 'yf2)/B{J/ip — > 7/2), B{J/ip ^f2)/B{ip' — > 7/2), and analyze the decay T — > ^-mr 
in the kinematic regime where the pions are collinear. 

4.2 Power Counting 
4.2.1 Inclusive Decays 

The first step is to match the QCD amplitude for a bb pair in a given color and spin 
configuration to decay to a photon and light particles onto combined SCETi and NRQCD 
currents. The SCET power-counting is in the parameter A ~ y^A/M, where M = 2mb, 
while the NRQCD power-counting is in v, the relative velocity of the bb pair in the T. 
Numerically, A ~ f ~ 1/3. The matching is shown graphically in Fig. 14.11 The effective 
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Figure 4.1: Matching onto operators in the effective field theory with one and two gluons 
in the final state. The currents on the left have a color-octet bb in either a ^Sq or "^Pj 
configuration. The matching for a color-singlet bb pair in a '^5i configuration is shown on 
the right. 

theory operators can be classified into those with the bb in a color-octet configuration (shown 
on the left-hand side of Fig. 14.1(1 and those with the bb in a color-singlet configuration (shown 
on the right-hand side of Fig. 14. 1|) . The leading octet operators can be further subdivided 
into a those where the bb is in a ^5*0 configuration and those where the bb is in a ^Pj 
configuration. The octet operators are [IH IHH] 

J^{8,'So)=Y,cf''{M,f^)Ti^xlpBl4^p, (4.1) 

i 

where 

= -^T^tfpM + 5,(AM (4.2) 

Ss 

The operator "P^ projects out the label momenta in the perpendicular direction The 
sum in Eq. (|4.1|) is over all possible Lorentz structures denoted by F^^ , and C- ' " (M, fi) is 
the corresponding matching coefficient for each structure. The octet ^Pj operators are 

J^(8, ^Pj) = J2 cf'iM, f,)r^^^sxlpBlA . p-A • T-^Vp , (4.3) 

i 

where A is a Lorentz boost matrix. Each of these color-octet operators scales as 0(A) in 
SCET. The NRQCD power-counting has the ^Sq octet operators scaling as 0{v^); however, 
this operator has an overlap with the T state beginning at 0{v'^). (See Appendix lOTi 
Thus the ^Sq operator contributes at order v^X to the T radiative decay rate. The ^Pj 
octet operator has NRQCD scaling 0{v^), but overlaps with the T at order v. Thus the 
total power-counting of the ^Pj contribution is 0{v^X), which is the same as the ^5*0 octet 
operators. The leading order matching coefficients for both are 0{^/aJ(M)). 
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The color-singlet operators are 



J^(l, 'S,) = Y, K/35,xlpA . *T^Vp1V{B^ (M, P+) 5? } , (4.4) 

i 

where V-^ = + V, with V = n ■ V . These operators scale as C'(A^) in SCET and 0{v^) 
in NRQCD. The leading matching coefficients are 0{as{M)). Thus the ratio of color-octet 
to color-singlet contributions in inclusive radiative T decay scales as: 

octet V , , , 

(4.5) 



singlet X^as{M) ^/ajM) ' 
4.2.2 Exclusive Decays 

The situation changes when one considers exclusive decays. The currents we just discussed 
are SCETj currents where the typical invariant mass of the collinear degrees of freedom is of 
order /ic = \/ MA. In order to have overlap with the meson state we must match onto SCETn 
currents where the typical invariant mass of collinear particles is 0{K). Furthermore, the 
interpolating field which annihilates the meson state in SCETn is defined to consist only of 
collinear fields in a color-singlet configuration [HSI- Given these considerations it is simple 
to match the color-singlet operator in SCETj to an operator of identical form in SCETn- 
However, the matching of the octet contributions from SCETi onto SCETn is more involved. 

Before we consider the matching of the octet contributions from SCETj onto SCETn 
we turn our attention to the scaling of these contributions in SCETj. In order to produce 
a final state consisting only of collinear fields in a color-singlet configuration we need an 
interaction which changes the ultrasoft (usoft) gluon into a collinear gluon (as shown on the 
left-hand side of Fig. I4.2|) . This term in the SCET Lagrangian is power suppressed by A so 
that the time-ordered product of the octet current with the collinear-collinear-usoft vertex 
scales as 0{\^) in the SCETj power counting. In addition, the exchanged gluon introduces 
an extra factor of the coupling constant at the matching scale: Qs(/ic)- Including these 
factors, the time-ordered product of octet currents with the subleading Lagrangian scales as 



0{as{^J.c)^ Cis{M)}?v^), and the ratio of time-ordered products to the singlet contribution 
is 

octet ^v^^osi^^^^^^^ ^^g^ 



singlet ^Jas{M) 
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for the bottomonium system. For charmonium the above ratio is about 0.2. This result is 
very different from the result for the inclusive decay given in Eq. I|4.5jl . In SCETj the octet 
contribution to exclusive radiative T decay is not only suppressed in the limit v, X ^ 0, but 
numerically suppressed by a factor of ~ 10 for typical values of the parameters. This is the 
same order of suppression we expect from higher order SCET and NRQCD corrections; thus, 
we should be able to safely neglect the color-octet contribution in SCETj. However, before 
we can neglect the octet contribution in our analysis we must show that the suppression of 
the octet piece holds after matching onto SCETji. 

We first turn our attention to the simpler calculation: matching the color-singlet oper- 
ator. In SCETi we perform the field redefinition |92j : 



which decouples usoft from collinear degrees of freedom. Under this field redefinition 



and the Y^s cancel in the trace of the color singlet operator given in Eq. ()4.4() . Thus we 
match the SCETi operator after the field redefinition onto an operator in SCETn of a form 
identical to that in Eq. 1)4. 4() : 



jy^)(l,35i) = j;rj,^,^xlpA-TVpTV{i?fi^cf' ''\M,P+-t,,)B^,J , (4.9) 



where fic = V MA is the SCETi-SCETn matching scale, and the subscripts indicate SCETn 
fields. From now on we drop the subscripts. In SCETn the power-counting parameter is 
r/ ~ A^, and the SCETn color-singlet operator in Eq. (|4.9j) is 0{v^r]'^). The short-distance 
coefficient is inherited from SCETj and is 0{as{M)). 

The matching of the color-octet current is more complicated. In order to match onto 
an SCETn operator with color-singlet collinear degrees of freedom we must consider time- 
ordered products where a usoft gluon radiated from the bb pair is turned into a final state 
collinear degree of freedom. An example of such a diagram is given in Fig. 14.21 Two 
collinear gluons are required for the collinear final state to be color-singlet. One of the 
collinear gluons comes from the octet current, and the other can be produced by pulling 



(4.7) 



(4.8) 
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Figure 4.2: Matching from SCETj to SCETn- This Feynman diagram is an example of 
a time-ordered product in SCETi that matches onto an operator in SCETn that has a 
nonzero overlap with the final-state collinear meson. 

a gluon out of the bbg Fock state of the T, and kicking it with a collinear gluon from the 
current. This requires a collinear-collinear-ultrasoft coupling which first appears at order A 
in the SCETj Lagrangian [TT^ ITTHl ITTTj : 

4^) = ^Tr [\iV'^,iDi^]\iV^,iDi]] , (4.10) 

where = Dc + n- Dusn^ /'^- The decay amplitude comes from a time-ordered product 
of the color-octet current and C^cg , or a time-ordered product of the color-octet current, 
C^cg , and a leading order gluon interaction. Though our result will hold for either type of 
time-ordered product we will, for the sake of concreteness, only consider the former: 

Ts = j d^xr{j(8,-)(0),4?(^)}, (4.11) 

where the dot stands for ^Sq or '^Pj. In the time-ordered product two gluon fields are 
contracted to form the internal propagator in Fig. 14. 2( which scales as 1/}? . We require 
two uncontracted fields (in a color-singlet configuration) so that we can match onto an 
SCETn operator in the form Tr[A^74^] which annihilates the final state collinear meson. 
In this example the leading contribution is an n-An gluon field in one of the Wilson lines in 
J(8, •) contracted with an n-An field in C^cg ■ After the contraction what remains in C^cg is 

^TT{[gT^,Ain[V,Ai]], (4.12) 

which scales as (A^)(1)(A) = A^. Note we now have the correct field content for the operator 
shown on the right-hand side of Fig. 14.21 there are two outgoing A^ fields, one from C^cg 
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and one from J(8, •), in a color-singlet configuration, and an incoming soft gluon field also 
from C^cg ■ 

Next we decouple collinear and usoft in SCETj through the field redefinition in Eq. I|4.7|) . 
This introduces factors of y and into our expressions. When matching onto SCETn these 
become soft Wilson lines S and S"^. Since these Wilson lines do not affect the power counting 
we ignore them. Now we can match onto a convolution of SCETn operators with SCETj- 
SCETn matching coefficients. Since these arise from integrating out the internal collinear 
propagators they scale as A^'^ for each propagator. In our example there is one propagator 
so the matching coefficient scales as 0{\~'^), which is 0{'q^^) in the SCETn power counting 
(remember r/ ~ A^). Since the SCETn operator has two fields each scaling as rj, and a 
soft field scaling as r], it scales as . Combining the scaling of the SCETn operator with 
the scaling of the SCETj-SCETn matching coefficient gives an 0{as{lJ-c)v'^) contribution. 
If we include the order NRQCD scaling from the heavy sector, and the 0{^JciJJ\d)) 
contribution from the QCD-SCETj matching coefficient the color-octet contribution to 
exclusive decays scales as 0{v^rj^ \/ as{M)as{lJic)) in SCETn- Taking the ratio of the color- 
octet to the color-singlet contribution to exclusive T decay in SCETn we find: 

octet v'^ttsinc) f. 
— — ; — ~ — , , (4. id) 

smglet ^/as{M) 

which is the same scaling we found in SCETj. Thus we can safely neglect the color-octet 
contributions at this order. 



4.3 Complete Basis of Color-Singlet Matching Coefficients 

Now that the color-octet contribution has been eliminated we determine a complete basis 
of Lorentz structures E^^^^ that can appear in the color-singlet matching coefficient in 
Eq. (|4.4() . At leading order in as{M) only one Lorentz structure was found to be non- 
zero j41j : 

C{ 'iM,uj)T^pSf, = ^j^9af39t.5- (4.14) 

However, at higher order other Lorentz structures may appear. These coefficients can be 
constructed from the set: 

{giiu,nf„nf„Vf,} , (4.15) 
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where v is the four- velocity of the T, under the restriction that T^^^^ satisfies the appro- 
priate symmetries. For example, the full theory ampHtude is parity even, as is the effective 
theory operator, meaning that the matching coefficient must also be parity even. As a result 
the epsilon tensor is not included in Eq. 1)4. 15() . 

We treat v as an object independent of n, n |117j . and use = = 0, n ■ n = 2, and 
= 1. Before we write down all the possible operators which can appear we note some 

simple properties that will make our task more manageable. First we note that T^fBS^ must 

be symmetric in a and /?. To see this consider the object: 

cI'''''^\m, co)r^ps,T^r {Bl6^,pyj , (4.16) 

which is the collinear part of the color-singlet operator where a sum over lo has been in- 
troduced. First interchange a and f3, and then use the cyclic nature of the trace to switch 
the two B± fields. Note, however there is a projection on these fields from the operator 
in the Kronecker delta involving V+. Since V+ = +V this operator projects out mi- 
nus the label on i?" and projects out the label on To preserve this relationship 
when the order of the fields is switched we must let 6^^p^ ^u),-v+- letting uj — > —uj 
we have = ^u,'P+-> operator goes into itself. However, the Wilson coeffi- 
cient is now Cj-"*^' ^^\m,—ijj). To demonstrate that Eq. 1)4. 16() is symmetric under a /? 
we must show that c\^' ^^\m,u)) is even in uj. We use charge conjugation for this. The 
heavy quark sector of the operator has charge conjugation C = — 1 as does the photon. 
As noted in Ref. [TTS^ two gluons in a color-singlet configuration must have C even. Since 
QCD is charge conjugation conserving the product of operator and coefficient in Eq. 1)4.16(1 
must also be C even. This is the case if the matching coefficient C^^'' ^^\M,uj)Tl^i^g^ is C 
even. Following Ref. |119j . under charge conjugation the above product of operator and 
coefficient goes to itself with oj — > —lo in the coefficient function. Thus charge conjuga- 
tion implies C^^' ^^\m, —uj) = C^^' '^^\m,uj), and as a result Eq. (|4.16|) is symmetric in 
a and (3. Second, any vs appearing in T^f^Sfi gives zero contribution to the operator, since 
V ■ A = 0. Third, na,na (and by symmetry nfj^np) appearing in the operator also gives 
a zero contribution since these indices contract with indices on the B± field. Finally, we 
use reparameterization invariance (RPI) of SCET [971 I115| . The terms satisfying these 



requirements are 
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CI '{M, Uj)T^pSii = CigalBgSfi + C25a/37 ^ + CsQaf^ (4.17) 



+ C4 



ga^l [gf3s 1 + [gas [gpi, 

n ■ V / \ n ■ V / V n ■ V / \ n ■ v 



So far we have allowed v to be an arbitrary vector. Now we restrict ourselves to a 
frame where Vj_ = 0. Furthermore we are interested in the case where the photon is real, so 
we can restrict the photon to have transverse polarizations. This leaves only two linearly 
independent terms: 

cf ' '^^\m, uJ)^„ps^^ = algif^gi^ + al {g^sgj^ + g^f^gjs - gapd,^ , (4.18) 

i 

where g'^ = g^" — (n^n'^ + n^n'^)/2. The first term projects out the trace part of the 
Tr v+'^a) operator, while the second term projects out the symmetric traceless com- 

ponent. Since the Lorentz symmetry in the perpendicular components is not broken in 
SCET these two terms do not mix under renormalization. The leading-order matching fixes 
the coefficients ai and 02 to order Os- 

a?(P+;^ = M) = ^|^, af(P+;/. = M) = 0. (4.19) 

Since there is no mixing, a2 = + 0{a1{M)) at all scales. Note this matching assumes the 
T states are non- relativistically normalized: {T{P')\T{P)) = 5^{P - P'). 

In addition to gluon operators we must consider the basis of all possible quark operators 
which can appear in radiative T decays 

Jl{l, 35i) = xLpA • Vp4,p, W^r;5(M, V+)W^in,p, . (4.20) 

i 

The basis of Dirac structures, {^,^75,^7^}, was given in Ref. jll9j . and the most general 
basis of quark operators can then be constructed out of these Dirac structures and the set 
{^af3fj.u, g^iu,n^,n^,v^}. Using the symmetries of SCET and RPI we find 

E^s = al^g^s + 4^1^4s + 4^1^ns . (4.21) 
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The first term transforms as a scalar, the second term transforms as a pseudoscalar, and the 
third as a vector. The matching coefficients at the scale ^ = M for the quark operators in 
radiative T decay are all zero at leading order in perturbation theory, but the scalar quark 
operator mixes with the scalar gluon operator through renormalization group running and 
can be generated in this manner. The pseudoscalar and vector term do not mix with the 
scalar gluon operator due to Lorentz symmetry and will not be generated at this order in 
the perturbative matching. 



4.4 Decay Rates &6 Phenomenology 

We now consider the phenomenological implications of our analysis for exclusive radiative 
decays of quarkonium into either a single meson or a pair of mesons which are collinear. 
The (n + l)-body decay rate is given by: 

- ^^") - 2^ / n ^^^^^^ 

X K7(g)pi...Pn|JM^|T(P))|2, 

where J is the QCD current, A is the photon field, and Fn denotes an exclusive final 
state consisting of n collinear particles. We consider only decay rates where the final state 
momenta pi are all collinear with combined invariant mass m,^ = (Yl^ Pi)"^ ~ ^- The 
effective theory decay rate is obtained by matching the current J onto the SCETn current 
given in Eq. (|4.9|1 plus a quark operator: 

(7(g)pi...Pn|JM^|T(P)) 

i 

i 

= (OlxLpA.^Vp |T(P)) {j{q)\A^^ \0)gi 



(pi...p„,|Tr{P2a?(P+;/.)P^}|0) + {pi...Pn\Cn,p,W^ "'^m'^^ ^W|0) 



(4.23) 
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In obtaining the second line we make use of the results in Eqs. H4.18() and ()4.19|1 . and use 
the properties of SCETjj to factor soft and colhnear degrees of freedom. In the last line we 
changed to a nonrelativistic normalization for the T state. 

Next we define the light-cone wave functions 

{pi . ..pn\PTv [B^S^^r^B^] |0) = M3-"0f"(x), (4.24a) 
{pi ■ ..Pn\U^,W^6^^p^W^(n,., |0) = Af3-«<^^"(x) , (4.24b) 

where states are relativistically normalized, and the discrete label uj is converted to a contin- 
uous one, X = oj /n ■ p, as explained in Ref. |inij . The wave functions are dimensionless. 
See Appendix 14. Bl for the relation of these SCET light-cone wave functions to those con- 
ventionally defined in QCD. Then the collinear matrix elements in brackets in Eq. 1)4. 23|) 
can be written as the convolution: 

{pi...Pn\Tr {Bl al{P+; fi) i3^}|0) + {pi...Pn\^n,p,W^ al{P+; fx) W^^n,p,\0) 

l-i (4.25) 
= M^-" j ^ dx (af (x; ^)<Af (x; fi) + a?(x; /z)<Af" (x; /x)) . 

The dependence on the scale // cancels between the long-distance matching coefficients and 
the wave function. We will elaborate on this point in a moment. First we expand both 
; /i) and 0^"^(x; /i) in Gegenbauer polynomials: 

a?(x;/.)= ^ aipHf,)Cl/\x), (4.26a) 

n odd 

?(x;/.) = a(")(/.)(l - x')ct\{oc) , (4.26b) 



n odd 



f (x;/x)= b^^\fi){l - x')C'J\x) , (4.26c) 

n odd 

f (x;/.) = Y - ^')Cn~i(x) ■ (4-26d) 



n odd 



Then the convolution becomes an infinite sum of products of Gegenbauer coefficients: 

/ dx{a'({x;fj.)(pg"{x;iJ.) + al{x; fj,)(j)g" {x; n)) 

(4.27) 

n odd 
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where 



(„) ^ n(n+l)(n + 2)(n + 3) („) ^ (n + l)(n + 2) 

^5/2 9(n + 3/2) ' ^3/2 ^ + 3/2 • ^^"^^^ 



We now return to the question of the scale. Here we pick /i ~ A which minimizes 
logarithms in the wave function; however, large logarithms of M/A then appear in the 
Wilson coefficients. These large logarithms are summed using the renormalization group 
equations in SCET. This calculation was carried out in Ref. jlUlj . and we only quote the 
results here. We find: 



dx{al{x; "(x; fx) + a?(x; ^i)(j)g"{x; ^)) 



-a?(M) ^ 



n odd 



\as{M)J \as{M)J 



+ 



fi"^) (n) 

■^3/2 Tag 

An) ^(n) 
■'5/2 



2A("V/3o' 



as{M) 



2A("'//3o 



as{M) 



2aL"V/3o' 



(4.29) 



where 



/3o = 11 - ^ 



(n 



(n 
7ijg 



7^" 



7^" 

'93 



'95 



2^/ 
3 



(n) 
793 



A 



(n) 



1 r 
2 



7 



Iqq 



(n) 



7ST + 47j?7£^ 



1 



(n+l)(n + 2) 2 

Cf + 3n + 4 
"3"(n + l)(n + 2) ' 

3n + 3n + 4 

"'■^n(n + l)(n + 2)(n + 3) 



1 "+^1 

--2y- 



1=2 



+ 



n(n + l) (n + 2)(n + 3) 6 



1 

i-2yi 



3 ' 



(4.30a) 

(4.30b) 
(4.30c) 

(4.30d) 
(4.30e) 

(4.30f) 
(4.30g) 

(4.30h) 



The quantities a|^^ which appear in the exponents in Eq. (|4.29|) are negative for any n > 1. 
Furthermore A^'' < 0, while X^P = 0. This property allows us to consider the asymptotic 
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limit M A, where as{M) 0. Then 



1 



lim 

M-^QO 



-1 



(4.31) 



16 Cf 



i{M) L(i)(A) + ^6«(A)] ^B^-a{{M) 



which defines a nonperturbative parameter B^" . However, for values of M around the T 
mass this is not a very good approximation, and for values around the J /tp mass a much 
better approximation is to assume no running at all. 

4.4.1 Two-Body Decay: T 7/2 

Having taken care of the technical details we can now use the above results to study the 
two-body radiative decay T — > 7F1 . The decay rate is 



where the full expression for the term in brackets is given in Eq. 1)4. 29(1 . After factoring, 
the soft matrix element involving the heavy quark fields was further simplified using the 
vacuum insertion approximation for the quarkonium sector, which holds up to corrections of 
order PSI. Note that the operator above overlaps only with the A = ±1 helicities of the 
T. Then using the rotation symmetries of NRQCD |12Uj we can relate the non-relativistic 
matrix element above to those conventionally used: 



For the final state meson Fi to have nonzero overlap with the operators in Eq. 1)4. 24() it 
must be fiavor singlet, parity even and charge conjugation even. One candidate with the 
correct quantum numbers is the /2(1270). Furthermore this decay has been measured both 
in T and J/V' radiative decay, which is why we consider it. An interesting point is that only 
the helicity A = component of the /2 contributes at the order to which we are working. 
To see this begin by considering the decomposition of the following gluon matrix element 



r(T ^ 7i^i)scET„ = r^ml:<y\x~p'xl^ol^v \^) 




(4.32) 



^IX-p'xLp^lV'p |T) = ^(T|0(l,35i) |T) . 



(4.33) 
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into all possible light-cone form- factors: 



(/2|Tr [BIB^J |0) = A{e{X))gf + B^ef{\), (4.34) 

where e"^ is the symmetric-traceless polarization tensor of the f2- We give the explicit 
form in Appendix I4.CI There are only two form factors above since the matrix element 
must be decomposed into tensors that have non-zero perpendicular components. The only 
structures available are g'^ and . For A = ±1, e^{\ = ±1) = 0, so this helicity 
component does not appear at this order. The coefficient A{e{\)) is a scalar function 
which can be constructed from Tr (ej,) and na,npe°^^ . Because the helicity-zero polarization 
tensor has the property that e^^ [X = 0) cx ^ and the helicity- two polarization tensor has 
Tr(e_L) = and fianpe'^^ = 0, we can fix the normalization of the coefficient A{e{\)) so 
that the first term on the right-hand side of Eq. ()4.34|) parameterizes the A = contribution 
while the second term parameterizes the A = ±2 contributions. The helicity-zero piece is 
picked out by the aig'^ term in Eq. 1)4. 18() . while the helicity- two piece is picked out by the 
02 term. Thus at leading order in perturbation theory, the dominant decay should be to 
the helicity-zero component of the /2. 

The NRQCD matrix element in Eq. 1)4.321) can be expressed in terms of the leptonic 
decay width of the T. At leading order, 

r(T - e+e-) = ^^(T|0(l,35i) |T) , (4.35) 

and the decay rate for T — > 7/2 can be expressed as: 

r(T ^ 7/2) = 16^«jW (j3/.)2r(T ^ e+e-) . (4.36) 

We can repeat the same analysis for the decay rate T{J/ip — > 7/2) and form a ratio of 
branching fractions, which in the asymptotic limit is: 



(4.37) 



where Mqq = 2mQ. Using rub = 4.1 - 4.4 GeV, = 1.15 - 1.35 GeV, B{T e+e") = 
(2.38 ± 0.11) X 10-2, and B{J/ip e+e") = (5.93 ± 0.10) x 10"^ p, we predict the ratio 



B{T - 


> 7/2) 




^ /4CF + 3^ 


.2 B{T^e+e-) 


B{J/^ - 


7/2) 




\4CF + ij 


1 B{J/^^e+e-) 
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of branching fractions to be in the range: 



^(T ^ 7/2) 
B{J/^ ^ 7/2) 



0.14-0.19. 



(4.38) 



As was mentioned earher the asymptotic Hmit is not a particularly good approximation for 
the T, and quite bad for the J/ip- We can improve this approximation by keeping more 
terms in the resummed formula in Eq. (|4.29l) . The dominant term is the part of the n = 1 
term proportional to b^i^\fi), and in this approximation: 



BjT ^ 7/2) 
B{J/ij ^ 7/2) 



as{Mi 



hh) 



7i"-7i"(a,M/a.(A%))"' * 



L7i" 

e^e 



(1) , '"'f~ 



B{J/^ 
0.13-0.18, 



(4.39) 



where n ^ 1 GeV. The range of values has not changed much from Eq. 1)4. 38|) . however, 
theoretical errors are reduced: corrections to Eq. ()4.39|) from the b'^^ and higher-order 
terms in Eq. (|4.29|) are estimated to be roughly 40%, while corrections to the infinite mass 
limit from higher order terms are estimated to be roughly 80%. In addition there are 
theory errors from neglecting higher-order terms in the perturbative expansion, as well as 
in the expansions in v and r]. Our prediction can be compared to the measured value of 
0.06 ± 0.03, using the measurements B{T 7/2) = (8 ± 4) x 10"^ and B{J/ip 7/2) = 
(1.38±0.14)xl0-3 p. Given the theoretical errors we can only conclude that our prediction 
does not disagree with data. 

Our predictions for the ratios of T and J/^|J branching fractions to 7/2 are consistent 
with those derived in Refs. |112| I113| I114j . which use an expansion in twist. In particular, 
we reproduce the suppression of the helicities |A| = 1,2 in the final state relative to A = 0. 
In contrast with Ref. jll4j . we extract the NRQCD color-singlet matrix elements from 
the leptonic decay widths of T and J/ip instead of the decay widths to light hadrons, for 
which corrections from color-octet contributions must be taken into account for a reliable 
calculation |121j . The leptonic decay width, however, receives large corrections at NNLO in 
perturbation theory |122|I123] . In either case, one hopes that the uncertainties are mitigated 
in taking the ratios of branching fractions. 
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We can also compare the decay rates of J/V' and iIj' to 7/2 predicted by Eq. (|4.32|) at 
the matching scale /i = M, where a[{x; M) = and a\{x; M) is a constant. Dependence on 
the integral of the wave function 4>g^{x; M) cancels out in the ratio of branching fractions: 

B{ip' 7/2) Bi^l)' e+e ) 

while the measured value is 6.57 it 1.42. We used B{il)' e+e~) = (7.55 it 0.31) x 10~^ and 
B{iP' 7/2) = (2.1 ± 0.4) X 10-^ m. 



4.4.2 Three-Body Decay: T — > 'yirir 



Next we consider a two pion final state in the kinematic region where the pions are collinear 
to each other with large energy and small total invariant mass mTryr A. In this case we we 
have a three-body final state where the two pions are collinear. It is convenient to define 
the variables: 

mL = (Pl+P2)^ ^=^- (4-41) 
In terms of these variables the differential decay rate is 



dx{a{{x; f^)4>^'' {x; /i) + a?(x; nWq^ {x; ^)) 



(4.42) 



to leading order in rn^^/M^. The properties of the meson pair light-cone wave function (j)^'^ 
have been investigated in Refs. jl241ll25] . which interestingly find that in the region where 
A ^ m-TrTr ^ -^T they are given by an integral over two single-particle wave functions. The 
ratio of the T and J/V' rates to 77r7r in the kinematic region of low is numerically the 
ratio in Eq. 1)4. 39() times an extra factor of m^/rn^ ~ 0.07 — 0.1, that is, 0.01 — 0.02. This 
suppression is due to the much larger total phase space available in T ^ ^ttti relative to that 
in T — > 7/2. No T jirir events have yet been observed in the region m^^ < (1.0 GeV)^ 

US 



4.5 Conclusions 
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We have systematically analyzed the exclusive radiative decays of quarkonium to energetic 
light mesons within the framework of soft-collinear effective theory and non-relativistic 
QCD to leading order in the effective theory power counting, as well as to leading order 
in the strong coupling. We show that color-octet contributions are suppressed by a factor 
of v'^as{fic)/ y^ctsiM) ~ 0.05 in exclusive T decays, and can therefore be safely neglected. 
This is different from the situation in inclusive radiative decays in the endpoint region where 
octet contributions must be kept. 

We then turn to the color-singlet contribution. The tree-level matching onto this op- 
erator is carried out in Refs. |411 139j : however, the authors do not consider the complete 
set of operators that could appear in this decay. We use the symmetries of SCET and 
NRQCD, including RPI, to show that the operator which is matched onto in Refs. [411 139j 
is the only operator that can appear for the decays in question. We also consider the set 
of possible quark operators which can arise. Again only one of the possible quark opera- 
tors can contribute to the decays we are interested in. This operator has zero matching 
coefficient, but it can be generated through running. We use the results of Ref. |101j for 
the renormalization group mixing of the quark and gluon operators, thus resumming large 
logarithms. Our results agree with an analysis in twist carried out in Refs. jll21 Ill^Jl I114j . 

Finally we study the phenomenology of quarkonium radiative decay to the /2, as well 
as to vrvr where the pions are collinear. We make predictions for the ratios of branching 
fractions B{T ^f2)/B{J/tp — > 7/2) and B{J/tp — > ^f2)/B{ip' 7/2), as well as for the 
differential decay rates of T and J/^|J to jtttt in the kinematic region of two collinear pions. 
Our predictions for the decays to 7/2 are consistent with experimental data, but with large 
theoretical uncertainties, while there is insufficient data for ^tttt with which to compare. 
Further theoretical work and more experimental data, especially for the light-cone wave 
functions of /2 , will improve the precision of these predictions greatly. 

4. A Appendix: Power Counting of States in NRQCD 

The color-octet operators introduced in Eqs. (|4.H) and (|4.3j) do not overlap with the T 
state until subleading order in v. We choose to represent the T state — a color-singlet spin-1 
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state — as being created from the vacuum with an interpolating field: 



\T{p))r^^PlA.a'x-p\0). (4.43) 

Thus, the NRQCD parts of the color-singlet operator 0{1,'^ Si) in Eq. ()4.4p overlap with 
the T state at leading order in v: 

(0|xlpA-^Vp|T)~O(^;3), (4.44) 

where the is the scaling of the operator itself. However, the color-octet operators do 
not overlap with the T until we put in insertions from the subleading part of the NRQCD 
Lagrangian. For example, for the ^5*0 operator, 

(0|O(8, 'So) |T) ~ {0\T{xlpT^iJp,^i{cT • B)^^} |T) ~ 0{v'') . (4.45) 

The field B scales as while the tp propagator in the time-ordered product scales as 1/f^, 
giving an overall v"^ suppression in the matrix element relative to the color-singlet matrix 
element. Thus the total scaling is 0{v^). The matrix element of the color-octet ^Pj operator 
is the same as the ^5*0: 

(0|O(8,3 Pj) |T) ~ (0|T{xlpT^A-r A-'tVp, V^^^q} |T) ~ 0{v'), (4.46) 

since the operator 0{8,'^ Pj) itself scales as v^, p • D scales as v^, and the ^p propagator 
scales as 1/f^. 

Thus these two color-octet operators contribute at the same order in v to the T-to- 
vacuum matrix element, suppressed by v'^ relative to the contribution of the color-singlet 
operator. 

4.B Appendix: Nonperturbative Matrix Elements and Light- 
Cone Wave Functions 

The matrix elements in Eq. (|4.24j) defining the SCET wave functions ipg"^ can be related to 
conventional QCD wave functions for flavor-singlet mesons. The two-gluon wave functions 
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for a meson with momentum q and net helicity A = 0, ±2 are defined as jl27j : 

{0\Tr G^,{z)Y{z,-z)G,x{-z)\q,X = 0)^^=f^q^qx j' ^dC e''-'i^<t>^siC, f^o) , (4.47a) 



{0\Tr G^,{z)Y{z,-z)G,x(.-z)\q,X = ±2] 



MO 



(4.47b) 



= fs [{^^vp - Quef,p)qa - iq,,e^a - (lue^a,)qp\ J ^ dC e^^''^^^^ (C, //) 
where Y{z, —z) is the path-ordered exponential of gluon fields: 



Y{z, —z) = Pexp 



ig I da ■ A{a) 



(4.48) 



Going to the hght-cone frame where = -^n^ and z'^ = ^n^, we invert these formulas 
to find 



'^5(C;^o) 



'/'s(C;^o) 



-—J— / dz+ e-<«"^+/2^0|Tr G^,{z+)Y{z+, -z+)G\{-z+) \q,X = 0) 



(4.49a) 
(4.49b) 



X (0|Tr G^,{z+)Y{z+, -z+)G^p{-z+) \q, A = ±2) 



where = n-z, and q = n-q. Now we match the QCD fields on the right-hand side to 
fields in SCET: 



n''G^,iz+) 
Y{z+,-z+) 



(4.50a) 
(4.50b) 



where 



with 



(4.51) 



(4.52) 
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Therefore, for example, the matching between the QCD light-cone wave-function (p^ and 
the SCET operator is 

(4.53) 

Integrating over z^ , and converting from the discrete index to to the a continuous iVc where 
C = u)c/q~ we obtain the matching relation between the QCD and SCET light-cone wave 
functions 

•^s'^lC;^) - -f7r(i - C'X'^'-'^^CC;^) • (4.54) 

The SCET wave functions on the right-hand side are given by [cf. Eq. (|4.24l) ]: 

(Ol^Tr [Bl5^,-p^Bj^] \M{q)) = (g-)2<^f , (4.55a) 
el„^(0|PTr [Bld^^p^Bi] \M{q)) = {q-f<P^^^^ . (4.55b) 

Relations between the wave functions for the quark operator in QCD and SCET can be 
derived as in Ref. |119j . 
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4.C Appendix: Spin-2 Polarization Tensors 

The spin-2 polarization tensor for a particle of mass m and momentum k can be built from 
spin-1 polarization vectors using Clebsch-Gordan coefficients to arrive at: 



/o o\ 

1 ±i 
±z -1 

\o oy 

Ik 



1 



e'^'^CA = ±1) = ^— 



2m 



e^'''{\ = 0) = 



/ 
|k| 
±i|k 

V 
/ k^ 






±z|k| 







-— 

2 








\ 

J 



^ 



El J 



(4.56a) 



(4.56b) 



(4.56c) 
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Chapter 5 

Enhanced Nonperturbative Effects 
in Z Decays to Hadrons 

"White!" he sneered. "It serves as a beginning. White cloth may be dyed. The white 
page can be overwritten; and the white hght can be broken." 

"In which case it is no longer white," said I. "And he that breaks a thing to find 
out what it is has left the path of wisdom." 

Gandalf, in The Fellowship of the Ring by J.R.R. Tolkien 

In this chapter, we apply soft-collinear effective theory for inclusive decays — SCETj — to 
the hadronic decays of Z bosons. The material in this chapter is based on Ref. |33j . 

5.1 Introduction 

Some of the most successful applications of perturbative QCD are to processes such as Z 
decay to hadrons or e~^e~ annihilation at large center-of-mass energy, in which a state with 
no strong interactions decays into final hadronic states. Here we will discuss the case of Z 
decay, but the results apply equally well to the other cases. Not only is the total hadronic Z 
decay width calculable, but so are less inclusive infrared-safe quantities like the Z decay rate 
into 2-jet and 3-jet events, the thrust distribution and jet mass distributions. Comparison of 
perturbative predictions for these and other quantities with experimental data on Z decays 
from LEP and SLD has led to a remarkably accurate extraction of the strong coupling 
constant as{Mz) |HE1 lEHl USHl EOl CHIl Ea IISHI- Although the extraction of the strong 
coupling from event shape variables is less accurate than from the total hadronic Z width, 
it is more model-independent since (neglecting quark mass effects) it does not depend on 
the values of the quark couplings to the Z. 
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For the totally inclusive hadronic Z decay width, the operator product expansion allows 
one to include in theoretical predictions nonperturbative strong interaction effects that 
are characterized by vacuum expectation values of local operators. The effects of higher- 
dimension operators are suppressed by powers of the strong interaction scale Aqcd divided 
by the center-of-mass energy Mz- Since the Z mass is large, these effects are very small. For 
example, if quark masses are neglected, the leading nonperturbative effects in the Z decay 
width come from the vacuum expectation value of the gluon field strength tensor squared, 
{G ^yG^'^) . This dimension-four operator gives rise to corrections to the total hadronic width 
suppressed by Kqqj^/M^ ~ 10^^. 

Less inclusive variables that characterize Z decay to hadrons give rise to nonperturba- 
tive effects suppressed by smaher powers of Aqcd/Mz HSl USHl USHl II3Z1 II3H1 
Furthermore, these corrections often become even more important in corners of phase 
space where hadronization effects are significant, such as in the thrust distribution very 
near T = 1. It has been conjectured that the enhanced nonperturbative effects to many 
event shape distributions have a universal form with a single nonperturbative parameter 
|34| I138| I139| I140| I14H I142j . These arguments are based on analysis of renormalon ambigu- 
ities in the QCD perturbation series and on the behavior of resummed perturbation theory. 
The conjectured relationship between the nonperturbative corrections to event shape dis- 
tributions has recently been tested experimentally |35j . 

In Ref. j32j . the enhanced nonperturbative effects that occur for the jet energy distri- 
bution in corners of phase space were studied using effective field theory methods. This 
approach uses the fact that very low momentum degrees of freedom which contain the non- 
perturbative physics couple to the degrees of freedom with energies of order Mz via Wilson 
lines. Nonperturbative effects have been extensively studied previously jl37j using factor- 
ization methods to divide the process into hard, jet-like and soft subprocesses |8H I143j . 
Nonperturbative effects are computed from the soft subprocess. The effective field theory 
approach is similar to the one based on factorization methods. In this chapter we elaborate 
on the work in [321 extend it to other shape variables. The enhanced nonperturbative 
effects are expressed in terms of weighted matrix elements of operators involving Wilson 
lines, where the weighting depends on the event variable being considered. Our hope in 
this chapter is to make the results of Ref. jl37j more accessible to the community of high 
energy theorists who are most familiar with effective field theory methods. 
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We study smeared distributions, allowing us to expand the nonperturbative effects in 
powers of Aqcd, and write them as matrix elements of Wilson line operators and their 
derivatives. The computations are similar to those of smeared distributions in the endpoint 
region in B decay — the point-by-point computation requires knowing the nonperturbative 
shape function, whereas nonperturbative effects in the smeared distributions can be written 
in terms of Ai,2 provided the smearing region is large enough. 

For pedagogical reasons we start with a detailed treatment of the jet energy Ej in Z 
decay to two jets, where the jets are defined as Sterman and Weinberg did in their original 
work on jets in QCD |144j . We spend considerable effort on this variable because the 
theoretical expression for its enhanced nonperturbative corrections is simpler than for other 
more phenomenologically interesting variables like thrust. At lowest order in perturbation 
theory, the Z boson creates a quark and an antiquark, each with energy Adz/2, and so the 
jet energy distribution is equal to 



where r2-jot total two-jet rate at lowest order in perturbation theory. This leading 

order theoretical expression for the jet energy distribution is singular at Ej = Mz/2. 
Furthermore, the leading perturbative and nonperturbative corrections are also singular 
at that kinematic point. However, a non-singular quantity that can be compared with 
experiment without any resummation of singular terms is obtained by smearing the jet 
energy distribution over a region of size A that contains the lowest order partonic endpoint at 
Ej = Mz/2. The leading nonperturbative correction to this smeared energy distribution is 
suppressed by Aqcd/A. So, for example, with A ~ 10 GeV the nonperturbative corrections 
are expected to be of order 10%, roughly the same size as perturbative corrections, and an 
order of magnitude larger than the order Aqqd/Mz correction expected in the complete two 
jet rate. We argue that for Ej very near Mz/2 it is not possible to capture the dominant 
nonperturbative effects simply by shifting, Ej — > Ej — fi^p, in the perturbative expression 
for dr2_jet/d-E'j (where //np is a nonperturbative parameter of order Aqcd)- 

In the next section, we derive an expression for the leading enhanced nonperturbative 
correction to the smeared jet energy distribution for two jet events using methods from 
soft-collinear effective field theory (SCET) [HHl EOl ESI 1^ • This correction is given by the 
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Figure 5.1: Sterman- Weinberg jets. An event is characterized as a two-jet event for specified 
cuts f3, 6 if two cones of half-angle 6 contain most of the energy in the event, with no more 
than (3Etot of energy outside the cones. 

vacuum expectation value of a nonlocal operator involving Wilson hues. Perturbative order 
as corrections to this variable are derived in Appendix 15.41 

Section f5.3l discusses the leading nonperturbative corrections for thrust, jet masses, the 
jet broadening variables, the C parameter and energy-energy correlations. In agreement 
with Ref. |137j we find that the correction to jet mass sum and thrust are related. How- 
ever, without additional model-dependent assumptions we do not find that the enhanced 
nonperturbative corrections to the C parameter and jet broadening variables can be related 
to those for thrust and the jet masses. We compare the level of our understanding of the 
enhanced nonperturbative effects in these variables. 

5.2 Operator Product Expansion For The Two Jet Energy 
Distribution 

The nonperturbative corrections to the jet energy distribution for Z decay to two jets, 
dT2-jet/dEj near Ej = Mz/2 are computed in this section. The perturbative corrections 
will be discussed in Sec. 15.41 The results are given for the Sterman- Weinberg jet definition, 
where a cone of half-angle 5 contains a jet if the energy contained in the cone is more than 
-E'cut = f3Mz , as illustrated in Fig. 15.11 We take the cone half-angle 6 and the dimensionless 
energy cut variable /? to be of order a small parameter A, and compute in a systematic 
expansion in powers of A. We are interested in the jet energy distribution within a region 
A of Mz/2, where > A > X^Mz- For example, A ~ XMz- 
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SCETi is the appropriate effective field theory for the kinematic region of interest, and 
will be used for the derivation of the nonperturbative corrections to dr2-jet /d-Ej near Ej = 
Mz/2. It is convenient to introduce two lightlike vectors n and n which satisfy n*^ = n*^ = 1 
and n = — ii. Four-vectors are decomposed along the n, n and perpendicular directions: 
V = {y+,V-,VA_) where V'^ = n ■ V , V = n ■ V and = V^' - - V-n^'/2. 

For the problem of interest, SCETi contains n-collinear, n-collinear and ultrasoft degrees 
of freedom |119j . The n-collinear and n-collinear degrees of freedom have typical momenta 
that scale as 

~M^(A2,1,A), pf^ ^ Mz{lA\\). (5.2) 
and the ultrasoft degrees of freedom have momenta that scale as 

P„~Mz(A2,A2,A2). (5.3) 

We take A ^ \J Aqcd I^z which implies that the typical "off-shellness" of the ultrasoft 
degrees of freedom, ~ M|A^ ~ ^qcD' ^y the QCD scale while the typical "off- 

shellness" of the collinear degrees of freedom, ~ M^A^ ~ M^Aqcd, is much larger than 
Aq(-,j-,. Hence the collinear degrees of freedom can be treated in perturbation theory. 

Cone algorithms for jets, like that of Sterman and Weinberg, are ambiguous at higher 
orders in perturbation theory |145l I14tjj . This arises when there is more than one way to 
assign a particle to a particular jet. However, in this section we work to lowest order in 
perturbation theory, where the events consist of two almost back-to-back jets plus ultrasoft 
degrees of freedom. Since the cones are well separated, there is no ambiguity in assigning 
partons to the jets. 

The nonperturbative effects we are after are characterized by matrix elements of opera- 
tors composed from the ultrasoft degrees of freedom. In Z decay into two jets, the jets are 
almost back-to-back, and n is chosen along one of the jet directions. The degrees of freedom 
in the two jets are then represented by n-collinear (for the antiquark jet) and n-collinear 
fields (for the quark jet). In this section we work to lowest order in perturbation theory in 
the collinear fields. Hence we match the weak neutral current in full QCD onto the effective 
theory at tree level, 

f = \lnWn\V^^lin\. (5.4) 
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where = gvj± + 9Al±lb involves the vector and axial couplings of the Z boson. The 
fields and ^„ are collinear quark fields in the n and n directions and we have adopted 
the convention 

U^) = Y.e-'^-%n,p{x)^ (5.5) 
p 

where the label momentum p contains the components of order 1 and X, n ■ p and p_L, and 
the order components are associated with the space-time dependence of the fields. The 
Wilson lines Wnfi are required to ensure collinear gauge invariance ^J. Since in this section 
we work to lowest order in QCD perturbation theory, they play no role in the analysis and 
can be set to unity. 

The typical momenta of the partons in the jets are of the order of the collinear momenta, 
Eq. (|5.2|) . where the overall scale of their momentum is set by Mz- However, it is possible 
for the jets to contain partons with momenta that have an overall scale that is much less 
than Mz- Because of the sum over all values of p in Eq. (|5.5|) . such partons can still be 
represented by collinear fields. The interaction of n-collinear fields among themselves is 
given by the full QCD Lagrangian, and so the hadronization of n-collinear partons into a 
jet is described by the full theory. 

The Lagrangian of the effective theory does not contain any direct couplings between 
collinear particles moving in the two different lightlike directions labeled by n and n; how- 
ever, they can interact via the exchange of ultrasoft gluons. It is convenient to remove 
the couplings of the collinear degrees of freedom to the ultrasoft ones via the field redefini- 
tion ins]: 



where An is an n-collinear gluon field and 



Yn{z) = Pexp 



^9 



POO 

/ ds n • Au{ns + z) 
Jo 



(5.7) 



denotes a path-ordered Wilson line of ultrasoft gluons in the n direction from s = to 
s = oo. This is the appropriate field redefinition for outgoing collinear fields, since if a 
factor of exp(— es) is inserted in the integrand to decouple the interactions at late times, one 
reproduces the correct ie prescription for the collinear quark propagator. For annihilation 
which contains incoming collinear particles Yn is from s = — ootos = and the daggers 
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are reversed in Eq. H5.6() . An analogous field redefinition with n ^ n removes the couplings 
in the Lagrangian of ultrasoft fields to the n-collinear fields. 
The differential decay rate for Z decay to two jets is 



final states t 

X (27r)^5^(pz - PJ„ - PJn - ku), 



(5. 



where the sum over final states includes the usual phase space integrations and e is the 
polarization vector of the decaying Z boson. Since after the field redefinitions shown in 
Eq. 1)5. 6|) . there are no interactions between the ultrasoft and collinear degrees of freedom, 
the matrix element factorizes, and at lowest order in perturbation theory in the collinear 
degrees of freedom, 



AV 1 d^p, d^p,- ..^:2 

1 /„ ai 



Xo, 



(5.9) 



In Eq. ((EH), lA^Syp is the square of the Z ^ qq decay amplitude averaged over Z polar- 
izations and summed over the quark and antiquark spins and colors, T(T) denotes time- 
(anti-time-) ordering, Nc is the number of colors, and we have explicitly displayed the color 
indices on the ultrasoft Wilson lines. 

The derivation of Eq. 1)5.9(1 in many ways parallels the use of the operator product 
expansion to compute the deep inelastic scattering cross-section, or the rate for inclusive 
semileptonic B decay. There is, however, one important distinction. The sum over final 
states in deep inelastic scattering and B decay is a sum over a complete set of color-singlet 
hadron states. In Eq. ()5.8|) . one is summing over a complete set of jet and ultrasoft states. 
These are a complete set of partonic states, and are not necessarily color-singlet states. 
In fact, unitarity would be violated if one separately imposed the color-singlet condition 
on each of \Jn), \Jn) and {X^)- The derivation of Eq. (|5.9|1 is valid to the extent that 
the sums over partonic and hadronic states are equivalent. In jet production, the color of 
the fast quark that turns into a jet is eventually transferred to low-energy partons during 
the fragmentation process. The low-energy partons communicate between the different 
jets, and make sure the whole process is color-singlet. The assumption is that this color 
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recombination does not affect the decay rate at order Aqcd/-^z- 

To calculate dr2-jot/d-E'j we integrate Eq. (|5.9|1 over the allowed values of the quark 
and antiquark three-momentum with the factor S{Ej — p^) inserted. This corresponds to 
choosing the quark jet as the "observed" jet. If one does not distinguish between quark and 
antiquark jets then Eq. (|5.9jl still applies since the value of dr2_jet/di?j when the "observed" 
jet is an antiquark jet is the same. It is convenient to work in the rest frame of the decaying 
Z, pz = (Mz, Mz, 0_l), and align n with the quark three-momentum pq. The decomposition 
of the quark's four- momentum in terms of label and residual momentum, pq = pq + kq, has 
the form p^ = p^ + with p^ = 0, Pq± = 0. (Note this means that pq± = kg_L = and 
k~ =0.) Hence the phase space integration over quark three- momentum becomes 



d^p, _ 1 



Y,P+ dk+. (5.10) 



(27r)32p0 4(27r)2 ^ 

For the antiquark's four-momentum the decomposition into residual and label momentum 
is Pq = k'^ , p'^ =pq + k^ and Pq± = Pq± + kq_|_ . One cannot set pq± = by a choice of n, 
since n = — n, and ii has already been fixed by the direction of the quark jet. 

Expressed in terms of label and residual momenta the phase space integration over 
antiquark three-momentum is 



d^Pq ^ f d% 



(27r)32p0- 



Pq 



r dk^d%^ 1 

^1 2(2vr)3 p^.- ^^-^^^ 



Pq 



Here the delta function 5[{pq + kq)'^) = 5{pqkq — p?_|_) was used to do the k^ integration 
setting kq = p^_\_/pq ■ At leading order in the SCET expansion parameter A the invariant 
matrix element ■M.f'f only depends on the label momenta p+ and p^. In terms of label and 
residual momentum the energy-momentum conserving delta function becomes: 



&^{PZ -Pq-Pq- ku) = 26{pz -Pq -k^ )6{p^ - pt -pt - kt)^'^{Pq± + kn±) 



(5.12) 



The relation kq = Pq±/pq and the Kronecker delta that sets Pq± to zero imply that kq = 0, 
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and so this variable does not appear in the penultimate delta function in Eq. ()5.12|) . 
Using these results gives: 



dEj SMz{2'k) 



T[Yna %11(0) 



X ^ (0| f[Ynd T(0) \Xu{ku)) (Xu{ku) 

X ^ (o| r[y„rf ^y^t^ '^](o) \xm) {xu{ku)\ T[Yna 'yI '^](o) 



(5.13) 



We write the remaining delta function as the integral: 



V 2 2 



du 
2^ 



exp 



(f 



(5.14) 



At this stage the collinear degrees of freedom have been integrated out, and the matrix 
elements above, which involve only ultrasoft degrees of freedom, are evaluated at leading 
order in the SCET expansion parameter (i.e. A — > 0). Recall that the Sterman- Weinberg 
jet criteria restrict particles outside the cones used to define the two jets associated with the 
quark and antiquark to have energy less than Ecut which we are taking to be order \Mz- 
In the limit A ^ this energy cut becomes much larger than a typical component of an 
ultrasoft four-momentum. Hence, for the matrix elements of these operators, -Ecut should 
be taken to infinity and does not restrict these matrix elements. Similarly the cone angle 
is taken to be of order A while the typical angle between components of ultrasoft momenta 
is order unity. Thus the cone angle should be taken to zero in the effective theory that 
contains only ultrasoft degrees of freedom and so there is no restriction on the ultrasoft 
states that are summed over in Eq. (|5.1.S|) . 

Using the exponential dependence on ku to translate the anti-time ordered product to 
the space-time point un/2, and then using the completeness relation to perform the sum 
over all ultrasoft intermediate states, we find for the jet energy distribution: 



dr^ 



^^ = ^'^US{Mz/2-EjI 



(5.15) 
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where the shape function S is defined by |136j 



S{k) = ^ l ^e'"^ (o| f[Y^, ^Yl '^]{un/2)T[Yna ''](0) |o) , (5.16) 
and the total two jet Z-decay width at lowest order in perturbation theory is 



r 



having implicitly summed over spins and colors. The n-directed and n-directed ultrasoft 
Wilson lines commute since (sin — S2n)'^ = — 4siS2 < 0, and the gauge fields in the Wilson 
lines are space-like separated. 

In this derivation we chose the jets to be composed entirely of collinear degrees of 
freedom. This is appropriate since jets are confined to narrow cones. For example, the mo- 
mentum of any massless particle in the quark jet satisfies P^, which is the appropriate 
scaling for collinear particles in the n direction. However, it is possible to repeat the above 
derivation allowing ultrasoft degrees of freedom to be inside a jet. Then instead of inserting 
5{Ej-p°g) into Eq. one inserts 6{Ej-p°g-klj), where k^j = (A;+^/2)[l+C'(A)] denotes 
the total ultrasoft energy inside the quark jet. Using the delta functions in Eq. (|5.12j) we 
obtain again Eq. (|5.1.S)) . with k^ in the final delta function now denoting the total ultrasoft 
momentum outside the quark jet. However, as mentioned previously, at leading order in A 
the cone angle of the jet shrinks to zero, and one recovers the previous result. 

It is possible to remove the time- and anti-time-ordering completely in the definition of 
the shape function S. Using the results from Appendix I5.A. II our expression for the shape 
function becomes 



(0) 2 



NcMz, 2 , 2 



S{k) = ^e^'- (o| [y^ dYl ^]{un/2)[Yna ^^n' c](0) |o) 



[yV dYl ^<5(A; + in ■ d/2)[Yna ""Yn c] 



(5.18) 



, 



where the overline denotes an anti-triplet Wilson line. 



Since in the kinematic region of interest Mz /2 — Ej is much larger than n ■ d acting on 
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ultrasoft gauge fields it is appropriate to expand the delta function above which gives 



S{Mz/2 - Ej) = 6{Mz/2 - Ej) + 5'{Mz/2 - Ej) (0| Oi |0) 
+ \5"{Mz/2-Ej){^\02\Q) + --- 



(5.19) 



where 



1 



^ n ne 



in ■ d 



Y Cy " 



Nc 



-Tr 



y-' 



in • 



y. 



(5.20) 



The simple form for the operators Om arises because the variable Ej is totally inclusive on 
the "unobserved" antiquark jet. 



The formula for dr2-jct/d-Ej is 



dr 



2-jet 



dEj 



T 



(0) 
2-jet 



5iMz/2 - Ej) + 5\Mz/2 - Ej) (0| Oi |0) + • • • 



(5.21) 



The delta function term in Eq. 1)5. 19() simply reproduces the leading perturbative formula 
for dLT^2-]ct/^-^J while the higher-order terms contain the effects of nonperturbative physics. 
The derivation presented here assumes the observed jet is the quark jet. A similar derivation 
in the case where the antiquark jet is observed gives operators 



Or. 



1 



Tr 



y. 



t f in - D 



y. 



(5.22) 



Since the vacuum expectation values of Om and Om are equal by charge conjugation, our 
results also hold in the case where one does not distinguish between quark and antiquark 
jets. 

We define the matrix elements using dimensional regularization with MS subtraction so 
that in perturbation theory the vacuum expectation values (0| Om |0) are zero. 

Note that O2 is a very different operator than Oi so it is not possible to capture the 
effects of nonperturbative physics for \Ej — Mz/2\ ~ X^Mz simply by taking the lowest 
order perturbative formula in Eq. ()5.21|) and shifting Ej by a nonperturbative parameter 



^More correctly the differential cross section dr2-jct/di?j smeared over a region A of energy (that contains 
Ej = Mz/2) with A of order X^Mz- 
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Unp, that is, Ej Ej — ^Unp. This ansatz results in the shape function 



S{Mzl2 - Ej) = 6{Mz/2 - Ej) + 8'{Mz/2 - Ej)fi^p + ^6"{Mz/2 - Ej)f,% + ■■■ (5.23) 



where the series of derivatives of delta functions has coefficients that are simply related by 
(0| Om |0) = (0| Oi |0)'", which is not correct. 

For \Ej — Mz/2\ ~ X'^Mz all terms in the series of Eq. 1)5. 19() are equally important. 
However for \Ej - Mz/2\ ~ A > X^Mz the vacuum expectation value of Oi provides 
the leading order Aqcd/A nonperturbative correction. In this kinematic region the shift 
Ej — > Ej — /inp, with ^np = (0|Oi|0), correctly captures the most important effects of 
nonperturbative physics. 

We have focused on nonperturbative effects that are enhanced in the region near Ej = 
Mz/2. If one considers a variable like the average value of the jet energy over the entire 
allowed phase space, then there are sources of nonperturbative corrections that we have not 
considered. 

Using the results of ADDendix l5.A.21 the operator Oi in Eq. ()5.2U|) can be expressed in 
terms of the gluon field strength tensor |136j : 



Oi in Eq. 1)5. 24() vanishes if the ultrasoft gauge field is treated as a classical degree of 
freedom. Then the Wilson lines in Eq. 1)5. 24() are unitary matrices and the trace vanishes 
since the gluon field strength tensor is in the adjoint representation. Note that the vacuum 
expectation value of Oi can still be nonzero because of quantum effects. Usually operators 
involving products of gluon fields require renormalization. However, it is straightforward 
to show that Oi is not renormalized at one loop. For example, we have shown that the 
two-gluon matrix element of the operator Oi shown in Fig. 15.21 is identically zero, even 
before performing the loop integration. That is. 




-TT[Y^{in ■D)Yn] 




(5.24) 



(0|Oi|^^(ei,j.i)4(e,p2)>=0 



(5.25) 



after extensive algebra. 



129 




Figure 5.2: Two-gluon matrix element of Oi at one-loop order. Matrix elements of Oi such 
as this are identically zero in perturbation theory; hence, Oi is not renormalized. 



5.3 Enhanced Nonperturbative Corrections to Event Shape 
Variables 



There are a number of event shape distributions that are commonly studied in the literature. 
Conventionally, one defines a general event shape distribution da/de, where e is an event 
shape variable defined such that the region e — > corresponds to the two jet limit. Examples 
are e = 1 — T for thrust, e = B for jet broadening and e = C for the C parameter. 
Any event shape distribution in Z decay contains both perturbative and nonperturbative 
contributions. The perturbative effects can be computed as a perturbation series in as{Mz)- 
At leading order, only two-jet (i.e. qq) events contribute. Events with more hard partons 
are suppressed by powers of as{Mz)- In general, nonperturbative effects are suppressed 
by powers of Aqqd/Mz, but in corners of phase space where e <C 1 these nonperturbative 
effects become enhanced. Here we consider the region Aqcd ^ Mze <C Mz and focus on the 
enhanced nonperturbative contribution suppressed only by a single power of AqcD/iMze). 

Perturbative expressions for the jet variables considered in this section have been ex- 
tensively studied in the literature jl471 11481 11491 115UI 11511 1152j . Our main interest is in 
nonperturbative physics. Working to leading order in as{Mz), the dominant nonpertur- 
bative effects are corrections to the two-jet distribution. Nonperturbative corrections to 
higher-order processes are suppressed by additional powers of a<j(M^). We will compute 
the enhanced nonperturbative corrections to some commonly measured event shape distri- 
butions, just as we did for the jet energy distribution in Sec. O Recah for the jet energy 
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distribution the dominant nonperturbative correction came from expanding 



+ ••• 



(5.26) 



in Eq. ()5.13() to linear order in k'^ . The delta function from Eq. (|5.12|) sets 
we therefore find 



and 



dr 



(0) 

2-jet _ -p(O) 
- 2-jet 



(5.27) 



where 



T 



^nd ^ ne 



(0) xM)(XuiK 



(0) 0)A;+. (5.28) 



The jet energy distribution has the nice property that one can write (A;^) as the vacuum 
expectation value of an operator involving Wilson lines of ultrasoft gauge fields [namely, 
Eq. (|5.24p ]. For some shape variables this is not possible. However, expressions analogous 
to Eqs. (|5.27H5.28|) can be derived. 



5.3.1 Thrust 

First we consider the thrust distribution dF/dT where the thrust T is defined by 

M^T = max^ |t • pi| , (5.29) 

where t is a unit vector that defines the thrust axis. The maximum is taken over all possible 
directions of t, and the sum is over all final state particles. To the order we are working 
the thrust axis t can be set equal to the spatial part of the lightlike four-vector n used 
to define the collinear antiquark field. It is convenient to call this direction the z-axis. 
The thrust distribution is calculated analogously to the two jet distribution except that 
the delta function 6{Ej — p^) is replaced by 5 [MzT — \pq \ — — X^^t l^uol)) where the 
sum is over all ultrasoft particles. We adopt the same conventions as in the jet energy 
distribution so that the phase space integrals are again done using the delta function in 
Eq. H5.12|) . Decomposing the total ultrasoft four-momentum, ku = + ku \ into the sum 
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of the ultrasoft momentum from particles in the same hemisphere as the antiquark (type 
a) and the same hemisphere as the quark (type b) the thrust T can be written as 



''ua I 



MzT=\p^^\ + \pI\+^\K 

a 



(5.30) 



Now the delta functions in Eq. ()5.12|) set = pt = Mz, 



kq = 0. Thus we find 



where we have also used ku = k^f' + ku"* . Thus, 



-A;+, and 



(5.31) 



dT _ (0) 

^j, ^ 2-jet 



p(0) 
^ 2-jet 



5{l - T) - 5'{l - T) 



5{l-T)-5'{l-T) 



Mz 



m 

Mz 



(5.32) 



The thrust axis and the hemispheres are determined by the jet directions, and can be defined 
in terms of the label momenta of the quark and antiquark. Thus t and the hemispheres a 
and h are label variables. Nevertheless, because of the hemisphere condition on the ultrasoft 
momentum in Eq. ()5.32p . there isn't a simple formula expressing the correction in terms of 
the vacuum expectation value of an operator involving Wilson lines like the one in Eq. (|5.24() . 

In a region |1 — T| ~ A'^ the higher order terms in the ultrasoft momentum that were 
neglected in Eq. 1)5. 32() are important. Eq. (|5.32() is appropriate for a region 5T near T = 1 
that satisfies 1 » 3> A^, for example, 5T A. 



5.3.2 Jet Masses 

The squared jet masses are the squares of the invariant mass of all the particles in 
the two hemispheres a and b, defined by the plane perpendicular to the thrust axis. Two 
commonly used variables are the sum of jet masses, M| = {M^+M^)/M^, and the heavy jet 
mass = max(M^, M^)/M^. The jet masses are = {pq+kt^f and = {pq+k'^'f. 



More explicitly, 
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Mi 



(a)\2 



(5.33a) 
(5.33b) 



Recall that pg is aligned along n so that p^^ = 0. Also, the delta function in Eq. (|5.12|) 

sets Y*ql. = ^-iid =Pq = Mz- Then, working to linear order in the ultrasoft momenta, 

, , , (b)- 

so 



= Mzfci"^+ and = Mzk^^^ 



where 



dr 

dM| 



dr 

dM| 



^(0) 
2-jet 



^(0) 
2-jet 



5{M, 



' qMsV 



5{Ml)e - 5' (Ml) 



Mz 



(5.34a) 
(5.34b) 



of-) = (max (fei'^H, 4") 



(5.35a) 
(5.35b) 



Note that in the kinematic region where expanding to linear order in ultrasoft and residual 
momentum is appropriate, the nonperturbative corrections to the M| and \—T distributions 
are given by the same nonperturbative matrix element. The nonperturbative corrections to 
the and M^ distributions are different. 

Working to higher orders in k^/Mz-, the definitions of thrust in Eq. 1)5. 31() and of jet 
masses in Eq. ()5.33|) become different beyond linear order. However, the corrections to 
event shape distributions at higher orders in AQCD/(-/^ze) come not from expanding the 
argument of the delta functions used to define these variables to higher orders in ku/Mz, 
but rather from expanding these delta functions as power series in the ultrasoft momentum, 
as in Eq. (|5.26|) for the jet energy. So even at higher orders, the enhanced nonperturbative 
corrections, i.e. of order [AQCD/(Afze)]", n > 1, come from the leading-order correction to 
the argument of the delta function, which are the same for thrust and jet mass sum. So the 
enhanced nonperturbative corrections to thrust and jet mass sum are related to all orders 
in AQCD/(Mze). 




Figure 5.3: Determination of the thrust axis. To the order we are working, the quark 
and antiquark have momenta \pq\ = \pq\ = Mz/2. The antiquark then makes an angle 
9q = 2 |kq_L| /Mz with the z-axis, and the thrust axis t makes an angle 6t = \^qj_ \ /Mz with 
both the quark and antiquark. 



5.3.3 Jet Broadening 

Jet broadening variables Ba^b are defined by 

where the hemispheres a and h are defined as before, and t is the thrust axis. The jet 
broadening variables at order ku/Mz require knowing the thrust axis to order k^/Mz- The 
thrust axis t maximizes |t • pj|. 

The angle between pg and the z-axis is given by 

^ - (5.37) 



^ Mz/2 ' 
and the thrust axis t can be written as 

t = (0,-sin6'j,cos6lt). (5.38) 

By symmetry, 

= (5.39) 

which is half the size of 6q (see Fig. 15. 3|) . 

Now calculate Ipj x t| for each particle. To linear order in ku/Mz we find for the quark, 



134 



\Vq X t| = —smOt = 



(5.40) 



and for the antiquark, 



|P5 X t| = Ik^^l cosOt ^sm6^4 = — |— 



(5.41) 



For each ultrasoft particle a, the cross product k^, x t is given by the determinant 



X y 



— sin 9t cos 9t 



(5.42) 



Since sin^^ is already of order ku/Mz the cross product is, to linear order in k^/Mz, 



ka X t — (fc^, fc^,0), 



(5.43) 



so jka X t| = Ikaxl. Combining the contributions of each particle to the sum in Eq. 1)5. 36() . 
and using the delta function in Eq. (|5.12j) to set kq_|_ = — k„_|_, we obtain for the jet broad- 
ening variables (to linear order in ku/Mz)' 



Ba 



1 / |k«j 



2Mz \ 2 



•al. I I ) 



(5.44a) 
(5.44b) 



where the sum on a is over the ultrasoft particles in hemisphere a or b. 



One conventionally defines two other broadening variables as 



-Bmax = max(Sa,Sb) , 

Bsuia = Ba + Bf). 



(5.45a) 
(5.45b) 



The jet broadening distribution is 



dT 
dB 



r(o) 

J- 2-jet 



6{B) - 5'{B) 



B\ 



Mz 



(5.46) 
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for Ba, fe, sum, max; where (Of^) = Mz (B) is the matrix element of the appropriate quantity 
in Eqs. (|5.44l5.45jl . Nonperturbative effects in the jet broadening measures are not related 
to the jet energy or thrust. 

In this paper, we have assumed that the nonperturbative physics is completely described 
by ultrasoft degrees of freedom. It is possible that some of the subtleties associated with 
nonperturbative corrections to the jet broadening variables that have been discussed in the 
literature jl41j can be attributed to nonperturbative effects in the collinear sector, which 
we have not included. 



5.3.4 C Parameter 

The C parameter is defined as 

C = 3(AiA2 + A2A3 + A3Ai), (5.47) 
where Aj are the eigenvalues of 

e^-' = ^yEM (5.48) 

Mz Y IP*I 

and r,s = 1,2,3 are the space components of the momentum pj of the i^^ particle. 

The largest component of 9'^^ is 9^^. The quark and antiquark in the jets have z- 
momentum = — p5 = Mz/2 to the order we are working. Then, to linear order in 
ku/Mz, the eigenvalues of 9^^ are given by: 

det{9 - XI) = (1 - A) det(X - XI), (5.49) 

where / is the identity matrix and 
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Here the sums over a are only over ultrasoft particles. (The contributions from the quark 
and antiquark to these components of 9^^ are suppressed by another factor of 1/Mz, since 
|pg| = |pg-| = Mz/2.) 

The largest eigenvalue is Ai = 1, and the other two eigenvalues satisfy 



A2 + A3 



— E 

Mz ^ 



(5.51) 



Thus, 



The C distribution is then: 



" Mz ^ |k,| 



dC 2-jet 



6{C) - 6'{C) 



Mz 



(5.52) 



(5.53) 



where (Of) = Mz (C) defined in Eq. (|5.52j) . Like jet broadening, the C parameter distri- 
bution is not local on the ultrasoft fields, and the nonperturbative correction is not related 
to that for any of the above distributions. 



5.3.5 Energy-Energy Correlation and Jet-Cone Energy Fraction 

The angular correlations of radiated energy can be characterized by the one-point and 
two-point correlations jl53j . 

"^'^ j dTY^^5{^-^,)5[^l'-^j), (5.54b) 



where the sum is over all particles, and includes the terms with i = j. They are normalized 
so that 

dn— = r, (5.55a) 
dn'^^ = — . (5.55b) 
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The energy-energy correlation function P(cos x) is defined by 

f d^S 

P(cosx) = / dO'dO (cos X - cos 0nn' ) , (5.56) 

where Oqq' is the angle between vectors in the and fi' directions. 

The angular energy correlations Eq. ()5.54b|) were defined in Ref. |15,Sj for e'^e~ annihila- 
tion, and the solid angle was defined with respect to the beam direction. For unpolarized Z 
decay, there is no preferred direction, so dS/dJ7 is a constant, and d^S/dildO' contains the 
same information as the energy-energy correlation function P(cosx). One can, however, 
define distributions analogous to Eq. H5.54b|) where the solid angle is measured with respect 
to the thrust axis t. The one-point function is called the jet cone energy fraction J. 

The energy-energy correlation and the jet cone energy fraction both are proportional 
to S functions if one considers the leading order process of Z decay into a quark-antiquark 
pair: 

P(cos x) = J (cos x) 

1 (5.57) 
= 2 ^0 (^°® X - 1) + (cos X + 1)] ■ 

Ultrasoft emission (in two-jet events) changes the distribution in two ways: (a) by changing 
the energy or (b) by changing the solid angle of the emitted particles. At order k^/Mz, the 
change in energy can be neglected, because it does not shift the angles of the partons; thus 
there is no contribution proportional to 6'{cosx i !)> a-s for variables such as thrust. The 
angle between the quark and antiquark is [compare Eq. ()5.37|) ]: 

cose,,- = -l + 2^, (5.58) 

z 

and the angle of the quark or antiquark with respect to the thrust axis is [compare Eq. (|5.39|) ]: 

COS = - cos = 1 - (5.59) 

z 

where k_L is the total _L momentum of the ultrasoft particles. The shift in angle is second 
order in k^/Mz, and so to first order, there is no enhanced contribution near cosx = if- 
There are nonperturbative contributions at second order. 
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5.3.6 Classes of Observables 

The different observables we have discussed can be divided into classes, based on the extent 
to which their nonperturbative corrections are inclusive on the ultrasoft degrees of freedom. 

A class I observable is the jet energy distribution. The nonperturbative correction to 
the jet energy depends on {k^), where ku is the total ultrasoft momentum, so the jet energy 
distribution is totally inclusive on the ultrasoft fields. The derivation of nonperturbative 
corrections to the two jet energy distribution is not quite on the same footing as the deriva- 
tion of nonperturbative corrections to the B meson semileptonic decay rate, because of 
the additional assumption about the equivalence of sums over partonic and hadronic states 
discussed after Eq. 1)5. 9() . 

Class II observables are thrust and the jet masses Mgjj. The nonperturbative cor- 
rections to these variables require the ultrasoft momentum to be broken up into two 
parts, ku = kij^^ + ku \ corresponding to the contributions from ultrasoft partons in the 
two hemispheres. The hemispheres are chosen based on the jet directions, i.e., based on 
the collinear degrees of freedom. The momentum in each hemisphere can then be de- 
fined by integrating the ultrasoft energy-momentum tensor over the hemisphere at infinity 
|1,S7| I154| I155| I157j . The class II variables are not totally inclusive on the ultrasoft 
variables, but require them to be divided globally into two parts. Whether our derivation 
of the nonperturbative corrections for class II observables (e.g., the relation between jet 
mass and thrust distributions) is valid depends on the nature of hadronization in QCD. 
The ultrasoft fields end up inside final state hadrons. The final hadron can contain ultra- 
soft partons from different hemispheres, so the hadronic energy flow in each hemisphere 
does not have to equal the parton energy flow in each hemisphere. If the hadronic and 
partonic energy flows differ by order unity, the derivation of nonperturbative effects in class 
II observables is invalid. If, for a smearing region of size A the mixing of ultrasoft momenta 
between the two hemispheres during hadronization is an effect of order Aqcd/A, then its 
impact on class II observables is the same size as k'^ effects, which are one higher order than 
the terms we have computed. 

Class III observables are the jet broadening measures -Ba,6,sum,max and the C parameter. 
These depend on knowing the individual ultrasoft momenta of each parton. This appears 
to be a notion that cannot be made rigorous in field theory. 
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5.3.7 Model-Dependent Relations Among Event Shape Variables 



Nonperturbative corrections to event shape distributions have been considered extensively 
in the Uterature in the past. For example, in the work of Ref. 11371 . nonperturbative shape 
functions were derived for thrust and jet mass distributions. The enhanced nonperturbative 
corrections to these distributions are given by first moments of these shape functions, and 
the results in sections 15.3.11 and 15.3.21 are in agreement with Ref. jl37j . 

The derivations of the enhanced non-perturbative corrections in this section have only 
relied on the fact that they arise from matrix elements of ultrasoft operators. It is insightful 
to understand what further conditions have to be imposed to reproduce other proposed 
relations amongst nonperturbative parameters for event shape distributions [341 11381 11391 

nm. 

As an example, consider the C parameter, for which the nonperturbative matrix element 
was defined as 

^Of)=3^E^)- (5-60) 

For on-shell soft gluons collinear to the antiquark or quark jet (i.e., in hemisphere a or b, 
respectively), A;('*)+ <C fe^")" and /c^^^" <C A;(^^+. This imphes that 




"13 

/9 



2 (4'^)+ 4^) 



This leads to 

(Of)coii = 6(Of). (5.62) 

To take into account that ultrasoft gluons can also be radiated at a finite angle, one can 
impose the condition that the matrix elements of Of and Of are given by the one-gluon 
contribution in perturbation theory, performing the angular integrals in the phase space at 
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a fixed value of |k_L|. Under this assumption, the matrix element of Of is given by: 

|k,xp\ _ / rh. . „ ik,|2 1 ' 



d6lsin6l- 



|ka| /i_giuon \^o (|k^|/sin0)sin^6'/ ^5^3) 

= vr(|k^|), 

where the factors of sin^ from the phase space, from the relation |k_|_| = |k| sinO, and from 
the squared amplitude for one gluon emission have all canceled out to give the final result. 
For the matrix element of O^, we calculate: 

. „|k| (1 — cosO) \ 



I l-gluon \ Jo 



sin^ Q / 

= J\\eMl^Z^\ (5-64) 

= 2(|kx|). 

This leads to the result 

Svr 

(Ol )l-gluon = ^(Oi )l-gluon- (5.65) 

Given the assumptions that have to be made to obtain Eq. (|5.65)) (or analogous relations 
based on higher orders in perturbation theory) , it does not seem likely to us that there is a 
simple analytic nonperturbative relation between (Of ) and (Of)- 



5.3.8 Comparison with the Data 

Predictions for event shape variables have been compared with experimental data in Refs. [351 
I128j . Nonperturbative corrections have been included using the ansatz that their effect on 
distributions for shape variables is described by shifting the variable by cj/inp/-E'cm in the 
perturbative formula for the distribution: 

Here c/ is a constant that depends on the kinematic variable /, ^i^^ is a universal nonpertur- 
bative parameter, and E^m is the center-of-mass energy. An analysis in perturbation theory 
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(similar to what was done in section I5.3.7|) provides simple relations between the c's for 
some of the event shape variables. We have found that, provided one is not in a kinematic 
region that is extremely close to the partonic endpoint (i.e., the shape function region), c for 
1 - T and M| are the same. However, we argued that c for other parameters like the heavy 
jet mass and C are not connected to c for thrust. Some experimental evidence for this can 
be found in the analysis of Ref. [HSl- For 1 — T and the jet mass sum^ a simultaneous fit 
for Os and /inp under the assumption that c takes on its conjectured values (see Fig. 15.4(1 
yields values of /inp that are close to each other, and values of that are consistent with 
other extractions of the strong coupling. However, Ref. [HSl finds that /^np for the heavy 
jet mass, C parameter, and jet broadenings are not related to /inp for thrust in the way 
that the analysis based on perturbation theory suggests, and, furthermore, a fit to these 
variables does not yield a value of that is consistent with other extractions. 

5.4 Perturbative corrections to dV/dEj 

Neglecting order Og corrections to the nonperturbative effects proportional to Oi, pertur- 
bative corrections to dr2-jet/di?j in Eq. (|5.9|) can be calculated in full QCD using standard 
methods. In this section we first review the computation of perturbative 0{as) corrections 
to the total two-jet rate r2-jet and then compute the jet energy distribution dr2-jet/di?j at 
order a^. We work in d = 4 — e dimensions to regulate infrared, collinear and ultraviolet 
divergences that occur in contributions to the differential decay rate. The jets are defined 
using the Sterman- Weinberg criteria which involve an energy cut l3Mz and a cone half-angle 
6. Corrections suppressed by UsP and asS are neglected. 

5.4.1 Two Jet Decay Rate 

Using the Sterman- Weinberg definition of jets, there are three contributions to the two-jet 
rate at 0{as): 

(a) One quark and one antiquark each creating a jet; 

(b) One quark and one antiquark each creating a jet, plus a gluon with energy Eg < f3Mz', 
^Ref. 121] advocates the use of a modified i?-scheme jet mass to reduce sensitivity to hadronic masses. 
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Figure 5.4: Experimental fit for model of nonperturbative power corrections to event shape 
variables. The DELPHI collaboration tested the Dokshitzer- Webber model for nonper- 
turbative power corrections to event shape variables in e~''e~-annihilation, fitting to the two 
model parameters uq and as{Mz)- The data suggest a fairly poor fit for most variables, 
except for the thrust and jet mass sum, which are most nearly consistent with each other 
while also agreeing with independent extractions of as{Mz)- These findings are consitent 
with the theoretical predictions in Sec. 15.31 
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(c) One quark and one antiquark each creating a jet, plus a gluon with energy Eg > l3Mz 
inside one of the jets (within an angle 6 of the quark or antiquark). 



Contribution (a) is simply the rate T{Z — > qq). The tree and virtual gluon graphs give 
the amplitude: 

Mz^gg = e^ipzKipg)T^'vaiPg) (^1 + , (5-67) 

where the color index a is summed over values a = 1, . . . , Nc, and Cp is the Casimir of the 
fundamental representation. Explicit computation of the one-loop vertex correction gives, 



1 - ^ ? ,n (Z%i?l) - 4 + - i In^ (=5Vi^) + ^„ f Z5VZ£ 



Integrating the square of the amplitude over the d dimensional two body phase space gives: 
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7^ 
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Ml 



(5.69) 



where il^ is the total solid angle in d dimensions. The 1/e poles will cancel out against 
divergences from the real gluon emission graphs. We do not need to expand the bracketed 
prefactor in Eq. ()5.69|) in powers of e because the identical factor will appear in the real 
gluon graphs. 

Contributions (b) and (c) come from integrating the square of the amplitude for real 
gluon emission, Z qqg, over the three-body phase space in d dimensions. We find for the 
terms that do not vanish as (3 and 5 go to zero. 



r 



Z^qqg 



glMl-'Nc Cp 
2567rS 



i'^^f'{-^y^2^^^^~^i9l+9\) 



r(-e) ^ 



(5.70) 



and 



Z-*qqg 



glMl-'NcCp 
2567r5 
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(2/3)- 



3e 13e2\r(-e)2" 



(5.71) 
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Adding together contributions (b) and (c), expanding in powers of e and converting to the 
MS scheme yields for the total rate for Z — > qqg in the two-jet region 



Z^qqg 



-{9v+9a) 



327r3 



-n 



3-e 



4 3 
e"^ e 



2, Ml 
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e /x^ 
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41n 2/3 In5-31n5 + 



13 llvr^ 



12 



(5.72) 



Finally, we add together the rates Tz^qq and Tz->qqg from Eqs. (|5.69|) and (|5.72|) . The 
e-dependent prefactors in brackets in the two equations are identical, as promised. The 
1/e-poles in the remainder of the expressions cancel out exactly (as do all the logarithms of 
M^//x), so we can set e = in the remaining finite parts, leaving 



r 



2-jet 



NcMz , 2 ^ 2^ 



1 + 



\2 



31n5-41n 2/3 ln6 



(5.73) 



which agrees with Sterman and Weinberg's original result |144j . 



5.4.2 Differential Decay Rate dFa-jet/d^^j 



We now turn our attention to the differential decay rate dT 2-jet /dEj. The contribution of 
Tz^qq to this rate is simply 



'-^^rz.q,s(E,-^\ 



dEj 



2 J 



(5.74) 



where Tz^qq is the total rate for Z ^ qq calculated to 0{as), which is given by Eq. 1)5. 69() . 

For the contribution of real gluon emission processes, we write the three-body phase 
space for this rate: 

^iff = 16Mz (2.1'^-- ^.-2^.~i^E,Et'dE2Et'd cos sm^"^ 6 

^ ' (5.75) 



X 5 



M^-2Mz{Ei+E2) 
2E1E2 



+ 1 — cos / 



5iEj-...)\M\\ 



where the 5{Ej — • • • ) defines Ej according to which partons actually go inside the jet. It 
is useful to split up the phase space slightly differently than for the case of the total rate: 
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(a) Gluon with energy Eg > f3Mz inside unobserved jet; 

(b) Gluon with any energy inside observed jet; 

(c) Gluon with energy Eg < (3Mz outside observed jet. 

These three regions exhaust the possible gluon energies and locations with respect to the 
jets. It is convenient to introduce the variable 

ej = ^-Ej (5.76) 

and focus on a region of ej near the origin with size of order l3Mz- 

For case (a), where a gluon with Eg > PMz is inside the unobserved jet, take Ei = Eg, 
E2 = Eq, so 6 is the angle between the gluon and antiquark, and Ej = Eg. Integrating over 
and Eg using the delta functions leaves an integral over Eg running between the limits 



and restricts ej to lie between 



S^PMz <ej< (5.78) 



Similarly, for case (b), where a gluon with any energy lies inside the observed jet, 
El = Eg, E2 = Eg, and Ej = Eg + Eg. Integrate over 6 and Eg using the delta functions. 
Then the limits of the Eg integral are 



and ej is restricted to the region 

Mz5'^ 



< ej < 0. (5.80) 



Before we proceed to case (c), note that the physical observable we actually want to 
calculate is the smeared distribution 



dr 



dEj 



dEjw^iEj)-^, (5.81) 
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where wa is a smooth function which smears the differential rate over a region of jet energy 
whose size is of order PMz- But the contributions to the rate from cases (a) and (b) have 
support only over a region of size S^Mz <C l3Mz near ej = 0. Consider smearing dT/dEj, 
or equivalently dF/dej, over a region near Ej = Mz/2 (ej = 0) of size of order /3Mz- Then 
w{0) - l/(5Mz, w'{Q) ~ l/(/3Mz)^ etc. Expanding, 



Since w'{^)/w{Qi) ~ 1/ f3Mz, and, for the contributions in cases (a) and (b), ej ~ 5'^Mz in 
the region where dV/dej is nonzero, the second term is suppressed by a power of 6"^ / (3 <^ 1. 
Thus only the first term is relevant.^ Keeping only the first term amounts to replacing the 
full dV/dej by 



However, integrating the contributions of (a) and (b) to dV/dej over all allowed values 
of ej simply gives their contribution to the total Sterman- Weinberg jet rate, that is, they 
will build up part of the term Vz-^qqg6{Ej — Mz/i) in dr2-jet/d-E'j. Since we have already 
calculated the total rate, we need not analyze cases (a) and (b) any further, as long as we 
can get the remaining contribution to the total rate from case (c). 

In case (c) we have a gluon with Eg < f3M anywhere outside the observed jet. Here 
El = Eg, E2 = Eg, and Ej = Eg. Writing out the formula for the rate explicitly. 




(5.82) 




(5.83) 



dr(^) 1 1 



nd-2nd-ie{ej)e[pMz{i - s^) - ej] 



dEj 16Mz (27r)2'i-3 
rl3Mz 

X dE„^ 




ic) 

Z^qqg 



2 



(5.84) 



The part of the amplitude that gives a contribution that survives as /3 ^ is 



M 




d-2 



{9v+9a) 



Mlpg ■ Pq 



(5.85) 



d-1 



{k ■ Pq){k ■ Pq)' 



^This argument assumes that the integral / dej ejdF/dej is finite, which can easily be shown. 
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Substituting Eq. 1)5.85(1 into the phase space, 



dEj 2567r5 \Mz/ d-1 ^ ^ 

X {gl + g\)ej^'^^e{ej)e[()Mz{l - 5^) - ej] In 



pMz - ej 



(5.86) 



The factor, {l/ej)ln[{(3Mz — ej)/(ej5^)], is singular as ej — > 0, and must be rewritten in 
terms of an integrable quantity. Use the "plus distribution": 

i3Mz fl3Mz 

dej f{ej)+ g{ej) ^ / dej /(ej)[g(ej) - ^(0)], (5.87) 



JO 

where / diverges at e j = and 51 is a test function finite at ej = 0. To replace / by /+, we 
would write 

/ dejfiej)g{ej)= dej f{ej)+ g{ej) + g{0) dej f{ej). (5.88) 
Jo Jo JO 

The second term amounts to replacing 

f{ej)^6{ej) I de'jfie'j). (5.89) 

But making this replacement in Eq. (|5.86() means writing a delta function 5{Ej — Mz/2) 
and integrating the differential rate over all allowed values of ej, which again just gives 
its contribution to the total Sterman- Weinberg jet rate. Together with the contributions 
from (a) and (b) this gives the one loop contribution to 6{Ej — Mz/2)T2-jct- Only the plus 
function piece gives a deviation of the jet energy distribution away from Ej = Mz/2. The 
final result for the differential rate to O(a^) is: 



dr2-jct _ i- 
~dE7 ~ 



(^Ej - ^)T2.jet 



+ ^^^^^^i9'v + 9\)0{ej)Bmz-ej) 



where the total rate r2-jet is given by Eq. (|5.73l) . 



1 , fpMz - ej 
— m 



(5.90) 
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5.4.3 First Moment of the Jet Energy Distribution 



As an application of the above result consider the first moment of the jet energy distribution, 
defined by 



Mi(/) 



Mz 



dEj 



1 dr 



2-jet 



dE 



J 



3l 

Mz 



(5.91) 



2 'fiBMz V-^2-jct 

Using the expression in Eq. (|5.21|1 for the nonperturbative correction and in Eq. (|5.90|) for 
the order ag perturbative correction to the jet energy distribution gives 



Mi(/) 



IT 



/log(7^)-(l 



/)log(l-/) 



(0|Oi|0) 
Mz 



for / < 1 and 



(5.92) 



(5.93) 



for / > 1. Note that the order Os contribution to Mi{f) is independent of / for / > 1. 
This occurs because the perturbative correction vanishes for Ej < Mz/2 — (3Mz- 

In Fig. 15.51 we plot Mi{f), for / < 1. For this figure the value of the energy cut is 
f3 = 0.15 and the cone half-angle is (5 = 15° and the vacuum expectation value of Oi is 
set equal to 500 MeV. We evaluate a<j at the scale (3Mz and find with these parameters 
that the order corrections reduce the two jet rate by about 16% from its tree level value 
lending support to the validity of perturbation theory for the values of the cone angle and 
energy cut used in Fig. 15.51 



5.4.4 Perturbative Corrections in the Effective Theory 

Although we have used full QCD to calculate the jet energy distribution it is possible to 
do the computation in the effective theory. Here we briefly discuss how that computation 
would proceed. 

The full theory amplitude for Z ^ qq is reproduced in SCET by the Wilson coefficient 
in the current matching: 

iqcD = [CnWn]Cifi,pg ■ p^)Tf'[WlU (5-94) 
where there is an implicit sum over label momenta and the matching coefficient C{iJ,,pq -pq) 
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f 



Figure 5.5: First moment Mi(/) of the jet energy distribution. The black sohd curve 
shows the perturbative contributions only, while the red dashed line represents the moment 
including the nonperturbative contribution. The figure corresponds to /? = 0.15, 6 = 7r/12, 
(0|Oi|0) = 0.5 GeV, and we have evaluated the strong couphng constant at the scale 
= fiMz. 

can be read ofF^ from Eqs. (|5.H7|) and (|5.68|) : 



1 + 



agCp 



27r 



-4 + 



vr 
12 



In^ 



-2Pg • Pq 



3, 
+ - In 
2 



-'^Pq-Pq 
^2 



(5.95) 



and the UV renormalization factor for the current in the effective theory is 



Zv = l + 



OsCf 

27r 



3 2, 
- + - In 



-2pg • Pq 
^2 



(5.96) 



Note that both the renormalization factor and the matching coefficient depend on the label 
momenta for the quark and antiquark. For outgoing particles, the collinear Wilson lines are 
defined as 



Wniz) = P exp 



ig I dsn ■ An{ns + z) 





(5.97) 



and one must include collinear gluons produced by a Wilson line in real gluon emission 



^The matching coefficient is just given by the finite part of the full theory matrix element {qq\j^\Q) 
because the full theory current has no anomalous dimension, so the 1/e poles are pure IR divergences, which 
must cancel out in the matching condition. The loop graphs in the effective theory contributing to this 
matrix element are zero in dimensional regularization, so the finite part of the matching coefficient is just 
the finite part of the QCD matrix element, given by Eqs. 115.671 and II5.58II . while the infinite parts become 
the UV counterterm in the effective theory |82m.158| . 
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to get the correct Z — > qqg amplitude. We find that the perturbative expressions for the 
two jet rate presented in the previous sections are reproduced by the effective theory if we 
cah any particles inside the observed quark jet n-collinear particles and all other particles 
n-collinear. In the effective theory, ultrasoft gluons in the final state contribute zero in 
perturbation theory and appear only in the nonperturbative shape function. A similar 
result holds for deep inelastic scattering |159j . 

5.5 Concluding Remarks 

We have studied nonperturbative effects in Z decay to hadrons using soft-collinear effective 
theory (SCET). The jet energy distribution for two jet events has enhanced nonperturbative 
effects when the jet energy is near Mz /2. These nonperturbative effects can be expressed 
in terms of the vacuum expectation value of operators involving Wilson lines. The Wilson 
lines arise from the coupling of ultrasoft gluons to collinear degrees of freedom in the jet. In 
Sec. 15.41 we derive the order perturbative corrections to the jet energy distribution and 
discuss the implications of perturbative and nonperturbative physics on the first moment 
of this distribution. 

For a region of \Ej — Mz/2\ that is of size A, the leading nonperturbative corrections 
to the jet energy distribution are of order Aqcd/A when A is large compared to Aqcq. In 
this region they can be characterized by the vacuum expectation value of a single operator 
involving ultrasoft fields which provides a contribution to the jet energy spectrum that is 
proportional to 5'(M^/2 — Ej). For multijet events, a similar analysis holds; however, an 
additional operator analogous to Oi but involving adjoint Wilson lines occurs for a gluon 
jet ISU. 

When A ~ Aqcdj one is in the shape function region, and the functional dependence 
on Ej is much more complicated. While we focused mostly on the kinematic region where 
Mz ^ A ^ Aqcd, it was shown that in the shape function region, it is not possible 
to capture the effects of nonperturbative physics by introducing a single nonperturbative 
parameter /i„p and shifting Ej Ej — ^np in the perturbative formula for the jet energy 
distribution. 

The jet energy distribution has the special property that it is totally inclusive in one 
of the jets, and hence expressions for nonperturbative effects can be derived using operator 
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methods that are similar to those used for the endpoint region in inclusive semileptonic B 
decay. Other event shape variables (e.g. thrust, jet mass, jet broadening) have nonpertur- 
bative effects that are enhanced in the partonic endpoint region. We discussed the extent 
to which these effects can be understood using field theoretic methods in QCD. 

5. A Appendix: Properties of Wilson Lines 

In this section we derive some useful properties of the ultrasoft Wilson lines introduced in 
Eq. (jEH). 

5.A.1 Relations Between Triplet and Anti-triplet Wilson Lines 



Consider the time- and anti-time-ordering of the Wilson lines in the shape function S defined 
in Eq. ()5.16() . For Yn, the path ordering is the same as time ordering and so T [Yn] = Yn- 
Consider writing Yn as the product of N infinitesimal integrals over path segments of length 
ds, 



J,gAids\ b 



(5.98) 



with the subscripts on the ultrasoft gauge fields denoting their space-time location along 
the path of integration. Taking its adjoint 



b_ ( ^-igAids\ bi l^-igAN-ids\ fejv-i /gi9Ajvrfs\ b 

'bN-i 



(5.99) 



Time ordering this expression, 



T 



yt b 

na 



^igA^ds 



Q-igA%ds 



b 



b ( Q-igAM-ids 



-igA\ds\ fei 



bN-l \^ I OjV-2 



(^-igAlds\" _yb 



(5.100) 



where the overline denotes an anti-triplet Wilson line. (Recall that the generators in the 3 
representation are minus the transpose of those in the 3.) Similarly, 



T 



Y„ 



— 



n 



yt h 

na 



yt h 

na 



(5.101) 



From these results, Eq. (|5.18|1 follows. 
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5. A. 2 Oi in Terms of the Gluon Field Strength 



We can express the operator Oi in terms of the gluon field strength tensor as written in 
Eq. (|5.24j) . It is convenient for this purpose to generalize the expression for the ultrasoft 
Wilson line to 



Yn{z;b,a) = Pexp 



f 

ig ds n ■ A{z + ns) 

J a 



(5.102) 



so that with a = and 6 = cxd we recover the standard Wilson line used above, Yn{z) 
l^(z;0,oo). Differentiating along the n direction, 



n • d Yn{z) = ig ds Yn{z; oo, s) [n ■ ■ A\{z + ns)Yn{z; s, 0) 

^0 

/•oo 

= ig ds Yn{z; oo, s) [n ■ dzU ■ A — n ■ dzU ■ A + n ■ dzU ■ A] (z + ns) 
Jo 

xYfi{z;s,0) (5.103) 

/•oo 

= ig ds Yfi{z; oo, s) [n ■ dzU ■ A — n ■ dzU ■ A] {z + ns)Yn{z; s, 0) 
Jo 

^d{n-Ay 



+ ig ds Yfi{z;cc,s) 

10 



ds 



{z + ns)Yn{z;s,0). 



Using the chain rule. 



/ ds — [Yn{z; oo, s) [n ■ A]{z + ns)Yn{z; s, 
Jo ds 



ds 



— yn(z;oo,s) 
ds 



n ■ A]{z + ns)Yn{z; s, 0) 
d(n • A) 



ds 



{z + ns)Yn{z;s,0) 
" d 



poo 

+ / ds Yn{z;oo,s) 
Jo 

f'OO 

+ / ds Yn{z;oo, s) [n ■ A] {z + ns) 
Jo 

/•OO 

■ —ig / ds Yn{z;oo,s)[n ■ A{z + ns),n ■ A{z + ns)]Yn{z;s,0) 
Jo 

d{n-Ay 



ds 



Yn{z;s,0) 



ds Yn{z; oo, s) 



ds 



{z + ns)Yn{z;s,0). 



Using the above equation to eliminate the last term in Eq. I|5.1U3|) yields, 



(5.104) 



n-DYn{z) = ig ds Yn{z; oo, s)n''n''Gf,u{z + ns)Yn{z; s,0), (5.105) 
Jo 



where 
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n ■ DYn{z) = n ■ dYn{z) - ign ■ A{oo)Yn{z) + igYn{z)n ■ A{z), 



and the gluon field strength tensor is defined by, 



Hence 



Oi = iTr[yJ(m .D)Yn] = \T, 



ig I ds Y^{z;s,0)n^n''G^uYn{z;s,0) 





which is Eq. 



154 



Chapter 6 

Summary and Outlook 

"There is no supersymmetry." 

Carlos Wagner, in the middle of his 2002 TASI Lecture, 
Introduction to Supersymmetry. 

Experiments to search with unprecedented sensitivity for particles and phenomena beyond 
the Standard Model of particle physics lie just over the horizon. These experimental ad- 
vances make urgent the task of theorists to prepare to extract evidence for new physics 
from the empirical data by improving our ability to calculate reliably within the Standard 
Model, and then to be ready to deduce the implications of the new discoveries for our 
existing models of physical phenomena. 

In this work we have chipped modestly away at these complementary tasks. First, taking 
note of the already-existing evidence of the inadequacy of the Standard Model in the baryon 
asymmetry of the Universe, we took steps toward a more reliable calculation of the BAU 
predicted by the scenario of electroweak baryogenesis in the MSSM, and comparing the sizes 
of CP- violating phases required by this scenario with the limits already placed on them by 
the null results of searches for electric dipole moments of electrons, neutrons, and ^^^Hg 
atoms. Making use of the closed-time path formalism for quantum field theory incorporating 
the effects of finite-temperature and nonequilbrium physics, we found regions of parameter 
space in the MSSM for which the baryon asymmetry would be enhanced over the predictions 
of semi-classical calculations. In these regions, electroweak baryogenesis is found to account 
successfully for the BAU with CP- violating phases in the MSSM as small as 10~^~^, which 
are still consistent with current bounds from EDM searches. Future electron and neutron 
EDM searches, however, should be sensitive enough to rule out this scenario for baryogenesis 
if no EDMs (or too small EDMs) are found, or to detect positively nonzero EDMs of this size. 
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Such a discovery, combined with discoveries of superpartners with masses consistent with 
our predictions, would constitute strong evidence that the baryon asymmetry could indeed 
have been produced during the electroweak phase transition. More likely, EDM searches 
and LHC results will further restrict the parameter space in the MSSM that could allow for 
successful electroweak baryogenesis, and point our investigations in the direction of other 
models (either extending the MSSM or favoring a different mechanism for baryogenesis) 
that could be the correct account of the particles which actually exist in Nature and the 
origin of those which ended up as the matter that makes up all of our observable (and 
improbable |16nj ) Universe today. 

On the Standard Model front, we employed effective field theories in QCD to untangle 
perturbatively-calculable and unknown nonperturbative effects in processes involving the 
strong interactions by organizing the theory in powers of some small parameter. We applied 
soft-collinear effective theory combined with non-relativistic QCD to the exclusive radiative 
decay of T or other quarkonia to light hadrons. It has been known that in inclusive radiative 
T decays the decay channel involving the bb pair in a color-octet configuration can be as 
important as the color-singlet channel in the kinematic region of near-maximal photon 
energy, ~ My/2. In exclusive decays, however, we showed that the color-octet channel 
could be safely ignored, thanks to the power counting in the effective theory of the color- 
singlet and color-octet contributions to the T — > decay rate. The power counting and 
perturbative matching produced an operator which dominantly contributes to the decay 
rate, leading us to predict that the dominant decay product in radiative T decay should be 
the /2(1270), which should be produced in a helicity-zero state. The limited experimental 
evidence so far lends support to this prediction, although the helicity of the produced /2 
has not yet been positively identified. Our analysis also suggests a larger branching fraction 
for T — > 7/2 than the data so far indicate. Further data on this decay and other radiative 
decays, if consistent with the effective theory predictions, would lend powerful support to 
the validity of the effective theory expansion, establishing its reliability for use in other 
processes, especially those which relate more directly to searches for new particles. 

Such processes are the production of hadronic jets in lepton collisions or, in our formu- 
lation, Z decays. SCET isolated for us the soft gluon matrix elements giving the leading 
nonperturbative corrections to event shape distributions in hadronic Z decays or, equiv- 
alently, e~^e~ annhilation. Reliable calculation of these variables in QCD is essential to 
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separate the Standard Model background in such processes in searches for new physics. 
Universahties among nonperturbative contributions to different event shapes would reduce 
the uncertainties in these calcualtions greatly, improving our ability to find deviations from 
Standard Model predictions while also revealing new information about QCD itself. We 
attempted to find such relations in the Z decay distributions in jet energy, thrust, jet 
masses, jet broadenings, C parameter, and other variables. The theory, however, did not 
fully acquiesce to our hopes, leaving us only one relation between the thrust and jet mass 
sum distributions. Experimental tests of models proposing more extensive relations among 
these nonperturbative corrections seem to confirm our findings. This suggests a future di- 
rection of research to define new event shape variables, other than those standardly used in 
the past, which may receive universal nonperturbative corrections. This task was begun in 
perturbative QCD in Refs. 116111162] . and could be analyzed also in SCET. Similar methods 
could also be applied to other classes of events with QCD-jet backgrounds, or to the study 
of processes with collinear hadrons in the initial state, as occurs at any hadron collider. 

Precision tests at low energy which constrain parameters of proposed models of physics 
beyond the Standard Model together with reduction of uncertainties in calculations of strong 
interaction phenomena prepare the way for the further scrutiny of the Standard Model to 
ever higher precision and perhaps the discovery of new particles and phenomena which lie 
beyond it. It is the task of the theorist to develop reliable methods to predict the observables 
in these phenomena and be prepared to understand the implications thereof on our models 
for the particles which make up our Universe and how they were generated in the very 
beginning. 
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